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ANNOTATION /2* 


The articles of this collection encompass a broad 
range of questions associated with theoretical analysis 
and design of equipment used in conducting space ex- 
periments. The information on theoretical analysis 
of the possible apparatus solutions used to generate 
information streams aboard spacecraft is covered most 
completely. The articles on coding methods reflect 
the urgent necessity for more sophisticated onboard 
processing of the information obtained. The problems 
of spacecraft antenna testing, radiometric equipment, 
and so on are examined. 


The volume will be of interest to specialists 


connected with the design and construction of radio- 
electronic and radiophysical space equipment. 


nF i 8 


TABLE OF CONTENTS 


Page 
Series-Length Coding under a Pricri Uncertainty 1 
Y. M. Shtar'kov and V. F. Babkin 
Coding cf Discrete Monotonic Functions 14 
A. B. Kryukov 
Simple Method of Jumbering Binary Sequences with Given 24 
Humbem of Units 
Yu. M. Shtartkov and V. F. Babkin 
Multipurpose Information Collection and Processing 34 


Systems 
hb, V. Kantor, S. M. Perevertkin and T. 5S. Sheherbakova 


Associative Compressed Information Output Stream 43 
Formation by the Statistical Trial Method 
A. V. Kanten, T. A. Tolmadzheva 


Analytic Study of Output Stream Formation Process in 57 
Multipurpose Information Compression Systems 
A. V. Kantor, S. M. Perevertkin and T. S. Shcherbakova 


Dispersion Space Radio Links 68 
L. G. Sapogin and V. G. Sapogin 


Spacecraft Antenna System Design 93 
A. P. Alekseyev, B. A. Prigoda and L. I. Skotnikov 


Low-Silhouette Spacecraft Antenna Systems 108 
B. A. Prigoda 


High-Sensitivity 3.5-cm Modulation-Type Radiometer 1°5 
A, Ye. Andriyevskiy, A. G. Gorshkov, V. V. Danilov, 
Vv. K. Konnikova, A. S. Lobarev, Vv. G. Mirovskly, 
V. V. Nikitin, V. I. Portman, Ye. A. Spangenberg, 
I. A. Strukov, N. Z. Shvarts and V. S. Yetkin 


If Amplifier Liniting F-.1uency Selection in Super- 132 
heterodyne MM- and CM-Band Radiometer 


Yu. A. Nemlikher, I. A. Strukov and L. H. Yudina 


iv 


Page 


Study of Schotti:ky Barrier Diode Frequency Converter 142 
in the Short Millimeter Wavelength Band 

V. F. Kolomeytsev, Yu. Yu. Kulikov, A. M. Kupriyanov, 

I. A. Strakov, L. I. Fedoseyev, Yu. B. Khapin and 

Vv. S. Yetkin 


Influence of “hase Shifter on Frequency 148 
Divider Characteristics 
Ya. E. Veyber 


Comb-Line Bandpass Filters 167 
Ye. A. Vlasov 


SERIES-LENGTH CODING UNDER A PRIORI UNCERTAINTY 


Y. M. Shtartkov and V. F. Babkin 


ABSTRACT. A modification of the method of 
series length coding of bounded length is proposed. 
It 1s shown that if the binary symbol occurrence 
probabilities differ markedly, the proposed method may 
improve the message compression coefficient. The 
improvement depends on the degree of difference of 
the probabilities. 


INTRODUCTION 


N 
* 


The problem of the most effective statistical coding of a 
sequence of independent random events with known probabilities of 
their occurrence was solved long ago [1]. However, practical 
application of the methods developed has been hindered by the 
necessity for cumbersome "coding" and "decoding" tables and also 
absence of exact knowledge of the message individual probabilities. 


If the source entropy H < 1, the first obstacle can be over- 
come with the aid of the statistical coding technique known as 
series-length coding. The first mention of this method is en- 
countered in [2], a considerably more detailed analysis was made 
in [3], and a modification which is very convenient for practical 
purposes was proposed in [4, 5]. Several additional theoretical 
and experimental studies [6, 7] were made later. The attention 
devoted to this method is explained by the fact that this is a 
nonmnemonic technique and therefore there is no need to store 
coding and decoding tables. Moreover, in many cases the results 
obtained differ very little from the optimal results. 


Considerably more serious is the problem of a priori un- 
certainty. Only a few studies have discussed particular statis- 


* Numbers in margin indicate pagination in original foreign text. 
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tical coding methods with unknown individual message occurrence 
probabilities. These methods include, first of all, substitutions 
of the universal type [8]. A quite general approach to the 
solution of this problem was proposed quite recently, which 
unfortunately did not permit obtaining a coding technique 
acceptable from the viewpoint of complexity. The latter problem 
has been solved quite recently [9]. However, for many applications 
the proposed coding method is still not sufficiently simple. 
Therefore,in the present paper we reexamine series-length coding 
and propose modifications which are adapted for operation under 
conditions of partial or complete a priori uncertainty. 


UNIFORM SERIES-LENGTH CODING OF FINITE LENGTH 


Without significant loss of generality we limit ourselves to 


examining the sequence of symbols 0 and 1, which are statis- 
tically independent and have occurrence probabilities equal to 
q and p = 1+ q, respectively. We also assume “hat q > p. 


In the original variant [2, 3] series-length coding con- 
sisted in breaking down the subject binary sequence into individual 
blocks ("enlarged symbc's") of the form 


01 (1) 

Ces 

Owl 
and coding of each such block by a nonuniform code. However, fA 
all the techniques for such transformation known at the present 
time are excessively complex. In [3] a ternary code is used which 
Js by no means always convenient, while in [6] the use of codes 
with separator symbol is proposed, which is in p. inciple nearly 
equivalent to the ternary code but leads to less simplicity of 
code combination formation. Moreover, examination of the infinite 
sequence of enlarged symbols (1) is impossible in practice. In 
order to find the next enlarged symbol it is sufficient to count 


the number of zeros preceding the one. SBut any real counter can 
count up to a definite limit. 


A different solution was proposed in [4, 5}. In place of 
(1) the system of enlarged symbols was examined 


1 tig 


(2) 


We 


where each succeeding block contains one more symbol 0 than the 
preceding block, the next-to-last and last blocks contain 7 
symbols each, and the total number of blocks is thus equal to 

27 + 1. Considering that (2) is a complete system of events, such 
partitioning is always possibile. 


At the same time,changeover to the finite system of enlarged 
symbols makes it possible to code using blocks of the same length, 
which is impossible in principle when using the system (1), 
which contains an infinite number of elements. Let the number r 
satisfy the condition 


Then we can establish one-to-one correspondence between events 
from (2) and 2+ 1 (from - .:/ .) by binary sequences of length 
r. It is simplest to match with any y the binary form of the 
number 1 (using all r places). Then when ccding it is sufficient 
to count in an r=place counter the number of zeros prior to the 
occurrence of the first one (this rule changes somewhat for Ay). 


For any fixed r it is well to select the group of events (2) 
as large as possible. Therefore we take 7 such that 


eS, (3) 


The larger the ratio of the input block average length n to 
the output block average length m, the higher is the effectiveness 
of statistical coding. In the present case the output tvlock 
length is constant and equal to r. Ana it is not difficult to 
estimate the quantity n, considering that the symbols 0 and 1 
are statistically independent and the probabiisities of the events 
Ags Aj; Ay» faring Aji? Ay are equal, respectively, tO P, GP; a‘p, 


t-1 Zz 
oeey GQ Ps a7: 


ip 2ap-- 3qtp— ee eight tp ety St on (4) 
Then 


KK: Ke A 1.-q-""! (5) 
ier car Sr 7) or) a 2 


As is frequently done, we shall call the ratio K = n/m the com- 
pression coefficient. Thus, K = K (q, r) and K = K (p, r). In 
she figure,the dashed curves show the dependence of K on r for 
p= 107+, 107°, 1079 and 1074, respectively. Each curve has a 
maximum,and the value of the variable r for which the maximum is 
reached depends on p. The upper bound of the compression co~ 
efficient is 


% ( 
KP) ara (6) 


where if (p) -- ploggp — (l= pila, = ») 4g the source entropy. In the 
figure, the values of K* (p) for selected values of the probability 
p are shown by the horizontal aash-dot lines. We see that 

max Ap. differs very little from K* (p) and the smaller p,the 


smaller this difference. Thus, for known p and q = 1 - p we can 
select a value of r such that the effectiveness of uniform 
series-length coding of finite length will differ very little 
from optimal statistical coding. This conclusion is valid for 


values of p < 1071, However, if p > 107! then even K*® (p) < 2 


£5 


and the advisability of any statistical coding becomes questionable. 


The fact that for dif- 
RIED para LTTE CREE ferent values of p the curves 







in the figure reach a maximum 


200) for different values of r = 
Ca ge —- r (p) means that the subject 
” . method cannot be optimal for 
0 Of a eee a wide range of p values. In 
Ji selecting a definite value of 
ko fe a= ae a a r when constructing a coding 
| rs a scheme we must remember that 
Effectiveness of uniform and non- "he value of r which is optimal 
uniform coding for series: for one p may be far from 
l- p=10°'; 2p 2.107; 3- optimal for another. 
p=10%; 4-p- 10" 


Generally speaking, there 
is nothing surprising in this. 

None of the statistical coding methods examined in [1, 6, 7] 
permit constructing a code which provides nearly optimal effective- 
ness for various values cf p. Nor is uniform series-length 
coding of bounded length such a method. Only in [9-11] was tne 
possibility shown and constructive techniques presented for con- 
structing codes which are equally effective for various values 
of p. They differ only in slower approach to K*® (p) with increase 
of n. 


From this viewpoint it is interesting to recall the results 
obtained with nonuniform ccding of an infinite system of enlarged 
symbols (1). According to [2], when using a quaternary alphabet, 
which 1s equivalent to the use in binary code of a separating 
symbol consisting of two symbols [6], for any value of p in the 
interval from 0 to 0.5 the values of K (p) satisfy the condition 
0.65 K* (p) < K (p) s 0.8 K* (p). Thus, at the cost of failure 
of K (p) to approach K*® (p),a finite relative difference between 
K (p) and K*® (p) 1s provided. The longer the length of the se- 
parating symbol (or the number of symbols in the code used), the 
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closer K (p) will be to K* (p) as p + 0, but at the same time 
the smaller the value of K (p) as p + 1/2, 


For small p, when K* (p) »1, the difference between K (fp) 
and K* (p) of a factor of 1,2-1.5 1s quite acceptable in practice. 
However for p = 1/2 such coding would lead to an error of a factor 
of 1.5 in comparison with the pasic binary sequence, which is 
completely unacceptable. It is obvious that with broadening of 
the range of p values of 1.0 the results deteriorate still more. /6 
Finally, as we mentioned in the beginning o.' the section, examina- 
tion of an infinite system of enlarged symbols, even as simple 
as (1), can be of theoretical interest only. 


Let us return to examination of uniform series-length coding 
of finite length. Under conditions of partial a priori un- 
certainty we can use two different approaches to selection of the 
parameter r. Let the range of possible p values be given Py < 
P £ Po. In many cases it is important to ensure the largest 
value of the compression coefficient K in the worst possibile case. 


We see from the figure that for any r the compression coef- 
ficient takes the smallest value for p = Pos Therefore we should 


select r = Pf, satisfying the condition A (pay 7) = er NU r), which 

ensures a value of Kn PW (age Be max A (py + For example, let 
r 

ps 107°. In this case, according to Figure 1,7 = 8 and K (1077, 


8) = 11.5. For p= 1973 and r = 8 we obtain a larger value of K, 
equal to 28.2; howevex,this value row differs markedly from K* 
(1073) = 87.6 and the smaller p, the larger will be K (p, r) and 
the larger its difference fron K* (p). These examples show thst 
uniferm series-length coding of hounded length can be used under 
conditions of partial a priori uncertainty. Naturally, the result 
obtained will be better, the smaller the a priori uncertainty. 


SIMPLE METHOD OF NONUNIFORM SERIES-LENGTH .ODING OF BOUNDED LENGTH 


We mentioned above that for known p wa can select r such 
that uniform series-length coding (in the following we shall omit 
the words "of bounded length" for brevity, since the system (1) 
will not be examined further) yields a value of the compression 
coefficient which differs very little from K* (p). It appears 
that this was the reason why nonuniform series-length coding was 
not examined anywhere and was not mentioned. However, under a 
priori uncertainty conditions nenuniform coding is a natural 
technique which makes it possible to improve the results for a 
given range of p values. We emphasize immediately that with 
transition to nonuniform coding the extremely important property 
of simplicity, which makes it possible to avoid the use of special 
coding and decoding tables, must be retained. 


In order that the nonuniform code improve the re.ults for 
some range of p values,it is necessary that the K (p, r) curves 
in the figure have "flatter" maxima. It is this feature “Unat 
made it possible to select r such that K (p, r) for various p 
differed very little from max K (p, r) = K* (p). Naturally, 
this broadening of the maximum region must not (or may only to 
a slight degree) be accompanied by reduction of the magnitude 
of the maximum itself. 


In this sense the problem becomes quite nontrivial. The 
only thing that can be said about conventional statistical coding 
of the system of enlarged symbols (2) for fixed p is that it 
brings max K (p, r) close to K* (p) precisely for this p. But 

r 


this effect does not play any significant role, since max K (p, r) 
r 
* K*® (p). As for the broadening of the maximum region, it can- 


not be completely guaranteed. Moreover, the suspicions that 
nonuniform statistical series-length coding for given p impairs 
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the characteristics of this nonuniform code in relation to other 

p seem completely Justified. In this sense uniform codins be- /7 
haves less selectively and from our viewpoint hetter than non- 
uniform coding for piven p. Consequently the problem is to find 
those factors which influence the quite rapid falloff of the 

curves on both sides of their maxima, which are characteristic 

for different values of p, and which can be eliminated immediately 
for all p. In order to solve this problem we need to understand 

and clarify the behavior of the curves in the figure. 


Let us examine the maximum of the K (p, r) curve for a fixed 
value of p. Since it differs very little from K* (p), .-iform 
coding is close to optimal statistical coding. FPut this is valid 
if and only if the probabilities of the enlarged symbols Ay» A 
ess A, do not differ markedly from one another. 


1? 


Now let us begin to increase r (moving, along the curve to 
the right). In this case Z increases and the occurrence probability 
of A, equal to q’, A,_1 equal to qt} p, and so on decreases. 
At the same time the occurrence probabilities of Lo» Ay; Ay» and 
so on do not change. As a result,increase of r leads to ever 
greater and greater difference between the occurrence probabili- 
ties of the individual enlarged symbols and uniform coding becomes 
ever less and less optimal. This then expleins the slow but 


steady decline of the K (p, r) curve with increase of r. 


Thus, in order to eliminate or, more precisely, to weaken 
the decline of the K (p, r) curve with increase of r it is neces- 
sary to use a nonuniform code for which the length o. the code 
combination corresponding to Ay is a monotonically increasing in 
the broad sense function of 1. This is characteristic for any p, 
however the concrete form of such a function will depend on p. 

It is difficult to say anything definite about the results of 
such a transformation from the viewpoint of broadening the range 
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of p values. The only obvious effect is a shift of the maximun 
of each curve to the right. Moreover it is difficult to expect 
retention of coding and decoding simplicity. 


Now let us examine the behavior of the curve to the left of 
the maximum. ‘With decrease of r the occurrence probability of 
the enlarged symbol A, approaches q, and all the remaining Ay 
(i < 2) have practically the same probability. In this case 
we should code all Ay (1 < 1) by words of the same length and A, 
by e shorter word. The simplest solution lies in coding the 
enlarged symbol A, by a single binary symbol (0, for example) 
and any other symbol by r + 1 binary symbols, of which the first 
is equal to 1 and the remaining r coincide with tne previously 
used words of the uniform coee of length r. In this case the 
total number of enlarged symbols can be increasei from 27 + 1 = 2 
to 2° + J, which is equivalent to transition frcm 2 to Z' = 
Z* 122", Im this case the coding and decoding algorithms are 
just as simple as before. When codins, the combinations of 
z' = 2 + 1 zeros are replaced by a single zero and the series of 
i < 2' zeros (terminating in a one) are replaced by a one with 
ensuing r-place binary form of the number “is f—1=2’, 
Decoding procecds just as simply. 


r 


We could rep .ace A, by a sequence of Po =2, 350 sory Pm 
symbols. When using lengths r + 1 for all the other code com- 
binations this would permit increasing the total number of en- 
larged symbols to “"* -4 +«!--1. This increase is always useful 
to some degree. On the other hand, however, for small values of 
r codin,;; becomes less optimal than with the choice ly = 1. And 
for large r, when the A, occurrence probability becomes smail, 


selection of the quantity r, does not play any role. We can 


f" 


obviously state that the proposed nonuniform series-length 
coding variant leads to shift of tne K (p, r) curve maxima to 
the left, towards smaller r. In fact, for uniform coding the 


maximum is reacned in the case when the Ay occurrence probability 
is of the same order as that of Ay (zp). But for the proposed 
nonuniform coding the maximum is reached with Ate] occurrence 
probabiiity equal to approximately 0.5. For p < 0.5 the second 
condition is always satisfied for smaller values of r than the 
first. This then means shift of the maximum in the direction of 
smaller r. 


Is 
@& 


wow let us examine quantitatively the behavior of the com- 
pression coefficient K (p, r) for nonuniform series-length coding. 
For greater convenience of comparison we immediately agree to 
take K (p, r) to mean the compression coefficient for the case 
when 2' = 2” and the length of all the code words other than 
that corresponding to Atel is equal to r+i1. This makes it 
possible to compare uniform and nonuniform coding with practically 
the same total number of enlarged symbols in (2) and the same 
maximal series length. Thus, K (p, r) = f/m. The quantity ni 
can be found from (4) with replacement of z7 by 2' =2+1 = ar. 
On the other hand, m#r. It is not difficult to verify that 


Me = (M— git ar) git (4 — gitt) — 1, (7) 
Now, using (4), (5), (7), we obtain 


~~ q—ett 
A (p,7r) =: a ig) = — 


= _ (8) 


(a =~ 49) r-; a 
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If we reglect the difference between 1 - q? and 1] = q’*t 
we obtain the simple approximate formula 
" A ri 
A(pyr)= fs 2 
near se aa (9) 


r 


The soiid curves in Figure 1 show the behavior of K (p, r) as a 
function of r for the same probability vaJues as before. 


iu 


Compariscn of the K (p, r) and X (p, r) curves for the same value 
of p shows significant broadening of the K (p, r) maximum region 
at the expense of marked increase of the compression coefficient 
for small r. Thus the posed problem is solved at least to some 
degree. 


We note that, aS expected, the K (p, r) curve maxima are 
shifted to the left relative to the K (p, r) maxima. This means 
that for any problem formulation the best value of r is less 
than for uniform coding. This leads to decrease of the maximal 
series-length, which may be very useful in the presence of indi- 
vidual nonstationary segments in the sequence. It is very in- 
teresting that max K (p, r) > max K (p, r) for any p. In other 

r r 


words, broadening of the maxima region is accompanied not by 
reduction but rather even by some increase of the maximum itself. 
Finally, we note that in the region of large r the K (p, r) 
values are somewhat lower than K (p, r), but this difference is 
not significant. In accordance with (10) this occurs for suf- 
ficiently large r, when agt - However, the deterioration 


introduced in this case does not exceed 1/r and is small in 
comparison with the one in the derominator. Thus the proposed 
method makes it possible to obtain immediately several advantages. 
Over the entire range of r values the compression coefficient 

K (p, r) either exceeds or practically coincides with K (p, r). 
The width of the K (p, r) maximum region is greater than for 

K (p, r) and the magnitude of the maximum itself is larger. The 
maximum is reached for a smaller value of the parameter r and 
correspondingly smaller maximal series-length. 


In conclusion we shall illustrate the noted advantages by /9 
examples. To this end we use the examples considered previously 
for K (p, r) and then compare the results obtained. For example, 
ifovsc'* lu- and we are required to maximize the value of K 
in the worst poss!ble case, then in accordance with the figure we 
11 


find that we should take r = 6 and this provides a value of the 
compression coefficient no worse than K = 12.1. This result 

nearly coincides with that obtained previously, and this is 

quite natural, since the present problem reduces to determine 

max K (p*, r), which as we noted previously differs very little 
r 


from max K (p, r). However, for p* = 1073 kK (1073, 6) = 44.4, 
i.e., it is considerably larger than for the selection (which is 
optimal in this sense)of the parameter in the uniform coding case 
‘K (1073, 8) = 28.2). The aifference will be somewhat greater 
for p = 107", K qo", 6) = 60.7 > K aon, 8) = 31.5. The 
examples show that in those cases when it is known in advance 

that the probability p will not exceed some quite small magnitude, 
nonuniform stress-length coding is better than uniform. However, 
if in reality the probability p is noticeably less than the indi- 
cated limit this a_fference may be very significant. 


We shall list some factors which can influence the choice 
between uniform and nonuniform serics-length coding under con- 
ditions of partial a priori uncertainty. First,the optimal 
series-lengths are different in the two cases. The optimal 
series-length is shorter for nonuniform coding, and therefore 
returning of the coder with the appearance of nonstationarities 
in the process being coded may be simplified considerably. 
Second, with optimal choice of the length of the series being 
coded uniform and nonuniform coding provide approximately the 
same gain, quite close to the limiting possible gain. In those 
cases when the a priori estimate of the occurrence probability 
of any symbo] (0 or 1) is not sufficiently exact, and in reality 
this probability may ¢ considerably less, nonuniform series- 
length coding may provide considerable improvement, which will 
be greater the inore the true probability differs from its 
estimate. 


le 
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CODING OF DISCRETE MONOTONIC FUNCTIONS 


A. B. Kryukov 


ABSTRACT. Two techniques are proposed for 
numbering all monotone functions which permit 
economic coding of any such functions by a uniform 
code. We analyze the effectiveness of this function 
representation and, in particular, it is shown 
that for a fixed number of quantization levels the 
compression coefficient increases without limit 
with increase of the number of function veadings. 

A characteristic feature of the techniques examined 
is the possibility of forming the code word in 
real time as the coded function readings arrive. 


INTRODUCTION 


Ns 


We examine a set of discrete functions, each of which is 
specified by a sequence of readings {s,, Sy, °°"? Spa sess S),}3 
where the readings 8, can take any integral values from 1 to m. 

The number of all such possible functions is equal to m; there- 
fore, for coding any of them by a uniform code, we need to expend 
‘ant, we‘) Q=nary units. (Here and hereafter <&> denotes the smallest 
integer larger than or equal to €.) In certain cases we can 

impase qualitative limitations on the nature of the function 
behavior. Then the number of possible functions reduces and, 
consequently, such a representation of the functions becomes 
redundant. For example, in a telemetry system with cyclic sampling 
of the sensors, as a result of the use of compression there is 
formed a sequence of significant readings whose addresses (in 

the limits of a single cycle) form a strictly increasing sequence. 
In [1] an effective technique for coding such address sequences 

is proposed, based on numbering all the strictly increasing 
functions. 


In the present study, we propose two techniques for numbering 
functions from a broader class -— functions which are monotonic 
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in the broad sense. An example of a case when such a limita- 
tion is imposed on the functions is the adaptive telemetry 
system in which the same sensor may be sampled several times in 
succession. In this case the significant reading addresses form 
a nondecreasing sequence. The techniques proposed below make it 
possible to encode economically by a uniform code nondecreasing 
or nonincreasing discrete functions and also both together. In 
the latter case it is necessary to expend one additional bit to 
indicate the sign of the increment (nondecrease or nonincrease). 


In the following, for definiteness, we shall examine only 
nondecreasing functions s5,;S%: ... <8: ..e Sey, 


CODING TECHNIQUES 


Technique 1. In this case, coding is accomplished in two 
steps. In the first step, the function is coded by a binary 
sequence of length m + n - 1 and in the second step, the obtained 
sequence is recoded into a different, more compact binary 
sequence or number. 


We assign to each nondecreasing function a sequence of 
zeros and ones By> Apo sees Anna This sequence consists of 
two parts. The first n-1 sequence elements are formed using 


the rule fl, 4f sa, 
a= {0 if ee fi=m1,2,...,r—1, 
’ Mea Ob 
and the remaining m elements are formed using the rule 


1,1f jes," 
Qy1a5 = prod,2,...,m 
Nwle? a 4f jEaS, i 7 ’ ® 
where S is the set of values which the function being coded takes. 


The formation of this sequence is shown in Figure 1. The /1l 
first n=-1 sequence elements are written in row form, so that the 
i-th element of this part of the sequence corresponds to the 


{(1+1)-th function reading. The remaining m sequence elements 





* Translator's note: There is obviously a mistak2 in the equa%ion 
in the foreign text. 
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are written in column form, 


fc eater bavi so that the (n-1+j)-th element 
aad ae ae Emre gece gee corresponds to the j-th quanti- 

iter a ee Ore ke zation levei. The order of 

5 OOO HQ te 

Si Olen Sah lies cose an sequence element numbering in 
een Sere ema ig the figure is shown by the arrow. 
d é a 7) 





yoQr??hte? pf) i In the row, the symbol 1 denotes 


those readings which do not 
Figure 1. Formation of binary . 
sequence of length mtn-1 (m= differ from the preceding, 


8, n=11) while in the column the symbol 

1 denotes the values which the 
function takes. Thus, to each nondecreasing function there 
corresponds a binary sequence of length mtn-1, and to the dif- 
ferent functions there correspond different sequences which 
differ in the location of the ones. We note that a characteristic 
property of these sequences is the fact that the number of ones 
in the sequence is always constant and equal to n; the number of 
zeros is correspondingly equal to m-1. In fact, to the first 
reading there always corresponds the symbol 1 in the column and 
each following reading either does not differ from the preceding 
{then a symbol 1 in the row corresponds to this reading) or it 
takes a new value (then a symbol 1 in the column corresponds to 
it). Proof of one-to-one correspondence between the functions 
being coded and the sequences may be based on the fact that the 
positions of the zeros in the row and ones in the column are 
simply the coordinates of the characteristic points (jumps) of 
the function being coded. 


In the second step recoding of the obtained sequence into 
a shorter binary sequence (or number) is accomplished. The 
number of binary sequences of length mtn-1 containing n ones is 


. They can be recoded by a special technique [1] 


equal to Catn-1° 


so that <log, Cnea+. Q-nary untts will be expended on each. To 
: , 
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this end, each sequence is assigned a definite number A, where 
OsACC as 1. For formation of the numbers A, we use the 


numbers of those sequence positions where ones (or zeros) appear. 


We note that, generally speaking, the numbers A can be formed by 
both ones and zeros of the sequence, since Ch-n-z: Cnina, but 

the number of arithmetic operations is proportional to the 
number of those symbols which are usei for this purpose. Thus 
the number of arithmetic operations is proportional to min (n, 
m-1). 


Technique 2. While in the first technique each function 
was initially assigned a binary sequence of length m+n-l and 
then this sequence was coded, in the second technique direct 
coding of the functions is accomplished, bypassing the binary 
sequence formation stage. To this end we number the functions 
being coded, i.e., each of them is assigned a definite number B, 
where O<B<Ciy,1-—1. Naturally, the numbers B corresponding 
to different functions are different. The rule for forming 
these numbers is as follows. 


As before, we assume the function is given by the sequence 
of readings {S15 Spo sees Seo cers st. Then the number B for 
this function is defined as the sum 


a ee (1) 
ea 


where ef is assumed equal to zero if B>a. It is not difficult 
to show with the aid of relation (2 (3,5 3 weet CL, that 
4n accordance with (1) different numbers B correspond to the 
different functions. 


Let there be given two nondecreasing functions {s. s....., 8, 


os} And {s. s, .60) 5-8, and let the largest number of the 
reading where the functions differ be equal to 1, and all the 
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readings with numbers larger than 1 be the same for both functions, 
1.€., %4 Sin Ss Sigeeee. 8, oS. For definiteness, we take 
“,.-s, and show that in this case /’_. #” always. 

The terms in (1) corresponding to readings with numbers 
larger than i are the same for both functicns: we denote their 
sum by R. Since readings which are larger in magnitude correspond 
to larger terms in (1), the "worst" case will be that in which 
Spe eee peg ANA HP Se! Le esi.) «Then the first 1-1 


terms in (1) for the function {s}} are equal to zero and 


B ; c ae R. 


< tig 
For the function {s}} 


BY -¢c', ots : 4 ¢ ~R-= 


(jrimarni | (ajeteaaciety 7 Gaiateapen 


aC 1-:- -- Bi —1. 


ried 


Thus, B' > B" even in the worst case, and this inequality is 
satisfied for any 1 (1 = 1, 2, ..., n). 


The decoding rule amounts to the following. In the fact 
step, we seek that value of k (k = 1, 2, ..., m) for which the 
inequalities are satisfied 


Chena B < Chena (2) 


Let (2) be satisfied for k = sf. Then in the next step we seek 
that k (k ® 1, 2, ..., sf) for which the inequalities are satisfied 


Cinta se B tr << O Ene -a (3) 
Let the condition (3) be satisfied for the value k = s#_,. [In 
these inequalities, as in (1), or is assumed equal to 0 if A>a.] 
This procedure is continued until obtaining the sequence (sf, 88, 
ee sh}. We note that if at any stage of the decoding pr’cedure 


18 


—- for example, when seeking the value of the r-th reading —— 

the corresponding inequality becomes an equality for k = at, 

the procedure may be terminated, since all of having i<r are equal 
to 1. Uniqueness of decoding is easily proved with the aid of 

the relation ¢;: C2), Ch ocr... + Cisa,. By virtue of one-to-one 
correspondence between the numbers B and the functions heing 
coded, {sf} = {s,}. 


The coding and decoding procedures described above can be 
illustrated with the aid of a special table composed from the 





numbers Ci.j.2 (6-1, 20... np kee de ee mt £13 
mM Chek oe Cs 1 Gh. 2 see Cinencs 
Teo ASS righ eh as he Nain Sel gd ndings ei Sew Steer Mau lantectae eget Nive’ “ahaa aes 
Aceh. 2 i oc 1 [es 13) Chis 
eed ee pee eee Cin 
Syed CRs F GEA as ch 
t | fo! u 0 0 0 
cra ge 


With the exception of the bottom row, consisting of zeros only, 
this table is a Pascal triangle in which the extreme left column 
of numbers (j. (i, ... Ch. is missing. Like the Pascal triangle, 
this table can be constructed recursively with the aid of the 
relation C2 2 Cy,-- Cx). 


Let us examine coding and decoding with the aid of this table 
for an example. Let m= 4, n = 4 and the function be given by 
the readings {1, 3, 3, 4}. By “superposing" the function on the 
table, we find the number B as the sum of the table elements cor~ 
responding to the function values: B= 0+ 3+ 4 +415 © 22. (In 
the table, these elements are shown bracketed.) This result 
naturally coincides with the value of B calculated using (1). 
We see from the table that the minimal value of B is equal to 0 
for the function {1, 1, ..., 1}, while the maximal value for 
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{m, m, ..., m} 18 equal to 
it Cpa Cheese! Wana Cocky 22; 

Now let us perform the decoding. In the first step, we 
seek in the fourth column the largest number not exceeding Be22. 
This is the number 15 located in the fourth row. The row number 
is equal to the reading value, therefore Ss, * 4, We calculate 
the difference 22-15 = 7 and in the third column we seek the 
largest number not exceeding 7. This is the number 4, located 
in the third row, which means that 83 = 3. We then calculate 
the difference 7-4 = 3, seek in the second column the largest 
number not exceeding 3, and find in the third row of this column 
the number 3. This means that 85 3 3 and s,° 1, since 3-3 = 0. 
Thus, as a result of decoding we obtain the original reading 
sequence {1, 3, 3, 4}. 


We note that, when forming the numbers B, the number of 
arithmetic operations is determined only by the value of n. 
Therefore, when n > m- 1, a smaller number of arithmetic 
operations is required for realization of the first technique. 
However, at the same time cuding using the first technique is 
conducted in two steps, which creates additional difficulties 
not characteristic of the second technique. It is significant 
that both techniques permit coding of the functions directly 
at the rate of reading arrival, since formation of the numbers 
A and B can be accomplished as the readings arrive and terminates 
with arrival of the last reading. In this case, storage of the 
preceding reading values is not required. 


ESTIMATE OF EFFECTIVENESS 
The examined techniques for coding nondecreasing functions 
make it possible to expend (log, mia.» qenary units on each of 


them. In the conventional coding technique, this requires 
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n<log, 4) qenary units. For "> 2Cmnicim" , therefore the 

proposed coding techniques make it possible to reduce the number ZA4 
of symbols expended. The magnitude of the advantage obtained in 

this case, i.e., the degree of reduction of the number of q-nary 
units, is usually characterized by the compression coefficient 


ang dugg, ne > 


K eom™ clog CR a hm 


The magnitude of the compression coefficient (4) depends signifi- 
cantly on m andn. In order to clarify the nature of this de- 
pendence, we shall examine two cases. 


Acom 
24) 





2! 2 2F a7 29 ma 


Figure 2. Compression coefficient as function of n and m (q=2). 


l- m=32, n variable; 2= n#32, m variable (smooth curves are 
drawn through che calculated points for better visualization) 


Case 1. m fixed, n variable. Figure 2 shows the compression 
coefficients fer various n (n = a 2°, srs’ 27) for the particular 
case m = 32 (curve 1). We see that the compression coefficient 
increases monotonically with increase of n. As n*™», neglecting 
the roundoff symbol <*> and using the relations 


ln(L +e) se—t pms, 


we can transform (4) to the form 


alos pn 


~ | ii ez Pate 7 2 
K oom TED ae [t era be (5) 
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where O(e€) denotes a function of ¢ such that 


* Ow) 
aera const. 
~ g 


It follows from (5) that the compression coefficient increases 
without bound as n>, Thus, the larger n, the greater the ad- 
vantage obtained using the proposed coding techniques. 


Case 2. n fixed, m variable. This case is examined to 
clarify the dependence of K oom on the number of quantization 
levels. The values of the compression coefficient for various 
m (m = 21, 2° apes 29) and fixed n = 32 are also shown in 


Figure 2 (curve 2). For fixed n and m+~, the expression (4) with 


the aid of the same transformations as in the preceding case 
can be reduced to the form 


Kia ee (6) 


In this case, in contrast with the preceding case, the expres- 
sion for the compression coefficient has the limit 


lim Kin = il, 
m?*o 
n=const 


The difference in the results obtained for cases 1 and 2 
is explained as follows. The number of nondecreasing functions 
(and the number of symbols expended for nonredundant coding) 
depends on n and m practically identically. (These relations 
are absolutely identical if we examine m-1 rather than m as the 
argument.) At the same time, the basic representation of the 
functions, i.e., conventional coding with the aid of n <1log, m> 
q-nary units, does not have such symmetry — namely, the number 
of symbols expended increases far faster with increase of n 
than with increase of m. In other words, for n>m this form of 
basic monotone function representation has greater redundancy 
than for m>n. Therefore for m»n nonredundant coding of nondec- 
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creasing functions may not lead to any appreciable advantage in 
comparison with the conventional coding technique, since the 
redundancy of the latter is not large in this case. 


This last remark makes it possible to conclude that non- 
redundant coding (numbering) of histograms by a uniform code may 
not lead to any significant reduction of the number of symbols 
expended. In fact, we can construct a distribution function 
(integral law) from each histogram. Since the distribution 
functions do not decrease,one of the techniques proposed above 
can be us=d for their coding. However, this dees not lead to 
an; marked advantage, since the relation m»n is characteristic 
for histograms. 
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SIMPLE METHOD OF J:‘JMBERING BINARY SEQUENCES 
WITH GIVEN NUMBER OF UNITS 


Yu.M. Shtar'kov and V.F. Babkin 


ABSTRACT. A method for voding messages 
with a priori unknown occurrence probabilities 
is proposed which has a very simple realization. 
The calculations presented show that the method 
effectiveness is close to the limiting possible 
value — the effectiveness of universal coding. 
Various modifications of the method are discussed. 


INTRODUCTION 


In direct application of the statistical coding ideas, it 
is necessary to use cumbersome tables, whose complexity increases 
exponentially with increase of the length n of the sequences 
(blocks) being coded. However, in many cases coding can be 
reduced to uniform numbering of all blocks with given limitations 
(properties). This approach was first used in [1] to code an 
independent source with unknown statistics. In this case, the 
technique proposed in [2] was used for uniform numbering of all 
binary symbol blocks containing a fixed number t of ones. If 
15> 15; stay 1, are the positions of ones in a specific block, it 
is assigned the number 


Osa, = Cir Cla te Cy SOA, (1) 


where cf is considered equal to zero if iej-+-i1. This numbering 
ensures one-to-one correspondence of all the a and blocks of 
length n with t ones. 


The labor involved in the numbering (1) increases only as 
a power-law function of n. Nevertheless, a specialized device 
which computes or stores the binomial coefficient values would 
be quite complex. Therefore, the problem arises of finding 
numbering techniques whose effectiveness is nearly optimal and 
whose realization is considerably simpler. Since the uniform 
numbering problem arises only in application to equi-probable 


2h 


~N 
~~ 
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sequences, in the general case the effectivencss can be defined 
as the average number of symbols expended in describing a single 
number. It is precisely on the basis of this criterion that we 
evaluate the effectiveness of the simple numbering technique 
examined below. 


NUMBERING ALGORITHM 


The only simple technique for statistical coding of binary 
sequences is known as series-length coding (its description and 
a detailed bibliography are presented in [3])}. The basic se- 
quence is broken down into "enlarged symbols" A,» containing 
i zeros and a following one(v-"-i<f 2-1: or 2 = 2" ~ 1 zeros 
(A,). Then each A, is replaced by an r-digit number B, > whose 
magnitude is equal to i. For any value of the probability 
m+; Of one occuring, we can select the parameter r to ensure 
high coding effectiveness. This method is realized very simply 
in practice with the aid of an r-place binary counter. We use 
series-length coding to number the blocks containing t ones. 
Beginning with the first symbol of the block, we partition into 
enlarged symbols A, up to and including the last, t-th, one (this 
is always possible) and replace A, by B The obtained code 
words satisfy the following conditions: 

1) a single code word corresponds to each block; 

2) a unique block corresponds to each code werd; 

3) the ensemble of all code words has the pr2fix property [4]. 


i° 


The first assertion follows directly from the coding techni- 
que described above. For proof of the second assertion, it is 
sufficient to consider that each code word consists of r-place 
numbers B,- Therefore, breaking any code word into groups of r 
symbols each, we define the sequerce of numbers B,. Replacing 
each of them by the corresponding enlarged symbol Ay» we restore 
uniquely the initial symbols of the original block up to and 
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including the t-th one. Adding to the restored part of the 

block a series of zeros complementing the already obtained number 
of symbols to n, we obtain the original block, which proves the 
second assertion. 


Finally, the prefix property means that no code word is the 
beginning of any other code word. In the present case, this 
property is due to the fact that each code word contains t 
"numbers" By» 42 and «<(2— 2 numbers By» and the last number 
is not equal to By (coding terminates after the appearance of 
the t-th one and identification of the corresponding A, ¥ A). In 
fact, let us assume that the beginning of the code word coincides 
with another code word. Then, just as the second code word, this 
beginning must contain t numbers B, (i # 2). But since the last 
number of any code word is not equal to Bi» on the whole the first 
code word must contain no less than t + 1 numbers B, (1 #2), 
which contradicts the code word structure. Therefore, the be- 
ginning of any code word cannot coincide with anotLer code word, 
which proves the prefix property. 


In accordance with the first property, any block or sequence 
of blocks is uniquely transformed into a code word sequence. Con- 
versely, according to the third property, from any sequence of 
symbols we can identify uniquely the ensemble of any known number 
of code words if the position of the first symbol of the first 
code word is known. According to the second property, from the ALL 
identified code words we can recover the original sequence of 
blocks. Thus, there exists a one-to-one correspondence between 
the. sequence of blocks and the sequence of code words. This makes 
it possible to consider the code words as numbers of the cor- 
responding blocks. The fundamental difference between this 
numbering technique and (1) consists in the fact that the number 
"length" m, 1s not a constant quantity. In accordance with the 
criterion formulated above, the effectiveness is determined by 
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comparing logs e: (the minimal fixed number length) and the average 
number Mm, of symbols used to describe a single number when using 
the subject technique. Since all cr blocks with t ones are 

equally probable, then 


ig 22 —Z— DM, (2) 


where the summation extends over all such blocks. 


It is obvious that m = Mm, (r). Therefore, for each t we 
must use the coding parameter r = r, which minimizes ii, « Since 
series (of zeros) length coding is effective for p ¢ 1/2, for 
t > n/2 we obtain r, = 1 and mi, = n. In order to eliminate this 
drawback, for t > n/2 we will code not the original block but 
rather a "complementing" block, obtained by replacing the zeros 
by ones and the ones by zeros. It is obvious that then rm = 
ryt and Mm, = mt: Correspondingly, after decoding it is 
necessary to replace the obtained block by the complementary 
block. 


Thus, in order to select r, (0 < t <¢ n/2) and evaluate the 
effectiveness of the subject numbering technique we need to obtain 
the expression for mi, (r). 


AVERAGE CODE WORD LENGTH 


Let us determine the number t + 5, of enlarged symbols obtained 


when using the algorithm described above, where Si - at is the 


number of enlarged symbols Ay: If 1,» los eeeg 1, are the 
positions of the ones in a specific block and 


Ug = ig — bgey De aE ey Oxsiy<l, Kd eed (3) 


(4, is assumed to be equal to zero), then for the subject block 
0 


a ze ( 4) 
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In fact, it 1s not difficult to see that upon partitioning 
there will be obtained the following sequence of enlarged symbols 


AA, eee aA, A, A,A, wes aly al, oe AA; eee Ay Aw. 


2 
Pi oa , a 


And since m =r (i + S,)> using (4) we obtain 


¢ 
acnrwcajer gs Bet), (5) 
kal 


where the overbar denotes ensemble averaging. 


In accordance with (3) 


(ue) fa2 


Ps =y = J Pea =) (6) 


xeayl 


where P (36%) = y) and P (a, = x) are the protabilities that 


a{k) = y and Oy, = X. Then, considering that P (a, *x>n-t) = 


0, we obtain 


~ 
Ke 
fos) 


e 
mn <1 : 1 ok 
5” m. > P (si) wi yyy e >» P (% . 7 yy “a 
wd we 
O—1 (ued) I-21 


veo 
Le »{ >: P (ay = apy i) »y P (4 =<), 
xed 


veel x=ul 


n=—f!} 


where 0- |---| and [+i is the whole part of z. Thus the problem 
reduces to determining the distributions P (a). = x), k= 1, 2, 
«+e, t. Since all the blocks are considered to be equiprobable, 
P (Co, = x) is equal tothe number of positions in which a, * Xs 


divided by Cr. 
We first of all show that 


Pla, a) = Play cy be 2, 8 ee te (7) 


For this, we examine an arbitrary arrangement of t ones specified 
by the vector (4145 dns eres 4.) or, what is the same, (ay > Oo» 
s+, 4). If a, # a, there exists exactly one vector (41. 4... %) 
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in which af = a3 aft = a, = a, (j #1, k). But then the number 
of arrangements in which a, = x is exactly equal to the number of 
arrangements in which OH, = xX. In fact, any arrangement with 

a, > Oo = xX is considered simultaneously in P (a, = x) and P (a, 
x) and to any arrangement with a, = Xx # On, there corresponds one 
and only one arrangement with af = 4,5 af = x and a, = af (j3 #1, 
k). The equality of the number of arrangements with > 
and a), = X proves (7), and it is necessary to determine only 

P (a, = x). If the first one in the block is located at the 

i) = a, + 1l-th position, the remaining t - 1 ones may be arranged 
among the remaining n-1, symbols in Cx, = Chua ways. This 
expression defines the number of arrangements with the specified 
a, and P (a, =x) = Cy4.4/C,. Then in accordance with (6) 


= x 


-1 wt) let 
gl Eo #131 2 Crh by 4 ret ia ofS cry] }i.0 Q=!A5- ’ ( 8) 
and from (5) and (7) 
t 
my =r (t Da) = ree eat) a ret + 3"). (9) 


Keo 
Relations (8) and (9) define completely the average code word 
(block number) length. 


EFFECTIVENESS OF THE SUBJECT NUMBERING METHOD 


Theoretical analysis of (8) and (9) is not possible. MThere- 
fore, we calculated the quantities Mm, for n = 63 and t = 0.1, 
++, 31 (as noted above, i, = mh t)> Figure 1 shows the results 
of these calculations. Curve 1 shows the behavior of mf = mf = 
log, cr corresponding to optimal uniform numbering, and curve 2 
shows the behavior of my, when using the algorithm examined above. 
As we would expect, me < mi, for any t > l, and max {f, = ny} = fy — 
My == fiyy — My = Migg ~ My = 42 binary symbols, and max (ih, — mim, = 
(fi, — m)'img = 1,28, Thus, considerable simplification of the coding /19 


procedure is achieved at the cost of effectiveness deterioration 
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by no more than 4,2 binary symbols or a factor 1.28. In many 
cases such losses are quite acceptable. Thus, when coding a 
sourze with unknown statistics,the coding redundancy of a block 
of length n cannot be made less than 0.5 log, n[5], which in 

the present case constitutes 3.0 binary symbols and is comparable 
with the losses of this simple numbering technique. 


80: In applications such as 
wl 4 space studies, certain other 

' aos 

Pe an characteristics are also impor- 
199° 


tant — for example, the maximal 
possible gain in describing a 
block, equal to the difference 

of the maximal numver length 

and the number of symbols in the 
block. In the optimal uniform 
numbering case, this difference 
is negative for any t (its 
maximum is reached for t = 31 and 


Comparison of effectiveness of 32 and is equal to 3.3 binary 


universal coding method and = 
universal coding method with symbols). When using the pro 


the use of series-length coding posed numbering technique, m,. 1s 
woth Di6ck, deneth net also always smaller than n, 

but the maximal value of m, may 
exceed mi, considerably and for large t may even exceed n. It 
is not difficult to see that 





¢ 0 02 16 10 a4 tO ¢ 


MAX My 2 Ty tg : (10) 
miu my: ryt. 

The horizontal dash-dot straight line 3 in the figure corresponds 

to the value n = 57, and the dashed broken lines 4 and 5 cor- 

respond to the upper and lower bounds of (10). The maximal 

value of m - nN is reached for t = 23, and is equal to nine 
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binary symbols. This loss is not so large that the danger of 
obtaining a small group of “inconvenient" arrangement would 
force us to reject the use of this numbering technique. 


There 1s a simple way to limit max im, - n}. We write 
ahead of each number an additional binary symbol and stipulate 
that if it is a zero the following number is constructed using 
the technique described above, and if it is a one the following 
symbols are the direct form of the "numbered" block. With 
this modification max {m, - n} = 1 at the expense of increasing 
the length of most of the numbers by a single binary symbol. 


Finally, we shall discuss briefly possible improvements of 
the subject method which permit approaching Mh, to m¢# and reducing 
the dispersion of mM, The value of vr, is a function not only of 
the number t of ones but also of the block length n. For every 
n', where 1< n' <n, and every t', where 0 < t' < n', we can de- 
termine the optimal value of the parameter r : mnorP (n! 4 t")s 
Then the previously described numbering technique may be altered 
as follows. At the initial moment we use ae (n, t). Let the 
first enlarged symbol be A, (i < 2). This means that in the re- 
maining part of the block of length n - 1 = 1 symbols there are 
t - l ones. But for such a "shortened" block with t - 1 ones there 
exists the optimal value of ee 7 (n-1i- 1, t - 1), which may not 
coincide with r (n, t). Therefore, the second enlarged symbol is 
identified and coded using the parameter r (n-i-+1, t - 1). 
However, if the first enlarged symbol is equal to Ay» then in the 
second step of the numbering we use r (n - 7, t). After the /20 
second step we again correct the parameter r, and so on until the 
appearance of the t-th one, after which the numbering process 
terminates. It is not difficult to see that in this case the 
prefix property is again satisfied and uniqueness of the decoding 
is ensured. The r (n, t) used in the first step is known a priori. 
Therefore we initially identify the r (n, t) first binary symbols 
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and the enlarged symbol Ay. In the next step we identify the 
r(n-ie-1, t = 1) binary symbols if 1< 2 or r (n = 2, t) 
binary symbols if 1 @ 2. Using the identified binary symbols 
we determine the second enlarged symbol, the new value of r, 
and so on until the appearance of the t-th one. 


The use of this proceaure makes it possible to reduce My in 
comparison with the basic subject technique. This follows 
directly from the fact that in the basic subject numbering 
variant we select the value of ry, which is optimal in the mean 
for numbering the entire block (up to the t-th one), while in 
the proposed modification of the technique we select the optimal 
value of r in each step of the numbering process (after identify- 
ing the next enlarged symbol). The only drawback of this modifi- 
cation is the necessity for realizing the function r (n, t) of 
two variables. For r, with fixed n this problem is resolved 
simply — it is sufficient to have a log, (n + 1)-place counter 
with the same number of coincidence circuits and a log, logs 
(n + 1)=-place counter, since r. is a monotone function of 
t (t sn/2). In spite of the fact that r (n, t) is monotone 
with respect to each variable,its realization is incomparably 
more complex. 


We note that it is possible to construct other modifications 
of the numbering methcd, for example, using nonuniform series- 
length coding [3] or modification associated with application of 
the subject technique in the universal coding procedure. However, 
these variants will not be examined here. 


REFERENCES 
1. Babkin, V.F. In the collection: Problemy Peredachi In- 
formatsii. Vol. 7, No. 4, 1971, p. 13. 


2. Lynch, T. Proceedings of IEEE. Vol. 54, 1966, p. 1490. 


32 


3. Shtar'kov, Yu.M. and V.F. Babkin. See article in present 
volumn. 


4, Fano, R.M. Peredacha informatsi1. Statisticheskaya teorlya 
svyazi (Transmission of Information. Statistical Theory 
of Communication). Mir Press, Moscow, 1965. 


5. Krichevskiy, R.Ye. In the collection: Problemy Peredachi 
Informatsii. Vol. 4, No. 3, 1968, p. 48. 


33 


MULTIPURPOSE INFORMATION COLLECTION AND 
PROCESSING SYSTEMS 


A. V. Kantor, S. M. Perevertkin and 
T. S. Shcherbakova 


ABSTRACT. We examine the construction and 
present a classification of multipurrose information 
collection and processing systems which are based on 
the use of information sorting with the aid of as- 
socistive memories. Multiaddress compression, 
multichannel compression, and multipriority com- 
pression processes are described. 


The problems of space telemetry impose urgent requirements for /20 
the development of the so-called multipurpose collection and pro- 
cessing information systems. Here we mean by multipurpose data 
collection and processing systems, those which compress multichannel 
telemetry information with subsequent formation on board the space- 
craft of several compressed data streams with prespecified charac- 
teristics. The requirement for development of systems of tnis sort 
is due, first of all, to the constantly increasing volume of meas- 
ured information, increasing from experiment to experiment, t‘hich 
Must be transmitted and recorded (with subsequent retransmission) 
over a few channels with limited handling capacity, the necessity /2l 
for integrating the telemetry system into the spacecraft control 
loop, and many other factors. In the present article, we consider 
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three possible forms of multipurpose information collection and pro- 
cessing system organization: multichannel, multtaddress, and multi- 
priority. 


Dependilup on the multipurpose system organization form adopted 
at their output, there are formed from the stream of non-equally- 
spaced significant multichannel telemetry information reading: one 
or more streams of equally-spaced readings with one of the [following 
characteristics: 


1) in the case of multichannel organi ation of tne muitipurpose 
information collection and processing systems, “he compressed infor- 
mation output streams are formed from those significant telemetered 
parameter readings which provide the specified values of the approxi- 
mation error in each transmission channel; the number of output 
streams corresponds to the number of different values of the original 
analog parameter approximation error; 


2) in the case of multiaddress organization of the multipur- 
pose systems, the compressed information output streams are formed 
from the significant telemetered parameter readings with definite 
destination address; the number of streams is equal to the number 
of addressees; 


3) in the case of multipriority organization, the significant 
reading sequence in the output stream is determined by the value of 
the priority labels assigned to each message source group; the out- 
put may be a single stream, or the formation of several output 
streams is possible. 


Various combinations of the multipurpose systems listed above 
can be used to create quite flexible information collection and 
processing systems. We note that for creation of multipurpose tele- 
metry information collection and processing systems based on the 
traditional adaptive discretization schemes, the number of parallel 
output stream formation channels would obviously be equal tv the 
given number of these streams. Such multipurpose system construction 


can lead to an inacceptably cumbersome scheme. In the multipurpose 
information gathering systems which use associative memories (AM), 
which will be exainined hereafter, the formation of a whole series of 
compressed information streans 18 accomplished in a single channel, 
1.@., using the same equipment. Individual questions related with 
the possibilities of using AM in telemetry information collection 
and processing systems have been discussed in {1, 2). 


The expanded block diagram of a multipurpose information col- 
lection and processing system using an AM is shown in Figure 1. The 
multipurpose telemetry information compression process is accomplished 







Collect.,| 
compress. 
calc. 
(readout ) 
of AL 










from MS 


R compressed 
data streams 


Figure 1. Block diagram of multipurpose information 
collection and processing system: 


MS — message source; DCCAL — device for compress~ 

ing and calculating associative labels; SRwAL — sig- 

nificant readings with associative labels; DFS — de- 

vice for forming streams; OSR — ordered significant 
readings 


sequentially by the device for coliecting and compressing informa- 
tion and calculating (selecting) associative labels (DCCAL), and the 
device for forming compressed data output streams (DFS). We note 
that, in the case of multichannel organization of the multipurpose 
information collection systems, the tclemetered parameter approxi- 
mation error is the significant reading associative label, in the 
multiaddress case the destination address is the associative label, 


and in the multipriority case the message source priority rank is 
the associative symbol. 
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The Output stream formation prececs includes both general 
aczsociative sorting of the compressed telemetry information and the 
brocezs cf sorting the ordered significant readings with respect to 
the output enannels with subsequent smoothing of the nonuniformities 
oi these streams in a smoother. Jn the associative corting process 
(Fisure 2), there is formed a significant reading sequerce, ordered 
in decreasing associative label order, termed hereafter the ordered 
significant reading table. Sorting of the ordered significant read- 
ing table with respect to the output channels leads to the formation 





Figure 2. Associative sorting classification 


of R (one for each output channel) ordered significant reading se- 
quences, each of which includes readings with corresponding associa- 
tive label values or ranges of values. In the following, these se- 
quences are termed transmitted significant reading tables. The com- 
plete transmitted significant reading table includes, without excep- 
tion, all readings from the general ordered significant reading 
table. The formation of the ordered significant reading table can 
be accomplished either in the process of entering the information 
into the AM, if an AM with input sorting is used, or in the process 
of reading the information from the AM, if an AM with ordered read- 
out is used. 


In the case of information sorting in the process of entering 


the readings into the AM with input sorting, we can speak of simple, 
combined, and combined-with-priority sorting of the readings. In 


37 


the case of simple information sorting, the readings are entered 
into tne AM as fcilows: 


1) in decreasing associative label order; 


2) without one-to-one correspondence between the cell group 
address and the associative label values entered in any cell of 
this group; 


3) with loss of readings corresponding to the minimal associa- 
tive label values in the case of AM overflow. 


In the case of combined information sorting, the significant 
readings are entered into the AM in the following sequence: 


1) in dcsoreasing order of the associative label values; 


2) with one-to-one correspondence between the cell group ad- 
dresses and the associative label values of the significant readings 
entered in any cell of this group; 


3) with loss of readings with given associative label value if 
all the cells of the group intended for entry of readings with asso- 


ciative label values lying in the limits of the specified range are 
filled. 


In the case of combined information sorting at the inlet of 
the AM with priorities, entry is accomplished: 


1) in decreasing order of the associative label values; 


2) with one-to-one correspondence between the cell group ad- 
dresses and the associative label values of the readings entered 
in cells of this group, provided overflow of the cell group intended 


for entry of readings with large associative label value has not 
taken place; 
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3) with loss of readings with minimal associative label 
values. 


Thus, in the simple information sorting case, we obtain float- 
ing distribution of the readings with respect to the associative 
label magnitudes in the associative memory. In the combined sorting 
case, we obtain a fixed distribution, and in the case of combined 
information sorting at the input to the AM with priorities, we obtain 
a mixed reading distribution which coincides with the fixed distri- 
bution in the absence of cell group overflow. ; 


When using an AM with ordered readout, associative information 
sorting is accomplished in the process of readout of the readings 
from the memory in decreasing associative label value order. In this 
case, we can speak of quasi-fixed distribution of the significant 
readings stored in the associative memory in the process of ordered fou 
significant reading table formating using the cyclic sampling method. 


These associative sorting modes are shown in Figure 2. Figures 
3, 4, 5 show schematically the ordered significant reading tables 
for the AM input or output associative sorting modes examined above. 


The ordered significant reading table is represented in the £25 
form of a column of rows. Each row corresponds to a significent reading 
with some particular associative label value, entered either in the 
corresponding AM cell (AM with input sorting) or in the AM output 
register (AM with ordered readout). 


Significant reading sorting with respect to the output channels 
(Figure 6) for multipurpose information compression in the genera] 
case can be accomplished using the corresponding program in accord- 
ance with the selected output stream formation criterion. This 
sorting is based on combining, in each output stream, significant 
readings with different associative label values. In the case of 
formation of only a single stream of significant readings which 
follow in a definite sequence, we shall speak of the simplest 
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transmitted significant 
reading table formation. 
In the general case, as 
mentioned previously, the 
principle of combining in 
each table of readings 
with a definite range of 
associative label values 
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Classification of output 


Let us examine some variants of transmitted significant 
reading table formation: 


1) each table includes significant readings with strictly fixed 
associative label value; the number of tables is equal to the number 
of possible associative label values (multiaddress organization of 
multipurpose information collection system); 


2) a single table which includes readings with the complete 
range of associative label (multipriority organization of multi- 
purpose systems); 


3) R compressed message output streams or R transmitted signi- 
ficant reading tables. The first table includes readings with the 
largest associative label values in each formation cycle, the second 
includes those appearing in the first table and also readings with 
associative label values which are closest to the largest in the 
same formation cycle, and so on. The last table includes the maxi- 
mal number of significant readings transmitted over the radio link 
during a time equal to the table formation cycle. The last table 
is intended for transmission over a channel with maximal information 
content. We call this type of output stream formation sorting with 
sequential merging; 


4) R output streams, where formation of the reading tables for 
these streams takes place following a definite algcrithm: in this 
case, we speak of algorithmic sorting. 


Smoothers with complete and partial register sets (Figures 7, 
8) can be used to smooth the nonuniformities of the compressed data 
output streams formed for transmission. In the first case, for 
readings with each associative label value, the matcher has as many 
output registers as there are output streams to which it is con- 
nected. In the second case, for readings with corresponding asso- 
clative label value, there is a single output register; however, 
the control schemes are considerably more complex. 
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Figure 8. Matcher with iucom- 


Figure 7. Matcher with com- 
Plete register set 


plete register set 


In conclusion, we note that preliminary calculations made by 
the present authors make it possible to hope that creation of multi- 
purpose systems using AM will provide definite advantages from the 
viewpoint of size and weight characteristics in comparison with 
similar systems constructed by traditional methods. 
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ASSOCIATIVE COMPRESSED INFORMATION OUTPUT STREAM 
FORMATION BY THE STATISTICAL TRIAL METHOD 


A. V. Kanton, T. A. Tolmadzheva 


ABSTRACT. We examine a technique for studying 
the process of forming several compressed information 
streams, and present some results of such a study. 

We obtain the dependences of the output compressed 
information stream characteristics on the input stream 
characteristics, and also the dependences of the asso- 
ciative memory operation characteristics on the input 
and out »ut stream characteristics. 


References [1 - 2] considered the general questions of construct- 
ing multipurpose systems for the collection and processing of multi- 
channel telemetry information and made an analytic study of output 
stream formation for a channel with maximal information content. Re- 
sults were presented in [3] of modeling of multipurpose information 
collection systems with actual telemetry information which showed the 
advisability of more profound study of compressed information output 
stream formation. 


Use of the statistical trial method [4] in studying multipurpose 
information collection systems makes it possi.le to model on a digital 
computer the process of output stream formation for a two-dimensional 
Significant reading input stream with arbitrary distribution laws of 
both the time intervals between readings and the magnitudes of the 
approximation errors corresponding to each of these significant read- 
ings. Of significant importance for analysis of the output stream 
formation process is the establishment of the relationships vetween 
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the characteristics of the incoming significant reading stream and 

the characteristics of the output streams; at the same time, there is 

considerable interest in obtaining relations characterizing the opera- 
tion of the associative memory, which is the basic element in realiz- 

ing multipurpose information compression. 


The present study is devoted to obtaining these relaticnships by 
the statistical trial method; we model on a digital computer one cf 
the multipurpose compression modes (see [1]) — multichannel compres- 
sion with number of output channels equal to the number of different 
approximation error values with floating significant reading distribu- 
tion in the associative memory, and with constant duration of the 
multipurpose compression system operating cycle. The basic character- 
istics of the incoming sample stream are two random quantities: 1) z 
is the magnitude of the time interval between neighboring significant 
readings (for the simplest stream the equivalent quantity is the quan- 
tity B, characterizing the number of significant readings per multi- 
purpose compression system operating cycle); 2) e€ is the value of the 
approximation error, characterizing the significant sample. The two- 
dimensional incoming significant sample stream model is used in model- 
ing the associative output stream formation processes [5]. The number 
B of significant samples per operating cycle is distributed in ac- 
cordance with the Poisson law with parameter A: 


P(t) =O, (1) 
where b is the value of the random quantity B; t is the operating 
cycle duration (Tt = const); At is the average number of significant 
samples during the time t; P, (b) is the probability that B = b. 
Therefore, the probability density P (z) of the time interval z 
between neighboring significant samples is 


P (2) = Ac-t, ( 2 ) 


In the modeling, we assumed that the approximation error €, charac- 
terizing each significant sample, is a discrete random quantity Ex 
given by the distribution series: 


“ 
NO 
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where Ya is the maximal number of different error values. The 


Studies were made for a distribution series of the form: 
Pr ~ Pye’, 1 S LN — Ainax: (3 ) 


With account for the normalization 


*inax 
> ’xat 
Kel 
Expression (3) is written as follows: 
4 


Kmax 


wee (4) 
2, 


Pr= eck, 


In the modeling process, the incoming stream is represented in the 
form of the sum of Kg one-dimensional streams. The characteristics 


of each of these streams are the quantities Nex» defining the number 
of significant readings with approximation error ¢€ = Ex written in 


the associative memory during the pth operating cycle, i.e., during 
the time t. Then, 


Finan 


> Nex 
Ant = M(n,} = 2, 


Fmax (5) 
where rx is the average density of the significant sample stream 


with error Exe It is obvious that 


Kimax 
ao. (6) 
kh=1 


The basic characteristics of each of the R output streams are: 
1) (f) pig — the clock frequency, equal to the significant 


sample repetition rate in the absence of dropouts; (R - $8) — the 
Stream number, 0 < 6 <¢ R - 13 
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2) (e)p_¢ — approximation error, characterized by the follow- 


ing random quantities: 


a) (Erb Rg — approximation pseudo-error, values of the error 


e of samples entering the output stream with number R = 8; 


b) (4) Rg — maximal value of the error e of samples entering 


the output stream with number R = 8 +1; 


c) (epg —- value of the error e of samples which are not 


transferred to the output stream with number R —- 8; 


3) (Koon R-8 — contraction coefficient of the clock frequency 


in ti. stream with number R - 8 relative to the stream with number R. 
The quantity (f)p_g is defined by the relation 


Ra 
= 





PAL 1 
Ura = i 7? 0<s<k—1, (7) 
where Rg is the number of associative memory cells from which 28 


information 1s added to the information of the stream R = 8 - 1 
to obtain the stream R = 8; 
R 


D a] 


=1 


is the maximal number of significant samples entering the pth output 
stream during each cycle (see below); 
| ree | 
(Koon) Rag = So = neat ® (8) 


7 
TL 


It 48 obvious that for the R°" stream, (K =. 


con)R 

Additional characteristics of the neh stream will be the quan- 
tities 5S, defining the number of ordered significant samples actually 
transferred to the RtM channel in the rt) operating cycle, i.e., 
during the time t. 
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The basic charactertstics and parameters of AM operation are: 


1) the characteristics of the random quantity Mey: defining the 


number of significant readings stored in the associative memory im- 


mediately prior to initiation of sampling in the pth 


it is obvious that 


operating cycle; 


Men = Nex eb Mess 


where Mex is the number of significant readings which are not trans- 


ferred in the (r - 1th operating cycle; 


2) the characteristics of the random quantity Mo» characteriz- 


ing the length of the reading sequence, where 


Kinax 
Me == Mes 
e 2 tx 
3) tT — operating cycle duration; 


h 


4) fp ~— clock frequency in the Re output stream; 


5) 1 — maximal number of significant samples entering the 
Rth output stream during each cycle; it is obvious that 


I, 
fe H form, fr=—>; (9) 
~([% for mI, 


6) Poy — Probability of AM overflow: 


Poy =i (me >Epy hs (10) 


where EAM is the AM capacity; 


7) (Poy Reg — probability of overflow of the associative 


memory zone, the information in the cells of which 1s used to form 


47 


the stream with number K - 8: 


Bmax n--4 
Pog ea Pf Ym) >( 2 oy)} (11) 

em Se | wal 

8) Pam 77 Probability of AM emptying: 
Poy= P {S< Ml} =1—P (8 > 11). 2) 


The sought relations are shown in the table. 


TABLE 


Characteristics of random 
quantities Nexis Mey End? €433 













Argument Constants 





Quantities: Po¥3 Pos (Poy) Rog 


The following notations are used in the table: Ch is the nomi- 


nal value of the coefficient C in (4); the nominal value of C is 
defined by the relation 


Kyn 
I one (13) 
K pa: aes =e K=1 


2 ——_—_——. 
con 1 PK max go Ninax yp) 


This expression is obtained from (4) and (8) for & = R = 1, and 
p = 1; 9 is the load, defined by the expression 


pate, (14) 
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We shall describe the tloeck diagram of the algorithm for 
studying the output satrean formation process. Figure 1 shows the 
expanded bloek diagram of the program algorithm for studying the 
output stream formation process. The program includes the following 
basic subprograms: the input stream formation program (ISP); the 
Significant reading zone information processing program (SRP); the 
ordered significant reading zone information processing program 
(OSRP); and the result output program (ROP). The following nota- 
tions are used in Figure 1 [in addition to those indicated in the 
text and the Relations (1) - (12) presented above]: 1 is the number 
of the random number in the ISP subprogram (Figure 2); th is the 


th 


r cycle termination time; r is the number of cycles. 


max 


In eact: operating cycle r the significant readings are gener- 
ated using the ISP subprogram (Figure 2), the time intervals Zs 


~~ 
lwo 


are determined by (2), and the values of the approximation error ex 
are determined by (4). The quantities Z, and €, are formed from the 
random numbers Ky with uniform probability distribution in the in- 


terval 0 = 1, generated by a digital computer using a standard pro- 
gram. Upon termination of the interval tre - the stream formation 


ends. However, the possibility of appearance of the last reading 
after the moment th ls not excluded. This reading is stored in the 


transient reading (TR) cell, and is used in the next cycle. The 


values of the reading occurrence time t = ty - trel (t, is the time of 


reading occurrence after the interval 2,4) and the error inagnitudes /34 
€, are stored in the machine memory in the order of their entry into 


the zone, called the significant reading entry zone. Then the in- 
formation is processed using the basic SRP subprogram (Figure 3) and 
the auxiliary subprograms in the SRP for calculating the mathematical 
expectation M {ney} and M {mex} (Figure 4), histogram H {ney} and 


H {mex}; dispersion D {nex} and D {me}; and also the probability 


(11) of overflow of the zone (R = 8), (P » mathematical 


ov R=8 


op Input 
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expectation M {m_} and dispersion D {m,} (Figure 5). The process-= 


ing preceeds sequentially by rows (cells) n. Upon termination of 
processing, formaticn cf the ordered significant readings (OSR) 

zone is accomplished, i.e., entry into the individual digital com- 
puter memory zone of the information from the significant reading 
zone (without destruction): the information is arranged in the 
ordered sirnificant reading zone in order of decreasj.igthe e values, 
and for equal ¢€ in order of increacing t. This proccss simulates 
associative information sorting in the AM. Then information 
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Figure 6. Block diagram of algorithm of subprogram for 
processing information in OSR zone 


th 


vacated are filled. This terminates modeling of the r operating 


cycle, and transfer to the (r + re ek cycle takes place. Upon ter- 
mination of information processing using the subprograms ISP, SRP, 


OSRP (r = Pak)? the results are processed using the basic ROP 


subprogram (Figure 7) and several auxiliary subprograms, specifi- 


cally, the subprogram for calculating (Poy Rg and Fxg (Figure 8). 


During processing using the ISP, SRP, OSRP subprograms, data are 
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Figure 11. Forma~ Figure 12. Forma- Figure 13. Forma- 
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48; c/C,, -w1,2 Keon 1” 323 p = 16; 3 — c/c = 1, 2; 
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accumulated for calculating the mathematical expectation, dispersion, 
and histogram of the random quantities Megs Mexs Enns Ens ©4373 Sy 


and data for (Poy Rog and Poy: The ROP subprogram is used to cal- 


culate the aforementioned characteristics based on the operating 
cycle tak values. 

The results of the digital computer study of assoclative output 
stream formation are shown in Figures 9 = 13. The data obtained re- 
flect quite completely the relationships between the incoming signi- 
ficant reading stream characteristics and the output stream char- 
acteristics, and also characterize the associative memory operation. 


In conclusion, the authors wish to thank I. I. Golovchenko and 


R. A. Sagitov for their assistance in developing the program and 
making the computer calculations. 
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ANALYTIC STUDY OF OUTPUT STREAM FORMATION PROCESS IN 
MULTIPURPOSE INFORMATION COMPRESSION SYSTEMS 


A. V. Kantor, S. M. Perevertkin and 
T. S. Sheherbakova 


ABSTRACT. An analytic study is made of the 
output stream formation process for a system with 
maximal information content. We take as the math- 
ematical model the single-channel mass servicing 
system with very simple input stream, with group 
demand servicing, and Erlangian servicing inter- 
val distribution. Analytic expressions are ob- 
tained for several system state probability dis- 
tributions and queue length mathematical expecta- 
tion and dispersion which can be used for certain 
simplified calculations of the formation process 
of all the output streams, but with the use of 
priority servicing and with a simplified servicing 
interval model. 


As was shown in [1], it is advisable to study the multipurpose 135 
information compression system by queueing theory methods. In this 
case, such a syStem can be considered a queueing system with waiting. 
A stream of significant readings with corresponding essociative 
labels enters the input of the associative memory (AM) [2] used as 
a sorter and buffer memory (BM) in multipurpose information com= 
pression systems. The significant reading stream is taken to be 
very simple (experimental studies confirm the validity of this as-= 
sumption for stationary segments of telemetry parameter behavior. 


57 


The output of readings from the AM, 1.e., servicing, is acconm- 
plished by groups of I] readings during the multipurpose information 
compression system operating cycle. 


In this study, the significant reading group servicing time is 
assumed to be distributed in accordance with the Erlang law with 
density 

phy Ant 


f(t) = F(R) eu, 


(1) 


1 
where =a is the average number of significant readings taken 


from the AM into the matct.er register per unit time; K is the para- 

meter of tae Erlang distribution. Specifically, for K =1/(t) = ue, 

which corresponds to the exponential distribution, for K = © we have 
the constant servicing duration case examined in [1]. 


In analyzing the states of the subject mass servicing system, 
we use the embedded Markov chain method [3] in which the system 
states are examined at strictly defined moments of time, namely, 
immediately preceding the moments of servicing termination. We note 
that rere, as usual, by system states we mean the number of messages 
(demands) in the associative memory at a definite moment of time. 


The Chapman-Kolmogorov equations are valid for embedded Markov 
chain points: 


F \1 A 
Pall) = BP — 1) Pn (2) 
where P,(/) is the probability that at the end of the jth cycle, 
the system will be in tne state n; ?P,(/—1) is the probability that 
at the end of the (j - 1)th cycle, the system will be in the state 17; 
Pn() is the probability of process transition from state 2 into 
state n, where 


Pin ()) = P(N () & WIN (f — 41) = U 


Let qn be the probability that in the gen eycle, the number of 
newly arriving readings will be b. It 1s obvious that for 0clell, 
ben, fori<cticu-ra,b--n—li-, and for t>ii+rn, b does not 
define the system state at the end of the jth cycle, i.e., 


Qn for 0<li<ii—1, 
Pin Qn-i + HT , for Usl<e i +n, (3) 
0 for M+ncl< w. 


The corresponding transition probability matrix has the form (the 
subscript j is omitted for simplification): 


Fo 91 Ja Ys | 
rr rr er er ee Wu—-i1 
Go M1 92 Is | 
H i = 
Pint Go 1 Is Is 
O go 1 92 


00 mn. et! 


se 2 © eo we ee 





With account for (3), Equation (2) is written as follows: 


n-1 IE-|-n 


Pu) = DPM — Daa + BD Pd — VN gneten(f) = 
t=9 tasdf 
n~-) +n (4) 
= 9n (i) Py Pj 1+ Py Pf —A) dnetsn (I). 
oo Dn >0. 


The probabilities Pa (Jj) satisfy the normalization condition 


Pa) ad 
2 (5) 
For the steady-state regime (joo), the system of equations (4) will 
have the form 7 

Po -= J > Py +Pudeo 


IQ 
m1 


Py h >y Py + Pug -} Pa + 190 
t=30 
ee ara ee aa re a (6) 


-1 Wen ~ 
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In solving the system of equations (6), we use the generating 


tion method. We introduce the generating functions: 


P(Z) = YP, 
nO 


Q(4)= Dae". 


a 0 


Substituting the values of P, from (6) into (7), we obtain 


~ . dist Wren 
PZ) = D2 in Br+ D PQunn| = 
n=O =0 txif 


a n=—l] w 
= 3 a2" Py Py + > 2 (Pndn + PrnQnea- ++ + PrtenQo) = 
n=0 =0 n= * Pe 
ti-1 fw 


= 0 (2) 2 Py ie D2"(Padn + Piterdner b+» + Prtenta) 
=) a 


fure- 


(7) 
(8) 


(9) 


The convolution of the two sequences {¢n} and (/..} appears in the 
brackets in the right side of (9). Unwrapping this convolution, 


we obtain: 
di--1 w o 
P(Z):= Q(Z) Y P+ (240) DY Z"Paen 
te-0 fim n=0 . 


where 
Son : 
> 2", == Q (2). 
n=0 


The sum )) Z°Piyn in (10) can be transformed as follows: 


nav 


x 1 
Y2"Pinn = 7" |P(2)-- Pz. 
n=Q tead , 


The Expression (9) can be written in the form 


1 


Met 
P(2) = (2) Yr FOS" [PH — Z Pz. 
i=0 . jag 
Solving (13) for P (2), we obtain: 


M1 
N pygZit — 2!) 
P(Z) .2 el) AL e 
om —! 
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(10) 


(11) 


(12) 


(13) 


(14) 


In the Erlangian servicing time distribution case, the Rxpreaston 
(8) for Q (Z) takes the form: 


9 (2) = [4 ABE (ty SEE AYE, (15) 


where p= 4A li is the load; K is the parameter of the Erlangian 
distribution; A is the incoming stream intensity; I is the number of 
sipnificant readings in th group. Then (14) takes the form: 





Mat 
Np, (aiiz') 
? 9 - . 
P (2) 2n(s 4 ea) 2s (16) 
Omitting the intermediate transformations, (16) may be written in 
the form: 
Cc 
P (2) = 75-1 
it @,-2) oy? 


jun 


where C is a constant; z5 are the roots of the denominator of (16) 
outside the unit circle. The constant C is determined from the con- 
dition P (Z) = 1 for Z = 1, hence: 

U+hK—1 


Cm [] (%-1) (18) 


jan 
Finally, we write the expression for the generating function P (4%): 


U+Ke1 , 


a= Ml aaa 





(19) 


To find the system state probabilities, we expand the expression 
for P (Z) into simple fractions: 





¢) ‘+ Ke’ 
PQ) a ging + eg tee ee te pity, (20) 


where the exvansiton coefficients BY are found from the expression: 


a 1 ad 
a oe 
By = ( 21 ) 
A im al | 
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Varying v in the limitam+K-1pveal, and substituting the corre- 
sponding expressions for (21) into (20), we have: 


“Before rihvteele- 
a real +h at(gee) te (22) 


We obtain the expression for Py from (21) as coefficients of the 


corresponding powers of Z. Thus, 


Bu 


Br Bris no 
y= = zy i rae +:- ae 


Tika (23) 





If Zp is the smallest-in-magnitude root of the denominator of (16) 


outside the unit circle, then the approximate expression for Pe 
has the form: 


Ba 


i= Zt” (24) 


The mathematical expectation of the number of demands in the system 
and the dispersion of the number of demands are found from (22): 


aP(Z 
M{N) = ap Ea (25) 
and 
d*p(Z) cil dP (2) 
DiN} = =o leak = od (ae si) (26) 


Direct calculation of the first and second moments using (25) and 
(26) is difficult. We use the semi-invariant method of [4]. 
Setting Z2=l, we have 


M (0) = 2 Pt", (27) 


where M (6) 1s the moment generating function. 
Taking the logarithm of (27), we obtain the generating function 
of the semi-invariants » (8). The mathematical expectation of the 


number of significant reauings in the system and their disnersion 
are obtained by differentiating the semi-invariant generating 
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~ 
"a> 


function, setting @ = QO: 


fies Wik a 
+f “ es ’ yo. Aye 
MIN} = _ es (4; «Ay 


Uhm 


DiNy 2M) OY yay ty. (29) 
oa Zul 


(28) 


Let us examine some limiting Erlang distributions: 

1) A: L1gH<o (servicing interval distributed exponentially); 
2) KK: ole to (servicing interval constant); 

3) # -1, 1S AS (single servicing); 


4) He xn, ISAKSN (group servicing with servicing after a 
single interval of an infinitely large number of SR). 


For the limiting Erlang distribution cases presented above, we 
obtain the following expressions for the generating function, dis-= 
tribution series, mathematical expectation, queue length, and queue 
length dispersion: 


2, —-1 
4. P(2) = ; (30) 
Se Fit) (31) 

Py, = Znat e 


Here, Zn is the single zero cutside the unit circle of the denomi- 


nator of (16): 


M{N} = ping, (32) 
; Zz 
DiN}= Gay (33) 
N-1 hed 
Hd - oe |] 
2. P(Z) = noe : (34) 
. (2) = ——~"pieni=ay + 


where Zs are the roots of the denominator insiue the unit circle and 


on its boundary. 
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Cane 2 Pr aeee 2 
ayn } Tr aH) m * 


inkelay 


1- dt 


wey my . 
M (Nps To » " wig 


Ht. -4 
"i — p)b ms 2814 — pt 5 Y, 
DN) = (1 4-20) + Gall (1 — p)*-- 8 (1 — p) —) 4“ 
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121 — p) 


(i Zi) * 


3. aligAg>. 


kK 
Z;-! 
PZ = “4 en 9 
(2) Als % 


“ 8 j= 
Pa = Dy zee , Be a 4 i ) 
a il eect lent 


h 
MN} = DY (2-104 


iat 

DIN} == D) 2)(Zj — 1). 
Pi ¢ 

b) W214, 4 =1. 
Zi 1 

P(Z)= =FT—Z' 


4 or | 
Preeti Set 
a an s] 





4 
MIN) 8 Ft 
4 
DN = aR 
ec) H-=4, Ka, 
z t= p)iZ—t) 
P(2) = Lea 4 
i-p for n=0, 
CaNes) for n=1, 
P= n 
e AY n-S Si perso 1 
(i —0) Bi Ayes eb (SO + or! for n>?, 
a— 
tie Row “p)* 
D{Np-s Mbt ol a 


12 qi --p) 


~“, Lc Kx “ 


(35) 


(36) 


(37) 


(38) 


(39) 


(40) 


(42) 


(42) 
(43) 
(44) 
(45) 


(46) 


(47) 
(48) 


(49) 


740 


6°70 = d — qQ fg°0 = dO — B 


:dnodd ut 
sdutTpead jo azaquinu snsdaA UOoTS 
-i9dstp uqguat anang ‘2 aunstg 


az g2 9+ 2t 8 4 a 
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Upon substituting into (19), in place of Z, the quantity 
(1 - 2/M), we obtain the expression for the generating function 
P (2) [5]: 
k 
{1 (Z; 1) 


Equating the denominator of (16) to zero for S = I {(Z = 1}, as 
Il + », we obtain: 
es ( ce ee 


With account for the roots of (51), we have for the generating 
function Pp (2.3 





K 

- s 
lim Py (Z) = 
pe aN) ar 


Zz » 


where Sy are the roots of (51) for Re (S) > 0. 


in 7 5 


— ane 
ca a, Be +r wv) 


wnere: 


lim MInt tN} y Si 


di oom ey 
Ly S;3. 

nA Ds 

i i= 


b) fi=w, — 


> a3 
tim Pu(Z) = sare 
lim —+ = s 
Hen il ‘ene 
PU | 
wn é 
tim S4S2  g73 
fen 
N 3 
lim 2132. g5%. 


Ww 


(50) 


(51) 


(52) 


(53) 


(54) 


(55) 


(56) 


(57) 


(58) 
(59) 


~ 


Nd 





c) il == ~, K = », 


lim P,(Z) = e7*, (60) 
ii-n 
. P, 1 npn oll ; 
Hea are es oo 
o MAN 
tim “1 =o, (62) 
- DIN 
lim Pt) = 0. (63) 


Calculations were made of the functions M{N}/? and D{N)}/JP as a 
function of 7 for two values of the Joading pe = 0.8 and 0.9 (Figures 
1, 2). 
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DISPERSION SPACE RADIO LINKS 


L. G. Sapogin and V. G. Sapogin 


ABSTRACT. We examine the dispersion method of 
information reception and transmission with use of 
the dispersive medium in which the radio waves 
propagate as an optimal filter. 


Highly directional ground-based antennas with large etfective 
areas and cooled high-sensitivity receivers are used for communica- 
tion with distant spacecraft. With the narrow antenna pattern and 
weak signal, spacecraft search becomes a complex problem because of 
the long time for establishing communication. The uplink communica- 


tion problem is simpler than the downlink problem for the following 
reasons: 


a) low onboard transmitter power; at the present time a trans- 
mitter power of 109 W can te considered the limit, since increase of 
the distance from the vehicle to the sun leads to marked reduction of 
the onboard energy supply; 


b) limited directivity of the onboard antennas because of the 
severe requirements on spacecraft attitude stabilization in space; 


ec) 1:mited ground antenna effective area; 
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ad) the tremendous distances and inexorable law of Squares 
lead to a situation in which the communication problem at the pre-= 
sent time is reso)ved using all possible approaches (combined method). 
Further technology development requires strenous efforts and large 
expenditures in all these directions at the present time. 


It is well known that the outer-space and near-Earth plasmas 
in which vadio waves propagate have dispersion, which has always been 
considered a serious obstacle for transmission of wideband and FM 
signals, since they are distorted strongly because of dispersion [1]. 
Such dispersion-type signal distortions increase with distance, and 
impose limitations on the transmitted information spectrum width. 
For lunar distances, the dispersion limitations require signal spec- 
trum width Af < 40 MHz, and at distances of billions of kilometers 
the required bandwidth is on the order of tenths of a Hertz. The 
existing long-distance radio communication techniques have ignored 
the dispersion properties of the medium and have considered only its 
attenuatiun. The present study shows the possibility of using the 
dispersion property of the medium in which radio waves propagate as 
an optimal filter. For a given medium, the signal spectrum can be 
selected so that the signal harmonic components, traveling with 
different phase velocities because of dispersion, combine at the 
reception point in the required phase and create a local spatial re- 
gion in which the peak signal power is very high. Here, the disper- 
sion limitations will have a different nature, and the proposed 
dispersion method of information reception and transmission may lead 
to some progress in the questions of ultra-long-range space 
communication. 


GENERAL QUESTIONS OF WAVE GROUP PROPAGATION 
IN DISPERSIVE MEDIA 


First of all, we shall clarify how wave groups (signals) 
propagate in dispersive media. This question is not new, and has 
been examined in connection with study of electromagnetic wave 
nropagation in feeders [2], plasma [3], amplification of light pulses 
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in radar [5, 6]. If we start from the problem of compressing a wave 
Group into a narrow pulse [7] (in the limit a 6-pulse), we can formu- 
late two problems: 


1) finding the form of the matched wave group (spectrum, modu- 
lation law) which is compressed by a system with given dispersion 
characteristics into a narrow 6-pulse; 


2) finding the required dispersion characteristics of the 
system from the given form of the signal which is compressed into 
a narrow 6-pulse. 


The first problem is of primary interest for the creation of 
dispersion-type space radio links, and the solution of the second 
problem is analogous. We shall in the following assume that the 
dispersive system is linear, i.e., its characteristics are independ- 
ent of wave group amplitude. In nonlinear systems, the question of 
wave group propagation becomes exceptionally complex, and has been 
examined in part in [8]. 


As is known [9], if a signal v (t) acts on the system input, 
the output signal will be 


u(t) _— \ K (a, 1)S, (a) et do, (1) 
where 


Sp(w) = J v(tye de 


is the input signal spectrum; K (w, 2) is the transfer function of 
the linear system. A linear di:persive medium or system can be 
represented as a symmetric quadripole with two independent parameters: 


wl 


i—- 
K(o, bl). sole wey i 


where @ (w) 1s the coefficient of attenuation or amplification of 
the medium; v, (w) is the wave phase velocity in the medium (the 
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dispersion characteristic of the system); and Z is the length of 
the considered segment of the medium. 


We shall examine two different dispersive systems K (6, 2) and 
Ky (6, 2), and two wave groups: broad X (t) and narrow Y (t) with 
spectra s, (6) and Sy (6), respectively. Assume that under the ac- 
tion of the broad group X (t) on the input of the system Ko (6, 2), 


there arises at the output the narrow group Y (t), and that under 
the action of the narrow group Y (t) on the input of the system 
Ky (6, 2), there arises at the output the broad group X (t) (reversi- 


bility). It is obvious that: 





S50) (3) 

Kg (, l) — Soy , 

S,(o) 
K, (0, l) = So) (4) 
Hence, it follows: 
1 K,(o, 0 (5) 
Ke(0, ly = K,@.0 a =] pl DP? 

% (0) = — Gp»), (6) 
bj), (w) = — vy,,(0). (7) 


We see from (6) and (7) that the characteristics of the "narrow- 
ing" and "broadening" systems are mirror symmetric relative to the 
w axis. The transformations performed on the wave group are not 
equivalent to reversibility of the propagation process (i1.e., re- 
placement of t by -t). While the system with normal dispersion and 
attenuation a broadens the narrow pulse, in the system with anomalous 
dispersion and amplification a, this broad pulse will be compressed 
into a narrow pulse. We note that all processes in systems with 
a # 0 are irreversible. If the dispersion law is not altered, the 
sign of the time for the pulse changes. For example, if a linear 
FM puise with time-increasing modulation frequency is subjected to 
compression, then, as a result of broadening of the narrow pulse, /45 
we obtain a broad pulse with time-decreasing modulation frequency. 


7] 


We note that (6) and (7) are valid for mutual transformation of a 
definite narrow pulse into a definite broad pulse and vice versa, 
but not for an arbitrary broad or narrow pulse. In this case, the 
same dispersive system may narrow one pulse and broaden the other. 


FINDING THE SPECTRUM OF A WAVE GROUP COMPRESSED 
INTO A §6=PULSE IN DIFFERENT DISPERSIVE MEDIA 


Assume that in one case a wave group 6 (t) arrives at the input 
of the system K (w, 2) and is broadened into Y (5), while in the 


other case the wave group X) (t) acts on the same system and is com- 


pressed by the system into the function 6 (t = tT). Then we can write: 
+8 ; 
j Yi) eC at 
KQo,) == ———-, (9) 
i} 5 ctyerio dt 


+r 
J de —neMae 
Oe (10) 
J AEG) oll ag 


—»~ 


After integrating and equating the right sides, we obtain: 


+2 +2 
iy Y; (t) eviut dt \ x, (t) eniwt dt =: e-tut ( 11 ) 


or the relation between the spectra: 


td (ed 


Su(2) = 35)" (12) 


Since the spectrum of the 6-function: 


Solo) = J S(therietde = 4, (13) 


-s 

by comparing (13) and (1), we obtain the condition under which a 
wave group in 6-function form is obtained at the output of the 
dispersive system: 


K(,1)S85(9) = Sen) =, Kn, 2) >a: (14) 


This formula makes it possible tc find from the specific form of the 
dispersion characteristic the spectrum o: the signal which is com- 
pressed into a 6-pulse and, conversely, to find for the case of a 
specific signal, which is compressed into a 6-function, the required 
dispersion characteristic of the medium. The exact solution of this 
problem for a specific dispersion law and wave group always involves 
considerable mathematical difficulty in calculating the Integral (1), 
which, as a rule, cannot be calculated exactly, and 1s evaluated 
approximately by the steepest descent method or on computers. There- 
fore, 1t is better to proceed in the opposite direction — make the 
mathematical study for approximate dispersion equations, when the 
Integral (1) can be calculated exactly. 


For an approximate qualitative evaluation of the signals which 
can be compressed into a 6-pulse, we shall use the following simple 
analytic equations, obtained by graphical approximation of the actual 
dispersion characteristics in some frequency region. The linear 
normal dispersion: 


rae(1— Si, Oo wo, sala aie (oa (15) 
o 
The hyperbolic normal dispersion: 
The parabolic normal dispersion: 
ite hy ee Rea ep (Se) (17) 


Approximation of anomalous cispersions yields the same results, but 
with inverse law of frequency modulation within the wave group. 


We find X (t) for the linear normal dispersion: 


(18) 








The sought intergral can be tabulated. The result will be: 
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NOs pe “Bex [-i jee Jexp [iar a (19) 


ye i 7 


i.e., 1f a wave group of very high amplitude, with carrier frequency 
varying in accordance with a definite law, arrives at the initial 
moment of time at the input of a system with linear normal disper-= 
sion, then a éd-pulse is obtained at the output. Let us examine the 
same case for hyperbolic normal dispersion: 


exp [i ae cinta = |/ Ph oxp(— iF exp "toy “ : (20) 
‘ > ‘ 


Any wave group of constant amplitude and linearly increasing modula- 
tion frequency is compressed into a é=pulse. For parabolic normal 
dispersion the analysis is difficult, since the resulting integral 
cannot be evaluated. Evaluation by the steepest descent method 
yields the same result as for hyperbeiiec normal dispersion, but with 
nonlinear modulation law within the wave group. The estimates made 
above are very approximate because of the unrealizability of such 
dispersion equations for the condition a = 0. Moreover, the 6- 
function yields an infinitely broad spectrum, while real wave groups 
cannot have such a spectrum. Therefore, it is necessary to solve 
the problem for an output signal in the form of a real A-~pulse. 


FINDING THE DISPERSION CHARACTERISTICS WHICH 
COMPRESS SIGNALS INTO A REAL A-=PULSE 


The correctly introduced real A-function for wave group repre= 
sentation must satisfy the Laplace equation and be sufficiently simple 
so that the resulting integrals can be evaluated. We take an auxili- 
ary function y (t, p), for which p > 0 depends continuously ont, 
and in the limit takes the values: 


| ~- for t>0, 
0 fort=0, (21) 


lim y(t, p) = 
e a for <0. 


This requirement is satisfied by u function whose graph coincides 
with the function y = 1/2 sen x, 
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o my ne St sk £ 
v(t p) le = dw = 7 Orele (22) 
@ 


The derivative of y (t, p) with respect to t, which we denote by 
A (t, p), will be: 


tp) Lf aoe: 
Or) A(t, p)= +) e-® os wl dio = —— oro (23) 


It is obvious that 


’ 0 for ¢4=0, (24) 
porta P= |. for t= 9, 

lim ( A(t, p)dt =4, (25) 
P-0 Soy 


&(t) = him A(¢t, p). (26) 
p—-0 


In the case of conventional integration of this function, the limit 
passage must be performed after calculating the integral. In other 
words, the small parameter must have lower order of smallness than 
the inercment At. An integral with such a A-function denotes cal- 
culation of the limit of the sum as At + 0, p + 0 and At/p + 0. 
Such integrals are improper integrals, and the function utilized is 
a generalized function. 


Since the introduced function A (t, p) is the real part of the 
analytic function 


1 1 1 
3 a pr’ (27) 


It is a solution of the Laplace equation and, consequently, can 
describe real wave groups. We find the spectrum of the function 
A (t, p): 


wo 


eae eee = 
Sate. p) () ay aT pia e-Pwl dt = e-IPl, (28) 





Now it is not difficult to find the connection between the wave 
group and the transfer function of the system which compresses this 
group into A (t, p): 
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ot 
S(w) ° 


K(a, )S(m) sel, Ko, 8) = (29) 


As p + 0, the Equation (29) becomes (14). We note that the form of 
(29) changes very little with change of the technique for specifying 
A (t, Pp). 


COMPRESSION OIF’ LINEARLY FREQUENCY-MODULATED WAVE 
PACKET WITH GAUSSIAN ENVELOPF 


The complex espression for the linearly frequency-modulated 
wave group with Gaussian envelope has the form 


S(t) = exp(— £/T*)exy [é (gt +: 101, (30) 


where T means half the pulse duration at the 0.37 maximal amplitude 
level. For a linear w variation law, the instantaneous signal fre- 
quency takes negative values in some region. In this connection, 
we assume that: 


for] Ping (31) 


and the region of negative instantaneous frequency values corresponds 
to negligibly small wave group amplitudes. We find the spectrum of 
such a@ wave group: 


~~ 


S(o) os : S (te dt = { exp — [(+s tr, 4 i(o~a,)t dt. (32) 


We denote 1 - iyT = u, and obtain: 


(w — w,)8/'2 


S(0) = “exp =) \ ee" dz =: CVS exp [—Ganee|. (33) 


We find the modulus of the complex number u: 


m= (eye Pa Pre 


Separating the real and imaginary parts in (33), we have: 
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where 2m og VE 1s the effective spectrum width. Conzequently, 


the amplitude spectrum of such a pulse has Gaussian form and 4s cone 
centrated near the center frequency 4. Using, (29), we find the 


transfer function of the optimal dispersive medium: 

A(») = exp [- Car loxp (= [rms — a arcig wih ; (35) 
where |4(»)| - oxp| Saar) is the amplitude-frequeney characterintic 
of the medium and #(») soxp{-— dirt om _ + ace:a rtf} is the dispersion 


ek ae 


vcnaracteristic of the medium. 


If a linearly frequency-modulated wave group is passed through 
a medium with these parameters, at the exit from the medium we obtain 
an amplitude-modulatcd wave group of beli-shaped form, for which 
carrier frequency modulation is absent: 
U (t) - ae ‘ S (0) Bo) edo .2 me Tm elise 
ex ov —s ' 


Comparing (36) and (30), we can see that the output signal duration 
is reduced by m times and the amplitude is increased Ly Ym times in 
cuazparison with the input group, i.e., the pulse has been compressed 
by m times. Thts also follows from the energy conservation law. 


The signal energy at the entrance fc and exit from the disper- 
sive mediun is defined by the expressions: 


ey ; 
Ean = x \ {S(o)P do = Te \ o*(w) dw, 
ty 1 ¢ 
Eox = te \ Sido ye J st(o)da, 


The signal compression factor is equal to the ratio of the entering /49 
wave group duration to the exiting wave group duration: 


2 
ance Ton = Tenfen ent 
Tay By en en? 
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where 
\ 35 (19) den 
ae Tq) 
W > => » dof) es 
og it 
| % (4) ln 
-~< 


s the normed spectral density. Since op (w) << 1, then n > 1. 


Thus, the compresston factor is equal (in order of magnitude )to the 
product of the signal duration by its spectrum width: 


m2 T * Aw. 
nm an w 


For frequency modulatea signals with large deviation frequency 


Tg, the occupied band width Aw = ea/t as then: 


mo 21ften » 


(37) 
but this is simply the modulation parameter of the frequency-modu- 
lated pulse. It is clear that for a given dispersive medium increase 


of the duration en requires corresponding increase of the trequency 


deviation and this makes it possible, purely theoretically, to ob- 
tain any compression factors. In practice, the entire problem lies 
in radiating the frequency-modulated pulse with very large frequency 
deviation in accordance with a strictly defined law, since the 
stability of the medium at large distances is quite high. 


ON DiSTERSTON OF MEDIA USED AS OPTIMAL FILTERS 


According to the latest experimental data [3], interplanetary 
space is a plasma with particle concentration N = 100 el/em3. For 
such concentrations, we can neglect the magnetic field influence 
and consider the plasma nonmagnetic and collisionless. The disper- 
ston equation of such a plasma has the form: 


,s € : Ant 
J 3 ——=—-— ? (Ng = ——< 


2) 


i ae) ) (38) 





For the interplanetary plasma, the frequency Wy is low, and for 


suffictently high frequencies, we can expand (38) into a series 


TR 


in powers of W/W and consider only the first two terms of the 


expansion. Then, 
2 
Vv, elt Sy, (39) 


which is greater than the speed of light. Formulas (38) and (39) 
lead to good agreement with the experimental data for the inter- 
planetary and near-Earth plasma. 


In the general case, the dispersion of outer space is determined 
by the magnitude of the integral electron eonecentration (1EC), which 
for near-Earth plasma varies continuoucly, depending on the time of 
day, season, and solar activity cycle in the limits (1 - 8) - 1023 
el/en@. Change of the IEC may lead to change of the compression 
factor m (Figure 1). In addition, the magnitude of the IEC for near- 
Earth plasma depends or. the magnitude of the zenith angle 8 follow- /50 
ing the law (Figure 2): 





N, s Ny acc 8, 
0 
where Ny is the IEC in the normal direction. 
0 
" 
Z 
Mb 
Figure 1. Compression factor Figure 2. IEC as function of 
&s function of IEC and dis= zenith angle 


tance 


Studies made using artificial Earth satellites have shown cone 
clusively that, along with the regular daily variation of Nas there 


are observed irregular variations due to tonospheric nonhomogeneities. 
The ionospheric nonhomogeneity dimension spectrum occupies the 
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interval from several hundred meters to several hundred kilometers, 
and the magnitude of the electron concentration “luctuation increases 
with inerease of nenhomogeneity size and lies in the range AN/N ® 
0.1 = 50%. Studies of the Mars and Venus ionospheres during 
flights of the Mariner 4 and Mariner 5 spacecraft make it possibie 
to evaluate the magnitude of the dispersion of the fonospheres of 
these planets. Using the electron concentration distribution pro- 
files in the Mars lonosrhere [10] and Venus ionosphere [11], we can 
find the IEC along the radio wave propagation path and the rate of 
its variation as the vehicle approaches the planet. Calculations 
show that N, along a ray in the Mars ionosphere may reach values of 


N, ~ 103 el/em*, and in the Venus ionosphere — several orders 


16 2 


greater, N, \ 1015 - 10 el/em*. According to the latest data, in 


the Sun's supercorona region (outer region of the corona at a dis- 
tance greater than five sun radii) the electron concentration N 
decreases, following the law [13]: 


N, (R) = 1.1 + 1073 Ro?, 


where R is measured in meters. For super-long distance spacecraft 
studying the depths of the Galaxy, the IEC may reach 5 - 


107? ~ § - 1073 e1/em? [14]. 


On the other hand, the question of electromagnetic wave disper- 
sion existence in a vacuum has been posed many times, and is of 
definite theoretical and practical interest. Rozenberg's data on 
electromagnetic wave velocity measurement in a vacuum for a very 
wide frequency band after averaging yield the following table 
[15, 16): 


Wave frequency f, Hz 109 1024 1019 1024 
Propagation velocity 299,787.4 299,781.7 298,300 297,400 
¢, km/sec 


We see that the difference of the velocities at the edges of the band 
studied amounts to 1%, which corresponds to velocity change per unit /51 
frequency equal to: 
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T.. = 2 U,Sa8TA 10, 


More precise measurements made in the RF band confirm this value 

of T. Rozenberg himself [16] questions the existence of dispersion 
in a vacuum, in contrast with other authors. Specifically, Teller 
(19] proposed exploding an atomic bomb in space and measuring the 
difference of the time of arrival at the Earth of the light and 
gamma rays. 


It is not difficult to find the shift in seconds for frequencies 
of different electromagnetic waves emitted at the same time from 
some object: 





M=t—= fe = Hee oe, 


C1 ats 


where Cy and Cc, are the electromagnetic wave velocities at the fre- 


quencies f) and fg. Kotelnikov and other authors [3, 12, 13, 17, 18] 
have drawn the conclusions on dispersion in interplanetary space. 


DISPERSION LIMITATIONS IN LONG-RANGE 
SPACE RADIO LINKS 


When organizing long-range space communications, it is neces= 
sary to consider the dispersion properties of the near-Earth plasma 
and interplanetary space, since signals are distorted markedly when 
transmitted over long distances. Therefore, it is important to 
clarify theoretically the limitations which arise. Tne group velo- 
city of signal propagation in a dispersive medium is: 


a 
Ver en, 


where n is the refractivity of the medium. For a collisionless 
plasma, 


where fp = 8: 10'; N 1s the number of electrons; e and m are the 


electron charge and mass. Then, the time of signal propagation 
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between the points xX, and Xo for arbitrary electron concentration 


variation law will be 


where f is the signal center frequency. For re >> f° (condition uf 


radio transparency), this expression may be expanded into a power 
series, and we can write in the first approximation: 


iv stor 
t=—\de+ a7 \ ¥ dz. 
z Cry 





The first term 1s tne signal propagation time in free space, and 
the second is the lag effect from dispersion, which depends on the 
frequency. Any transmitted signal occupies some frequency band Af, 
which is connected with the signal duration t by the known relation: 


Mf=t. (40) 


Since the signal lag time is inversely proportional to frequency 
squared, we will observe a difference in its magnitude at the edges 
of the signal spectrum. Let the signal spectrum upper frequency be 


Patt eM » and the lower frequency f, = f— yA . Then, 
Cty 
by —t at a eee aed \ Ndz. 
ae : la rap 
If fS» +; which is always the case, then 
ator [y= faut i) Te g.sutaye 
At = scan lume t ome Ndz = = >!\ Naz. 
Zz u as Lyy 
If the time At is commensurate with or greiter than the transmitted 
signal duration, undetectable signal distortions and information 
loss occur. Expression (41) with account for (40) may be written 


in the form: 
—T— 
Af ; F 
a er (42) 


(41) 
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This inequality is the condition for undistorted signal transmission 
in long-range space radio links. For circumlunar distances [20], 
the quantity 


Xe 
\ Ndz = N, = 5-10" e1/om@. 
™ 


Distortions will not be observed for a signal frequency f = 1 GHz, if 
Af << 40 MHz. For vehicles located beyond the Sun's supercorona or 
used for studying the depths of the Universe, the quantity Af will 

be on the order of fractions of a Hz. 


Changeover to the optical band (use of lasers) would make it 
possible to increase the signal spectrum width for long distances, 
but it is difficult to count on the appea:sance of such communica= 
tions systems in the near future [21, 22]. 


INFLUENCE OF PLASMA NONHOMOGENEITIES ON WAVE GROUP 
PROPAGATION IN OUTER SPACE 


Random nonhomogeneities influence wave propagation in a medium: 
they cause amplitude and phase fluctuations of the harmonic wave 
which depend on the frequency, i.e., there occurs a sort of dispere- 
sion which shows up in frequency dependence of the harmonic wave 
amplitude and phase fluctuation. Therefore, during propagation in a 
medium with random nonhomogeneities, the wave group will alter its 
form. Let us assume that weak large-scale nonhomogeneities are pre= 
sent only in the right-hand halfspace (x > 0). From the left-hand 
halfspace (x < 0), there is incident a wave group of arbitrary form: 


+o 


Us(a—ct) =m | C(w,)exp lt, (2 — ct)} da, (43) 


where the harmonic components 2f the wave group have the spectral 
amplitude 


+ 


B (tn) = (Uae —etpoap i= tm, (0 — cf) d(e — th (4H) 
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In accordance with [23 = 25], the electromagnetic field intensity 
can be expressed in terms of fluctuations of the wave level B (k, x) 
and phase S (k, x): 


U (w, 2) = C (a)expl — BF (ay, 4) — 1B (a, 2S (wr, 2) + 8B (os, 2) + 
+ IS (w,, 2) + Ia, (x — cf)l, , (45) 
where C (w, ) 4s defined by (44). After integrating (45) over all 
frequencies, we obtain the general expression for the field inten- 
sity in the wave group: 
U= \ € (1) exp | — B?(,) — 14 (@,) S (0,) + B(@,) + 
+ iS (w,) + iw, (x — et) dm, ( 46 ) 


hence, 


UU = Wewycisyexp {—? (w,) —B* (o,) — t (B (0S (w)—F (3) 5 (oyi+ 


“+ B00) ++ B (org) 3+ E |S (my) — S (ng) -+ 1(, — 0) (x — ct)} doydeag. 


This quantity can be found as a result of statistical averaging of 
the factor 


exp {B (oy) + B (10g) -i- 11S (on) — $ (4)}} 


in the integrand. Considering a normal distribution law of the 
quantities 4 (w,), #(m). S(o), S$(w), 4t 18 not difficult to obtain the 
expression for the average value of this factor 


exp {ap Bom) + Bon) + Bay) BGs) — pS) — yp STs) + 
+ S05 (oy) + (F(a) So) + Blas (or) — Flay S(oy)}} 


then, 


UU’ = 


tos 


\ C (w,)C (om) exp {- + 5*(,) — + 5%) - + B* (14) — 


=F S* (0) + B (04) BC) + S ()S (04) + 618 (04)S (ery) + 


+ B (04) S (ig) + (ny — 4) (2 — et))} dw, diy, (47) 
We see that the mutual correlation and auto-correlation functions of 


the level B (Ww) and phase S (w) fluctuations at different frequencies 
appear in the integrand. If we assume that the correlation 


a4 


coefficient for the refractivity fluctuation has Gaussian form: 


wt 7,3 — ft: 


nem exp (--S*; ’ 
with account for the method of [24, 25], we obtain 
[D+] ~ \\ Compe (oper {— 20 (HE Hh) + 
x DY” ps 


sg Va2V pe = 
Va ah » tap ants Dim Dy _ 


16 Vaz V wine 1 (Di — Day 
= ia aot + Pi De + 
ute ENG (10g) — (4) = 2 (he, — hs) — ( — wei] dur, dug, (48) 
where ig is the mean square of the refractivity fluctuation n, * /54 


n (w,) at the frequency Was D, 2 x/k,L*; the correlation function 


of the field fluctuation is assumed to be Gaussian, with nonhomo- 
geneity dimensions L. We examine the propagation of an arbitrary 
pulse in a nonhomogeneous medium. We expand the exponent in (48) 


into a series in powers of Vi Wy - Yo, and neglect terms of 
third order: 


j04| = K € (i + V4) C (ory + v4) oxp {— G(v, — v4)® + 
+11Q(vt =o) + P (vt — IH dy, | (49) 
where 
G= yeen (+ 7 Di(ny + wang) + [2em.+ ella) > Welty ey 


Q a-ak, +—- yore im +b ilglis ee = + ng 4. tui) +t 


Poskge—t+ Fist, 4- Oya) — Po. 


In these expressions, all the quantities with zero subscript 
relate to the frequency Wos and the differentiation is made with 


respect to w. We shall examine an input signal with amplitude 
spectral density: 


B5 


{a-mP _ 
Cw) = Ce TF alge ®, (50) 
This spectral density distribution selection makes it possible to 
examine all the physical phenomena of interest and, at the same time, 
simplifies the claculatior of (49). 


Substituting (50) into (49) and integrating, we obtain: 


\0*| = 05 exp |— - ; 
[orton alerdreayl 


This expression describes the influence of plasma nonhomogeneities 
on wave group intensity. The nonhomogeneity factor is taken into 
account by the coefficient G. By definition, G@ depends linearly on 
distance. We see from analysis of the pre-exponential factor in 
(51) that the wave group intensity decreases inversely as distance 
because of the defocusing action of the nonhomogeneities. This con=- 
clusion is independent of wave group form and nonhomogeneity varia- 
tion law. The only requirement is the normal nonhomogeneity distri- 
bution law, which, generally speaking, follows from the Lyapunov 
central limit theorem [26]. 


DISPERSION-TYPE SPACE RADIO LINK 


When organizing optimal long-range space communications using 
the proposed method, it is necessary that the spacecraft remain con- 
stantly in the region of maximal signal power. To this end, the 
maximal power region is displaced to follow the spacecraft by selec- 
tion of the radiated signal. Information transmission is accomplished 
using Pulse Code Modulation (PCM) methods. The dispersion-type space 
link consists of the transmitter with variable signal deviation 
frequency, the medium, and the receiver whose passband can be varied 
on command from the Earth and must increase with increase of the 
communications distance. The frequency-modulated transmitter pulse 
with decreasing modulation frequency has the duration t, which in- 
creases with increase of the distance, while the signal center fre- 
quency fy does not change (Figure 3). The rate of change of the 
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modulation frequency B = 2 w/t and 
the center frequency f5 in this 


signal are constant at all times, 
while the spectrum width increases 
linearly with distance. The re- 
ceiver in this radio link must be 
boradband, and special study should 
be devoted to its detailed analysis. 
The filter for the receiver is the 





space medium itself, while the re- Figure 3. Types of radi- 

ated radio signals for 
ceiving element may be a broadband different communication 
maser with controllable band, fol- distances 


lowed by a detector and amplifier. 

After the video amplifier, the signal and noise enter 4 resolver 

with suitable observation criterion: this may be a Zigert-Kotel'nikov 
ideal observer, Neyman=-Pearson minimal risk observer, or a minimax 
observer. For transmission of telemetry information from aboard the 
vehicle to the Earth, the average risk must be minimized by using 
group codes [27]. 


We shall analyze the variation of the signal compression factor 
with distance. With increase of the distance, the signal compression 
factor increases because of the dispersion of outer space, and de- 
creases because of the defocusing influence of *ne plasma nonhomo-= 
geneities. If the compression factor decrease because of nonhomo-= 
geneities is greater than the factor increase, this will hinder 
organization of long-range space communications by the dispersion 
method. We have already seen that the signal compression factor is 
equal to: 


m=2tn fa tT. 


Substituting herein (41) and f, % Af (which 1s valid for large 


modulation factors),we obtain 


162.107 /2.V,2 
mao ——pi-t- . 
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With increase of the distance, the required pulse duration will in- 
crease linearly. In order that the spacecraft be at all times in 
the maximal signal power region, it is necessary that the radiated 
signal duration and frequency deviation increase linearly with 
distance 


F, = AF,z, 


da 
where Fy is the initial frequency; A is a coefficient of propor- 


tionality; then, 
1030-4 -407F, 


pppoe ad 

ie (52) 
i4.e., the compression factor increases as the cube of the distance. 
With increase of the signal compression because of dispersion and de= 
focusing by the nonhomogeneities, the signal power will increase as 
the square of the distance. Since the energy flux density at the 
reception point is inversely proportional to the square of the dis- 
tance, we find that the peak signal power at the reception point 
will be independent of the distance. This conclusion is valid only 
in the case when, with increase of the distance to the spacecraft, 
the transmitted signal duration and frequency deviation increase and 
the deviation magnitude is not too large. 


In the case of compression of a rectangular frequency-modulated 
pulse, there arise adiitional side lobes on which signal energy is 
expended uselessly. Detailed examination of this question shows 
that the signal which is optimal from this viewpoint must have a 
Gaussian envelope. If we select a signal frequency ‘ieviation from 
100 to 500 MHz, the maximal pulse duration for communications with 
Mars amounts to about 34 msec, and with Jupiter — about 0.5 sec. 
The compression factors obtained are 5 + 107 and 109, respectively. 
The dispersion=type space radio link can operate in two regimes: 
compensation for compression factor fluctuation by means of the 
plasma nonhomogeneities, and antenna "noise compensation" by varying 
the antenna zenith angles. The first regime is provided by select- 
ing the operating point uf the radio link at the maximum of the 
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compression factor versus IEC curve. As is known, antenna noise 
increases with increase of the zenith angle, which reduces the 
signal/noise ratio. On the other hand, the IEC also varies with 
zenith angle change. With selection of the radio link operating 
point in the region 2 (see Figure 1), we can obtain an approximately 
constant value of Po/P, for various zenith angles. 


In order to establish contact with extraterrestrial civiliza- 
tions and automatic interstellar stations of the future, it is neces- 
sary to resolve the problem of interstellar communications, which 
can be broken down into three basic aspects [14]: 


1) is it possible to transmit a signal over very long dis- 
tances and, if so, how is this best done; 


2) how can we attract the attention of other civilizations, 
or vice versa; 


3) how and in what "language" can meaningful information be 
interchanged with a completely foreign civilization. 


The most important aspects are the first two, since the large 
distances make information transmission unusually difficult, and 
searching for signals at various frequencies and elevation angles 
1s quite impractical. On the basis of the present level of the 
development of radiophysics, we can expect, in principle, the crea- 
tion during the next 10 = 20 years of antennas with effective area 


10? n°, and recelvers with noise temperature Th = 1° K, Estimates 


for the isotropic case [28] show that such a receiver can record a 
signal from any point of the Universe from a transmitter with power 
PA 1034 erg:sec"2, At the present time, the total amount of energy 
radiated each day by mankind amounts to about 4 + 1029 erg, while 


the Sun's energy output per second is 4 : 1033 erg, which is an 
order of magnitude less than the required transmitter power. 
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The dispersion information transmission method makes possible a 
new approach to the problem of seeking extraterrestrial civilizations 
and providing interstellar communications. It is advisable to or- 
ganize the search as follows. With the aid of phased antennas lo- 
cated at different points of the terrestrial sphere, we can create 
in the distant zone a spherical wave radiated by several frequency= 
modulated transmitters. By varying the transmitted signal duration 
and deviation magnitude, we can "probe" various spherical layers of 
the Galaxy, increasing the search radius in the course of time. In 
order that the signal be compressed into a S-pulse in the region 
of UV-Cetus with frequency deviation from the optical to the RF band, 
a pulse duration of about 5 sec is required, and, in this case, the 


compression factor m 102? - 1023, It is obvious that reception of 
extraterrestrial civilization signals must be accomplished with the 
aid of exceptionally wideband receivers of the Kotel'nikov receiver 
types. Only after estavlishing contact with the civilization is 
dispersion-type communications with directive antennas organized. 

The described technique of dispersion-type reception and transmission 
is known in nature. Thus, certain forms of bats [29) apparently 
utilize the properties of acoustic dispersion of the surrounding air 
space and radiate ultrasonic frequency-modulated signals which are 
compressed into a é=pulise as they propagate. This technique in- 
creases the range of action and resolution capacity of the ultrasonic 
radar. 


In conclusion, we note that the described principles for the 
construction of dispersion=type space radio links provide potential 
noise immunity. 


The authors wish to thank Pref. P. A. Agadzghanov, Dr. Tech. Sei., 
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SPACECRAFT ANTENNA SYSTEM DESIGN 


A. P. Alekseyev, B. A. Prigoda and 
L. I. Skotnikov 


ABSTRACT. We examine methods for grourd 
development of unmanned spacecraft antennas. 
Considering the scope of this question, the ac- 
cent is on the methods and devices for evaluat-= 
ing the characteristics of nondirectional an- 
tennas with ambient condition simulation which 
is as close as possible to the actual conditions. 
Special emphasis is given in the article to the 
questions of measuring directivity characteris- 
tics and studying antenna corona formation. 


The problems of designing antenna systems having nearly optimal 
characteristics are becoming increasingly important in connection 
with the rapid growth of space technology and, particularly, the 
development of unmanned space stations. In the present article, we 
examiiie some aspects of ground=based development of antenna systems 
associated with the methods and equipment for ground-based develop- 
ment with the closest possible simulation of the actual operating 
conditions of spacecraft of the subject types. 


In cont.ast with the conventional antenna systems intended for 
stationary and mobile radio equipment designed to operate under 
ground=based conditions, the development of antenna and feeder sys- 
tems (AFS) for spacecraft requires that several specific conditions 
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be met. Basically, these conditions reduce to the following two 
requirements: isolation from the influence of the Earth and surround= 
ing objects; simulating the ambient medium conditions characteristic 
for the spacecraft during its operation. The validity of the results 
obtained during development of spacecraft AFS, and the quality and 
reliability of operation of the onboard radio equipment and of the 
spacecraft as a whole, will be determined by the degree to which 
these requirements are met. 


Research and development associated with spacecraft antennas 
are usually carried out on special antenna test ranges, equipped with 
stands which permit measuring the spatial directivity characteristics, 
and chambers of various sorts in which the AFS are tested with par- 
tial or complete simulation of the actual ambient conditions. As a 
rule, either full-scale prototypes, the real spacecraft, or their 
models are tested. Antenna development is carried out on models con- 
structed to some definite scale when full-scale operations are im- 
possible because of the large size of the spacecraft, or because of 
the impossibility or avoiding the influence of the Earth and sure 
rounding objects —- for example, in the shortwave band and longwave 
part of the UHF band. 


The capabilities of the modeling technique are limited to the 
cases when the structures of the antennas and of the spacecraft it- 
self are relatively simple, and the external contours can be de= 
scribed by simple surfaces. In this case, modeling is accomplished 
comparatively easily and the electrodynamic properties of the model 
are very close to those provided by the surface of the actual space=- 
craft. However, in the case of complex spacecraft shapes, the crea- 
tion of an electrodynamic model which simulates completely the sure 
face impedance distribution of the actual vehicle is practically 
impossible. Therefore, the study of spacecraft AFS on a full-scale 
prototype is most acceptable. In this case, particular attention 
must be devoted to simulation of those components which are not 
structurally rigid, and may change their position or shape in the 
course of the flight. For example, this applies to evacuated 
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barrier thermal insulation (EBTI), which is a multilayer composition 
of fiberglass cloth and metallized film used to protect certain por- 
tions of the spacecraft surface, and also to the extendible booms, 
hatches, and so on. As @ result of gaseous component release, the 
EBTI may expand in vacuum and its shape may change. Such changes 
can lead to redistribution of the surface currents excited by the 
antennas located on the spacecraft, which in turn causes distortion 
of the antenna pattern and disruption of communication between the 
Spacecraft and tne Earth. This effect is particularly undesirable 
if there are on board Doppler radio equipments, which are very sensi- 
tive to parasitic fluctuation of the signal received by the receiver 
due to the presence of surface segments with varying impedance char- 
acteristics, which may lead to the appearance of a false signal and 
interruption in system operation. 


Thus, the first and foremost requirement for AFS development 
is complete simulation of the external electrodynamic properties of 
the prototype, and assurance of the required stiffness of all the 
prototype components which participate in the formation of the basic 
parameters and characteristics of the spacecraft antenna system. In 
the following, »s present some examples of how spacecraft AFS devel- 
opment is carried out. Because of the limited space available, the 
discussicn is restricted to only certain AFS characteristics: spa- 
tial directivity pattern and electrical strength of the antennas 
under nearly realistic conditions. The study is made for broad- 
beam antennas, for which the spacecraft hull is an active element 
which influences these characteristics. 


STUDY OF SPATIAL DIRECTIVITY CHARACTERISTICS 


A large number of test facilities intended for the measurement 
of spacecraft antenna directivity characteristice have been described 
in the literature. In spite of their variety, they can be broken 
down into quite clear classifications. First, they can be sub-= 
divided on the basis of the degree of automation of the test object 
displacement process and recording and protessing of the results 
into fully automated facilities, partially automated facilities, 
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and manually controlled facilities. Second, they can be classified 
on the basis of structural characteristics: with stationary measuring 
radiator, and with mobile radiator which displaces relative to the 
test vehicle center of mass. Other classification criteria can also 
be used, for example, based on the nature of the radiation source 
used: ground, airborne, or galactic sources. 


In any case, regardless of the construction and nature of the 
test stand facilities used, they must be subject to the basic require- 
ments on which the measurement accuracy and validity of the results 
obtained depend. These requirements amount to the fact that the dis- 
tance between the test and auxiliary antennas must satisfy the dis- 
tant-zone conditions: the signals reflected from the ground and 
foreign objects must not affect the measurement results, and the 
mutual coupling between the test and auxiliary antennas must be negli- 
gibly small. As is known from the literature, the first condition 
is satisfied when the distance Rain between the spacecraft being 


tested and the auxiliary antenna: 

Ran = Fe. (1) 
where A is the wavelength, D is the maximal antenna dimension; k is 
the coefficient defining the measurement error. If the tolerable 
error AE in determining the field intensity is no more than 1%, 
this expression takes the specific form: 


Ron ==. (2) 


The measurement accuracy will increase with increase of the distance 
between the test and auxiliary antennas. This is shown in Figure l. 


In order to reduce the distortions caused by interference be-= 
tween the direct rays and those (secondary) rays reflected from the 
ground and local object, the area between the antennas must be open, 
the antennas should be located on towers, special reflecting panels 
or mats made from radiation absorbing material should be installed 
between the antennas. The antenna height above ground and the 
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minimal acceptable distances from the nearest interfering objects 
which create secondary waves are determined from the formula: 


hb Ain 
i c= y —, 
T Dad (3) 
where Dad is the linear dimension of the auxiliary antenna. 


The antenna test stands are valibrated 


af, % 

to exclude their influence on the test an- 100 
tenna characteristics. + 

If the auxiliary antenna dimensions 40 
are less than those of the test antenna, 
i.e., the auxiliary antenna has a broad Lars G4 40 404 
directivity pattern, then the Condition Figure 1. Measure- 

. ment error versus 
(3) becomes: distance between 
HW oe : antennas 
rad 


where D is the linear dimension of the test antenna. We see from 
what we have said that it is best to work with highly directional 
auxiliary antennas. 


Figure 2 shows schematics of test stands used to study the 
directivity characteristics of spacecraft antennas. In the variant 
of Figure 2a, the test vehicle is mounted on a horizontal pivoting 
platform, which provides azimuthal rotation relative to the vehicle 
center of mass with any orientation of the vehicle axes in space. 
The variant of Figure 2b differs in that the vehicle is mounted on 
a system of flexible cables, rather than on a rigid base. By choice 
of the cable support system and using a drive located in the cable 
support system, it 18 possible to make measurements of the spatial 
directivity characteristics. This scheme is particularly convenient 
when conducting impedance measurements, and also in studies of the 
polarization characteristics, since it permits altering the vehicle 
height above ground in wide limits, and thereby provides optimal 
isolation from the ground. In order to achieve optimal isolation, 
it 18 best that the towers and area between them be covered by 
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Figure 2. Antenna directivity characteristic meas- 
urement stands: 


1 — test vehicle; 2 — auxiliary antenna; 3 — 
drive; 4 — radiation absorbing material; 5 — re- 
flecting panels; 6 — cable support 


panels and mats made from radiation absorbing material. The test 
stands in the Figure 2a and 2b variants are used in studying space- 
eraft antenna systems intended for operation under free-space 
conditions. 


Special stands are used to conduct studies of the antenna sys- 
tems of lander-type spacecraft which are intended for operation from 
the surface of planets. An example of such a stand is shown in 
Figure 2c. The landing vehicle 1s mounted on a surface simulating 
the assumed planetary soil model in the landing region. The auxtl- 
iary antenna displaces along a circular arc at a distance A>»D*/\, 
from the lander's antennas. The signal level is recorded on the 
tape of a recorder which displaces synchronously with the measuring 
antenna movement. Studies of the directivity characteristics of the 
landing vehicles of the unmanned Luna, Mars, and other space stations 
were made in this way. 


STUDY OF SPACECRAFT ANTENNA ELECTRICAL STRENGTH 


One of the specific ‘aracteristics of spacecraft antenna 
operation is that of functaoning under conditions of a highly rare=- 
fied gaseous medium. Observations made during flights at altitudes 
on the order of severel tens of kilometers above the Earth have 
shown that antenna breakdown may occur even with low input power, 
measuring only a few Watts. Arcing occurs in those cases when the 
electrical field intensity reaches a definite level above the criti-e 
cal value. The use of dielectric coatings is advisable in order to 
reduce the field intensity level and reduce the possibility of arc- 
ing. The field intensity in the dielectric decreases in proportion 
to the magnitude of the relative dielectric permeability of the 
coating material. 


If these measures are taken, breakdown of the gaseous medium 
surrounding the spacecraft will exert the primary influence on the 
electrical strength. Breakdown may show up either in the form of 
a spark discharge between individual parts of the structure, or as 
a corona discharge having the form of a plume which develops on some 
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part of the antenna. The nature and 
intensity of high-frequency breakdown 
are determined by the pressure and 
temperature of the gas, the presence 
of fonizing radiations, the frequency, 
average power, and duration of the . 
high-frequency signal, the dimensions Ge MRE anh’ He 
and shape of the antennas, and so on. Figure 3. Breakdown power 
As a rule, in studying the discharge versus pressure 
mechanism and developing practical 

recommendations on compensating for the discharge, the primary atten- 
tion is devoted to experimental rather than theoretical studies, 
since the latter constitute a quite formidable task. 





Breakdown power 


The minimal electrical strength of an air medium is obtained 
when the antenna electromagnetic field circular frequency w is equal 
to the frequency v of electron collision with neutral molecules, 
which is expressed in terms of the gas pressure p in mm Hg as 
follows: 


v me 5,3-10%p. 


Consequently, we can write for w = v: 
pa = 36, 


where 4 1s the wavelength in cm. Figure 3 shows the characteristic 
curve of breakdown power dependence on pressure for signal frequency 
420 MHz. We see quite clearly the pressure region where breakdown 
shows up for minimal values of the power supplied to the antenna. 
With increase or decrease of the working frequency, this region will 
displace along the pressure scale to the right or left, respectively. 


Special chambers are used to conduct such studies. Figure 4 
shows a schematic uf a chamber in which experimental studies were 
made of the whip antennas of the Mars 3 unmanned space station 
landing vehicle. The hemispherical fiberglass dome is equipped with 
viewing portholes. The vacuum pump system provides a vacuum of 
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Figure 4. Radiotransparent vacuum chamber for developing 
antenna and feeder systems: 


1 — vacuum chamber; 2 — test article; 3 —- lonizger lamp; 
4 — vacuum evacuation system; 5 — control console; 
6 — radiation absorbing material 


1074 mm Hg, which is quite adequate for simulating the pressure both 


in the descent trajectory and after the landing vehi-le touches down 
on the surface of the planet Mars. A special gas supply system 

makes it possible to create in the chamber cavity a gas composition 
simulating various Mars atmosphere models. Simulation of the solar 
radiation which creates the initial ionization was accomplished by 
tubes of the PRK (mobile X-ray univ) type located inside the dome. 

By altering the pressure, gas composition, and level of the power 
supplied to the antennas, we can evaluate quite completely the pos- 
sibility of the occurrence of high-frequency breakdown on the antenna 
under foreign planetary atmosphere conditions. 


Studies made in recent years have shown that gas composition has 
very little effect on breakdown level. The primary influence is 
that of pressure. The formation of ionization which transitions 
into corona is a very undesirable phenomenon which increases the 
high-frequency iosses, changes the antenna input resistance, ine 
creases the standing wave ratio (SWR) in the channel, distorts the 
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directivity pattern, leads to partial or complete shielding of the 
antenna, failure of the onboard transmitter, and so on. The photoe 
graphs in Figure 5 show how the corona formation distribution pattern 
on a whip antenna changes with change of the pressure of the medium 
in the chamber. Beginning at the end of the whip, at a pressure of 
0.073 mm Hg, the corona formation extends over the entire length 

of the whip, reaching maximal intensity in the 1.0 mm Hg region for 
frequencies of ~ 400 = 500 MHz. 





Figure 5. Corona formation on whip antenna at 
different pressures 


Ic is obvious that the critical pressures in the atmospheres 
of the different planets do not show up at the same heights as for 
terrestrial conditions. Therefore, the problem of corona formation 
during entry into the atmospheres of these planets will arise in 
different parts of the entry trajectory. In this connection, it is 
important to select the transmitter frequenctes and power of the 
capsules and landing vehicles intended for descent to the planet 
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surface, since the critical pressure at which breakdown begins is 
very sensitive to frequency. In the Earth's atmosphere, the point 
corresponding to the minimal breakdown voltage value appears in the 
region where the product of the pressure (in mm Hg) and wavelength 
(in em) is between 10 and 40a. 


One of the problens in spacecraft antenna design is that of 
gas release from the materials forming the spacecraft shell, antenna 
insulating washers, disks, and so on. An experiment has been de- 
scribed in the literature in which a slot antenna filled with poly- 
styrene foam operated satisfactorily to pressures of 0.046 mm Hg 
with 200 W power. During repeat tests, after three-day exposure 
in a vacuum chamber, breakdown was observed at 80 W power. Gas ree 
lease from the materials in high vacuum complicates this problem. 
The high temperature conditions on the surface of the Mars and Venus 
probes and landers make the problem associated with outgassing 
more critical. 


INFLUENCE OF SURFACE HEATING AND PRESSURE 


The influence of temperature and pressure at the spacecraft 
surface is important not only from the viewpoint of the possibility 
of distortion of the antenna characteristics as a result of corona 
formation. Heating of the spacecraft surface may lead to change of 
the properties of the dielectrics used to coat the antennas, change 
of the geometry of the antenna itself, and so on. In this case, 
there may be marked changes of characteristics, such as the antenna 
input impedance, which influences directly the channel SWR and radi- 
ated power level, and distribution of the high-frequency currents 
over the vehicle surface, which basically determines the spatial 
directivity pattern. It is obviously best to study the spacecraft 
or their individual components which are subject to heating or cool- 
ing, together with the antennas and high-frequency channel compon- 
ents which are located in these areas. Two measurement facility 
schemes are used for this purpose. In the first case, the heat 
source and radio emission (radio reception) source are separated in 
space. Figure 6a shows a schematic of such a facility. The exhaust 


103 





fo) recorder - 


rn ae 


Figure 6. Schemes for heat testing lander antennas 
installed beneath thermal insulation layer: 


1 — jet engine; 2 — lander thermal protection coat- 
ing; 3 ~— test antenna; 4 — auxiliary antenna; 5 — 
mounting booms 


of a jet engine is used as the heat source. The radio emission 

source is located at some distance from the tent antenna. During 

the heating process, a record is made of the signal received by the L867 
test antenna, and the change of this signal, and also the SWR varia- 
tion, is used to evaluate the degree of effect of heating on the 


104 


antenna. Anot*ser facility scheme (Figure 6b) 1s also possible, in 
which the Sur. is used as the heat source and radio emission source, 
Tnis is accomplished witr the aid of a parabolic mirror c.ncentritor 
at the focus of which the test vehicie with antenna is located. 
Heating 48 accomplished by the radiaiut energy of the solar spectrum 
and the Sun's electromagnetic flux in the RF band is used as the 
source of information on the nature of the change of the radiotechni- 
cal properties of the antenna. The advantage of this facility con- 
figuration is that the Sun's RF emission spectrum is practically con- 
tinuous, and check of the antenna parameters can be made simultane- 
ously over a very wide frequency range. 


In conelusion, we shall discuss briefly the facility used to 
evaluate the operation of certain antenna units under high pressure 
conditions which occur, specifically, on the surface of the planet 
Venus. Such a facility is shown in Figure 7. It was used to test 
the Venus landers and, specifically, to develop the mechanism for 
firing the ejectable antenna of the landing module of the Venus 8 
unmanned space station, which is shown in this same figure. The 
ejectable antenna was ejected at a gas pressure on the order of 100 
atm, and heating of the ejection mechanism to 400 = 500° C. The ef- 
fectiveness and reliability of the ejection mechanism was evaluated 
on the basis of the cage ejection height. The spacecraft antenna 
and feeder system can be equipped with various mechanisms, with the 
aid of which rotation, opening, and other displacements of the an- 
tenna feed system are accomplished. In these cases, AFS reliability 
must be verified on the Earth in order to increase the reliability of 
AFS mechanism operation under actual conditions. To this end, the 
spacecraft AFS mechanisms are tested in vacuum chambers anil in ther- 
mostats, with heating or cooling of the mechanisms and the AFS, 


Examination of the techniques for developing spacecraft AFS 


makes it possible to draw the following conelusions: test stand 
facilities intended for development of the spatial directivity 
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characteristics of spacecraft antennas should provide maximal 1s0=- 
lation from the Earth, and introduce minimal distortions into the 
measurements; the vacuum, thermal, and other special test stands 
must simulate, as far as possible, the actual conditions of space- 
craft AFS operation: temperature, pressure, gas composition, Sun's 
radiation spectrum, and so on. The successful flights of the Luna, 
Mars, Venus spacecraft confirmed the effectiveness of the methods 
used for ground development of spacecraft AFS. 
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LOW-SILHOVETTE SPACECRAFT ANTENNA SYSTEMS 


B. A. Prigoda 


ABSTRACT. We examine ways to reduce the 
Silnouette and overall dimensions of spacecraft 
antennas. 


The modern spacecraft is connected with the external medium by /67 
& system of sensors, optico-mechanical sensitive elements, antennas, 
and other special devices, each of which has a definite zone of 
action (spatial angle operating sector). In solving the general 
problem of studying outer space and the properties of the atmos- 
pheres and surfaces of the celestial bodies, we often must make com- 
promises to provide the maximal capabilities for operation of a cer= 
tain group of onboard devices and artificially narrow the zone of 
action of the remaining instruments and devices, which are of less /68 
4mportance in the particular experiment. This approach is definitely 
undesirable, since it leads to a situation in which there is reduc-= 
tion of the potential capabilities of the scientific, telemetry, 
radio command and other systems, located both aboard the spacecraft 
and along the communication line between the spacecraft and the 
Earth. This frequently takes place because the zones of action of 
the various sensors, antennas, and other sensitive elements overlap 
spatially. When configuring these elements aboard the spacecraft 
within the limits of the specified overall weight and definite 
limitations on the dimensions, it is not always possible te separate 
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these instruments in space as they should be and mutual shadowing 
takes place, which leads to a situation in which — if there are 
aboard, for example, optical orientation sensors -—~ certain parts of 
the active scan zone are cut off, and if the celestial body falls in 
these regions, interruptions will take place in astro-orientation 
system operation. This leads to lengthening of the orientation 
seances and, consequently, additional expenditure of the service 

life of the onboard equipment and power supplies, which in turn leads 
to increase of the overall weight, and so on. 


In the present article, we examine some questions associated 
with reducing spacecraft antenna system dimensions, or, as ve usually 
say in aviation, reduction of the antenna silhouette or design of 
low-silhouette antennas. 


Up to ten or more antennas, each of which performs a definite 
task in regard to information reception or transmission, are in- 
stalled aboard the modern spacecraft. The following approaches can 
be used in resolving the problem of reducing the antenna system over-= 
all silhouette and dimensions: 


1) creation of frequency-independent antennas [1], 1.e., com- 
bining the functions of several antennas operating in different 
bands and having analogous electrical characteristics into a single 
antenna unit; 


2) creation of controllable antennas [2, 3], 1.e., combining 
the function of several antennas with different electrical charac-= 
teristics into a single antenna unit; 


3) optimal configuration of the onboard antenna and feeder 
unit with scanning and connection of the antenna with optimal char-= 
acteristics to a common feeder channel at each specific instant of 
time (4); 


109 


4) use of time diversity between the information reception and 
transmission seances in the case of overlap of the sensor spatial 
operating zones and antenna directivity patterns in the presence of 
their partial mutual shadowing; 


5) use of magnetodielectric coatings which make it possible 
to obtain the required antenna characteristics with reduction of 
their overall dimensions [5]. 


The first technique has been quite thoroughiv studied at the 
present time, and is widely used in practice. Thus, in spacecraft 
engineering various log-periodic structures (planar, cylindrical, 
hemispherical, and conical) are quite frequently and successfully 
used, which make it possible to obtain the same electrical character- 
istics in a wide frequency tand with overlap on the order of 10 or 
more, 1.e., f/f, > 10, where f, is the lower edge of the frequency 


band, ry is the upper limit of the frequency band. 


At the present time, the second technique is gradually beginning 
to be introduced into the engineering of spacecrart antenna units. 
This approach involves changing the configuratic- of the antenna or 
certain of its mechanical or electrical characteristics, so as to 
vary the required antenna parameters, depending on the specific task /69 
imposed on the antenna at the given instant of time. The antenna 
System of the Venera unmanned space station ean serve as an example. 
At certain times during the transfer trajectory, the station oper- 
ates in the single-axis orientation regime, and rotates about the 
direction of the sun, close to or coinciding with the axis of prin-e 
cipal moment of inertia. In order to provide communication with a 
ground station, it is obviously necessary, in this case, to have 
either an antenna unit which tracks the direation to the Earth or an 
antenna with radiation pattern of funnel-shaped form with maximum 
oriented toward the Earth. Since in the transfer trajectory, the 
Sunevehicle-Earth (SVE) angle varies following a definite law, it is 
obvious thet the Earth will gradually leave the zone of action of 
the antenna pattern maximum. In this case, in order to ensure 


110 


optimal communication along tre 
entire transfer trajectory, it 


is necessary to have several A sun 
antennas with funnel-shaped ge : 
conical radiation pattern, but 7 

with spatially separated pat- Vehicle ' 7 

tern maxima covering the zone SVE 

of possible SVE angles for . 


the given vehicle. 


The figure shows the curve 
of SVE angle variation and the 
set of radiation patterns cover- 
ing the SVE angle variation 


zoue (Gove min ‘ SSvE max?* As 
the antennas which provide pate 
terns of funnel=-shaped form 
with elliptical polarization 

of the antenna radiation field, 


it is customary to use for P f SVE : aed 4 

urve 0 angle variation an 
spacecraft log-periodic two- set of radiation patterns cover- 
and four-turn conical and hemi- ing the SVE angle zone 


spherical antennas for which 270 
the width and orientation of the radiation pattern maximum is deter- 
mined basically by a single parameter -—— the spiral wrap angle. The 
radiations patterns shown in Figure 1 are provided by three four- 
turn conical log-spiral antennas with spiral wrap angles equal to 
45°, 52°, and 65°, respectively. For the given antenna class, the 
problem of providing with a single antenna all three of the radia- 
tion patterns mentioned above reduces in practice to gradual variae 
tion of the wrap angle from 45° to 65°. A similar effect can be 
achieved by locating passive elements of the metallic disk or ring 
type at the apex of the spiral [3]. 
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The third and fourth techniques have much in common with one 
another, since they are based on alternate use of individual antenna 
elements with utilization of retractable booms, deploying mechanisms 
and other mechanisms which make it possible to reduce the dimensions 
of the system as a whole. In acdition to the necessity for intro- 
ducing certain mechanisms into the system, another drawback of these 
techniques is the presence of electronic analyzing devices, which 
are bound to have some influence on the reliability of the entire 
system. Studies are being carried out at the present time to opti- 
mize the onbaord system and maximize is reliability. 


Reducing the size of the antenna radiators proper, even with 
account for the possibility of using the aforementioned techniques, 
is still advisable. Since simple reduction of the dimensions of 
any antenna leads directly to distortion of such parameters as the 
input impedance Z2,,, directive gain (DG), radiation pattern shape, 


polarization properties of the antenna, and others, it is necessary 
to find techniques for reducing the geometric dimensions of the an-= 
tenna such that its electrical characteristics will remain unchanged. 
One such technique is to locate the antenna in a magnetic-dielectric 
medium with parameters which differ from those of free space. For 
example, the aboveementioned multi-turn log-spiral concical antennas 
are remarkable in that their upper and lower working frequency band 
limits are determined by the truncated cone apex and base diameters. 
Since the wave propagation phase velocity, which defines its length 
in the given medium, depends on the properties ce and uw of the mediun, 
we can asaume that, other conditions being the same, the wavelength 
in the more dense medium (e > Eg» u > My) will differ from the wave- 


length in free space and, consequently, the effective dimensions of 
the radiator located in this medium will also differ. Studies which 
have been made show that the dimensions (a) of the active region of 
antennas coated with a layer of material with parameters differing 
from those of free space (€ = €), u A uy), and of the uncoated an- 


tenna (located in free space) are connected by the expression! 
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where e€ is the dielectric permeability of the coating material; uw is 
the magnetic permeability of the coating material; Ey» Hy are the 


a(#.3) 


= (05, fen 


(1) 


corresponding parameters of free space; 285 is the spiral apex angle; 


~ 18 the spiral wrap angle. Expression (1), obtained for conical 

spirals, can also be used successfully for bifilar spirals if we set 

6 = 0. For example, for a bifilar spiral with wrap angle » = 6.5°, 

coated with a ferrite layer having « = 3.77 and u = 2.2, the theoreti- /71 
cal value of the size reduction coefficient is n = 0.58, in accord- 

ance with the formula shown above. This value of n is in good agree- 

ment with the values obtained in the laboratory in studying specimens 

of such antennas. Formula (1) shows that greater effectiveness can 

be achieved by using materials with large value of e« than materials 

with large u for the coating. 


The effective reduction of the antenna dimensions is determined 
not only by the properties of the coating material, but also by the 
coating thickness, its homogeneity, and certain properties of the 
specific antenna. The effectiveness of coating application will be 
different for different antenna types. The use of coatings with 
large value of e€ or yw involves some undesirable effects, such as, for 
example, distortion of the input impedance, reduction of the channel 
traveling wave ratio (TWR), reduction of antenna efficiency, and so 
on, Which must be considered during antenna development and compen= 
sated for in some definite fashion. 
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HIGH-SENSITIVITY 3.5-cm MQDULATION-TYPE RADIOMETER 


A. Ye. Andriyevskiy, A. G@. Gorshkov, V. V. Danilov, 
Vv. K. Konnikova, A. S. Lobarev, V. @. Mirovskiy, 
Vv. V. Nikitin, V. I. Portman, Ye. A. Spangenberg, 

I. A. Strukov, N. Z. Shvarts and V. S. Yetkin 


ABSTRACT. A brief description of the circuit, 
characteristics of the functional elements, and 
parameters of a highly sensitive modulation-type 
direct-amplification radiometer operating in the 
8000 - 9000 MHz band are presented. The radio- 
meter was used together with the RT-22 antenna 
(Crimean Astrophysical Observatory of the Academy 
of Sciences of the USSR, located at Simenz). The 
sensitivity of the radiometric receiver was 
AT ¢ 0.02° K. 


A highesensitivity modulation-type radiometer, operating at the 
3.5 em (8.55 GHz) wavelength, was developed in 1967 = 1969 at the 
Institute of Space Studies of the Academy of Sciences of the USSR, 
the Moscow Pedagogical Institute im. V. I. Lenin, and the State 
Astronomical Institute im. G. K. Shternberg of Moscow State Univer- 
sity, and has been used for radioastronomical observations on the 
antenna of the RT+22 Crimean Astrophysical Observatory (located at 
Simenz, Crimea). Fxtensive observational astronomical material has 
been accumulated, and the required information on the technical and 
operational parameters of the receiver has been obtained during the 
time of radiometer operation on the antenna. The results of the 


115 


ari 


first stage oi the observations were published in the Astronomical 
Circular of the Academy of Sciences of the USSR [1], are now being 
yublished in the Astronomical Journal of the Academy of Sciences of 
the USSR and in News of the CAO of the USSR, and will be published 
in the future as the data obtained are reduced. The objective of 
the present article is a brief technical survey of the components 
of the radiometric receiver. 


The modulation radiometer is constructed using the direct am- 
plification receiver scheme with v 1000 MHz microwave circuit pass- 
band, and has sensitivity AT ~ 0.02° K, with storage time constant 
t = 1 sec. Lobe modulation with lobe switching frequency F = 1000 
Hz is used in the radiometer. A simplified block diagram of the 
radiometer's microwave circuit is shown in Figure l. 


Radiometer input circuits. The radiometer feeds provide opera- 
tion in the lobe modulation regime and permit measuring the degree 
of circular polarization of the received radiation. Structurally, 
they are built in the form of two joined conical horns with circular 
waveguides (polarizers) and smooth transitions from the circular to 
rectangular section (Figure 2). The feed horns are located in the 
horizontal plane passing through the optical axis of the antenna 
system, and are displaced from the axis by + 9'. The phase center 
of the feed horns is aligned with the focal plane of the antenna 
paraboloid. We note that for antennas with small aperture, the in- 
fluence of attenuation fluctuations in the atmosphere on the radia- 
tion flux measurement accuracy decreases only slightly when using 
lobe modulation. For the RT=22 antenna, the sensitivity gain amounts 
to only a factor of four. But when using two horns, we obtain a 
characteristic output signal form which facilitates separating the 
source signal from noise when analyzing the output signal recordings. 


The conical feed horns provide 15 dB irradiation intensity drop 
at the edge of the mirror. In this case, the antenna utilization 
factor is 45%; the first radiation pattern (RP) sidelobe has an ine 
tensity of less than 18 dB. The detailed nature of the side lobes 
was not investigated. 
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Figure 1. Block diagram of radiometer microwave circuit: 


B — balancer; M — modulator; C — circulator; ML— matched 
load; F — filter; PA — para:ietric amplifier; TA — tunnel 
amplifier; WS — waveguide short; D — detector; PAB —= power 
amplifier block; RFG — reference frequency generator; ¢ — 
phase shifter; AT — attenuator; G — gate; PG — pumping 
generator; TB — thermobattery; TS — temperature sensor 


The device fcr measuring 
the degree of incoming radia- 
tion circular polarization 
consists of polystyrene 
quarter-wave dielectric plates 
{e = 2.5), mounted in a circu= 
lar waveguide, with one horn 
providing reception of radta- Figure 2. Schematic of antenna 
tion with right-hand circular peed noth (poterzaer: 
polarization, and the other 1 Sekt ge anceeh OR 9a am 
horn providing reception of 
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radiation with left-hand circular polarization. The difference in 
the record of the signal received by the first and second horns 
yields the degree of circular polarization and its sign. A drawback 
of the described polarizer is the difficulty in matching the polari- 
zer with the receiver and feed circuit; as a result, it was not 
possible to obtain a feed SWR better than 1.7 in the band % 1000 MHz. 
However, the quite high value of the feed SWR did not lead to marked 
deterioration of the fluctuation sensitivity of the radiometer, 
since the required isolation (*% 40 dB) was provided at the input of' 
the first parametric amplifier. This feed system, together with the 
RT-22 antenna, made it possible to realize a radiation pattern of 
(6.3' + 0.1) x (6.6" + 0.1), with 18' angle between the main lobes. 


Modulator. The modulator is a ferrite waveguide switch devel- 
oped for this purpose in which change of the direction of 2irculation 
of the Hig linearly polarized electromagnetic wave is accomplished 
in a Yecirculator. A cylindrical insert made from type ZSCh-15 
ferrite is used in the circulator. In order to reduce the eddy 
current losses in the circulator body and electromagnet magnetic 
structure (which is the basic problem in the design of a switch of 
this type if the switching frequency exceeds 150 - 200 Hz), the 
waveguide tee was fabricated by the galvanoplastic method, and has 
waveguide wall thickness no more than 0.05 mm; the magnetic structure 
is fabricated from type 2000NM1 low-frequency ferrite, the solenoids 
are made in the form of two allemachined flat springs having 20 turns 
each. The entire modulator structure is pressed from plastic, and 
does not exceed in size the dimensions of a 3=cm band circulator. A 
schematic of the switch construction is shown in Figure 3, and its 
electrical characteristics in the working frequency band are shown 
in Figure 4. A current amplifier is used to control the ferrite 
switch; the circuit provides relay commutation to block the switeh 
in either of the extreme positions (control of microwave ciroult 
amplitude-frequency characteristic, and so on). 


Salibrator. The modulation radiometer with three-dimensional 
scanning of the antenna radiation pattern has a dual-channel micro= 
wave system ahead of the antenna switeh (modulator). Yor 
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technological reasons, 
the components of the 
dual cnannel system cane 
not be made absolutely 
identical, and usually 
have attenuation scatter 
in the range of 0.1 = 

0.2 dB. Due to this, 

the radiometer sensiti- 
vity decreases slightly 
(by 3 - 6%); however, 

an additional unbalance 
Signal appears which 1s 
commensurate in magnitude 
with the antenna tempera- 
ture increase during use- 
ful signal reception (up 
to 10 - 15°) from power- 
ful cosmic RF sources. 

A special component is 
introduced into the 
microwave circuit to come 
pensate for the loss un- 
balance in the two radio- 
meter input channels. 
This component consists 
of two parallel rectangu- 
lar waveguide segments 
into which an absorbing 
plate is inserted with 
the aid of an electric 
motor; this plate equal- 
izes the losses in the 
microwave circuits. ‘The 
maximum attenuation which 





Figure 3. Modulator: 


1 — magnetic structure; 2 — waveguide; 


tee; 3 —— microwave ferrite; 4 — 
matcher 
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Figure 4. Modulator electrical char- 
acteristics: 


1 — direct losses; amg isolation; 
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can be introduced into either of the input circuit arms is %~ 0.2 GB. 
On the common axis of rotation of the electric motor and the absorbe 
ing plate, there is mounted a type PTP=2 precision potentiometer 
from which an emf proportional to the attenuation introduced into 
one of the arms is taken. This emf 1s recorded by a pointer-type 
instrument on the radiometer control console, and makes it possible 
to calibrate the entire radiometer using an artificially introduced 
noise signal. An estimate of the accuracy of this calibration 
technique was made using the observational data, and the accuracy 

is about 3%. 


The high-frequency amplifier block consists of two parametric 
and two tunnel amplifiers with wide passband 1000 MHz and overall 
gain G@ = 50 aB. An estimate of the required UHF gain was made fol- 
lowing the calculation made by Korol'kov [2]. The sensitivity of 
the modulation type radiometer is defined by the expression: 


at wa rae. 
aT V2 Tay 7 tt Guia? (1) 


where Th is the receiver input equivalent noise temperature, Af is 


the UHF circuit passband; AF is the LF circuit passband (AF = 1/4 7, 
where t is the storage time constant); To is the video detector 


temperature; @ is the UHF circuit gain; k is the Boltzmann constant; 
M is the video detector quality factor. The first term of the radi- 
cand defines the fluctuations caused by the detected HF noise of all 
the receiver elements preceding the crystal detector; the second 
term defines the contribution to the fluctuations introduced by this 
detector and the subsequent LF stages. In a well constructed radio- 
metric receiver, the condition must be met that the output noise 
introduced by the detector and LF section should be no more than 

10% of the noise introduced from the microwave circuit: 


ToAF 
20h ar = 10!4+ ater) =) 
or 
100 4 fT 
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In this case, the sensitivity of the compensation type radiometric 


receiver is determined basically by the first term of the radicand 76 
in (1): 
ar =, Be (4) 
V2 Var" 


According to [2, 9, 10], the sensitivity of the modulation type 
radiometric receiver with degenerate parametric amplifiers at the 
input, rectangular modulation, and sinusoidal demodulation is 


Pee ms 
AT= Vivar’ (5) 


Repeated measurements of the sensitivity of the developed receiver 
under laboratory conditions and on the antenna (using calibration 
sources) showed good agreement of the experimental results with cal- 
culation using (5). The self-noise of the receiver without antenna 
did not exceed Th = 150° K, and with individual parametric diodes 


decreased to Th = 130° K. In connection with this low noise UHF level, 
we developed a special detector section with quality factor M = 


150 = 200 wl? | According to (3), with these basic radiometer ele- 
ment parameters, the optimal UHF gain should be G = 109 (50 dB). 
Further increase of the radiometer UHF gain is not advisable for the 
following reasons: 1) with increase of the UHF gain, the stability 
of regenerative UHF stage operation decreases; 2) the probability 
arises of saturation of the last tunnel amplifier stage; for Th ant 
30° K and gain G = 38 dB in the first three UHF stages, power P = 

k (T, + T,) GAf e 2. 1078 Wis applied to the input of the last 


tunnel amplifier. 


The radiometer video detector maintains squareness of its char= 


acteristics reliably up to input signals not exceeding 1076 W. 
High UHF gain will reduce the dynamic range of the radiometer. 
Therefore, when using the radiometer for observations of quite 
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sources (with antenna temperature Tv, ~ 500° K, provision must be made 


in such UHF schemes for disconnection of the last UHF stage (second 
tunnel amplifier). 


Parametric amplifier block. Degenerate parametric amplifiers 
with gain G@ = 14 dB in the band Af = 1 GHz are used as the first two 
UHF stages. The parametric amplifier passband lies in the range 
8050 = 9050 MHz. The gain nonuniformity in the operating frequency 
range of each stage amounts to 0.2 - 0.3 dB. In the amplifiers, we 
used specially developed parametric diodes having short time con= 
stant, which makes it possible to realize quite low amplifier self- 
noise temperature. Thus, on some diode specimens, the parametric 
amplifier self-noise, measured at the input of the circulator and 
filter 1 (see Figure 1), amounts to T, v 90 - 100° K. The circu- 


lators used had the required pussband with isolation of 20 = 25 dB 

in each arm, and direct losses cf % 0.3 dB in the arm. We were able 
to select circulator spccimens with low direct losses % 0.18 = 0.2 4B 
in the arm for the radiometer input and for PAs which made it pos= 


sible to reduce somewhat the receiver self-noise (by 7 = 15° K). 
The circulator Ch» providing the required isolation between the modu- 


lator and PA,, was connected at the input (after the modulator M). 


This makes it possible to reduce the parasitic signal caused by in- 
terference of the amplifier self-noises with change of the modulator 
output impedance. 


In order to prevent penetration of the pumping power from the 
PA to the modulator, which in turn may lead to the appearance of a 
parasitic signal, a LF filter (F,) is provided at the input of PA, 


to reduce the pumping power by 40 = 50 dB. The losses at the signal 
frequency in the filter are quite small, and do not exceed 0.2 dB. 
The filter is made with a two-sided “waffle” structure, and low 
losses are realized by maintaining high surface precision and surface 
finish. Prior to silvering, the filter surfaces were subjected to 
electrogalvanic polishing, which made it possible to improve 
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considerably the surface finish after machining. The filters Fy and 
F, were fabricated similarly. The filter Fy isolates the Pi stages 
from the pumping power of the neighboring PA. The filter F pro=- 


tects the tunnel amplifier block against the possibility of over- 
loading by the pumping power which does leak through. The overall 
losses in the signal circuit from the input of circulator Cy to the 


input of PA) (i.e., the arm C,+ F, + the arm Co) amount to about 
0.45 @B, 


Pumping circuit. The pumping generator is a klystron located 
in an 011 bath to increase the heat rejection of the instrument. 
The pumping power from the klystron flows through the gate a, to the 


power divider, based on a dual T-bridge. The bridge was specially 
tuned to the pumping frequency with low active and reactive losses, 
and for optimal isolation of the arms in this case. As a result, 
VSWR parameters in the pumping generator arm equai to 1.5 were ob-= 
tained; the isolation between the pumping arms of each PA is 14 dB; 
the pumping power flows along the two paths through the gate Go,» 


waveguide short WS,» and attenuator AT, (similarly, through Q35 WSo,5 
AT.) to the parametric amplifiers PA, and PA,. The gates G, and Gs 


isolate the PA stages at the pumping frequency, which simplifies 
tuning considerably and improves the operating stability of the 
Stages. Gates with the following parameters were used in the pump= 
ing circuit: VSWR @® 1.12, isolation 30 - 36 dB. The phase shifter 
¢ is installed in the PA, circuit for selection of the optimal pump- 
ing phase. The VSWR of the parametric amplifiers in the pumping 
waveguide at the pumping frequency was 5 = 7. The electromechanical 
waveguide shorts WS] and WSo are used in the PA pumping circuits for 
convenience in stage-by-stage tuning of the PA, and for rapid check 
of operation of the radiometer stages. 


Characteristics of degenerate PA cascading. In the described 
radiometer, both PA stages are supplied from a common pumping gener- 
ator. As a result of this, the tuning and operation of the 
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parametric amplifier block have certain peculiarities: 1) the 
pumping phase of one of the cascades PAY must be optimized by means 


of the phase shifter $ in order that the overall gain of the two 
cascades be maximal throughout the entire amplification band for a 
given amplification level in each PA cascade; 2) the signal, after 
passing through the first PAY 4s modulated by the pumping frequency. 


The noise signal begins to carry "information" on the pumping phase 
and is processed by the second stage PA, as a synchronous signal. 


As a result, with optimal difference of the pumping phases on the 

two PA, the overall gain of the two stages is 3 dB higher than in the 
case when each of the amplifiers operates with its own pumping gen- 
erator. This makes it possible to reduce the gain of the independ- 
ently tuned PA to 11 dB and thereby improve somewhat the stability 

of operation of the radiometer as a whole. In this case, the overall 
gain of the two PA remains 28 dB. 


Theoretical examination of the operation of degenerate para-= 
Metric amplifiers in radiometers [2, 9, 10] shows that the interac- 
tion of the two bands, primary and unloaded, leads to reduction of 
the effective received frequency bandwidth by a factor of two. The 
expression (4) for the sensitivity of the compensation radiometer 
with degenerate input amplifier has the form 


ar Vie 


for sufficiently high regeneration. Here 2 Af is the actual PA pass- 
band, observed with the aid of a sweep generator, i.e., in the ex= 
pression for determining the radiometer sensitivity, there appears 
the factor 2, indicating deterioration of the sensitivity. This 
effect was confirmed experimentally by two methods. The first 

method involved connecting at the output of the parametric amplifier 
block a band filter with passband equal to one PA sideband — sf 
{either primary or unloaded). Then the PA block passband was reduced 
by a factor of two in comparison with the sensitivity. The losses 
introduced by the filter were compensated in the signal circuit by a 
calibrating attenuator connected ahead of the detector. The 
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measurements were made at a single detector working point (to avoid 
change of the transmission coefficient of the latter) and showed 
equality of the receiver sensitivity values with and without the 
filter. 


The second method involved receiving the signal at the degener- 
ate amplifier block output by a narrowband superheterodyne tunable 
receiver whose passband was on the order of 10 MHz, and was swept 
over the entire parametric ampiifier band. One sideband of the 
BSuperheterodyne receiver was suppzessed by a special preselector. 
The measurement results were analogous to those of the first method. 
In both methods, the noise signal from the detector output was re- 
corded on magnetic tape and then processed by an electronic computer. 
The expected qualitative difference was obtained in the probability 
density distribution functions of the noise signals processed by the 
degenerate and conventional parametric amplifiers [10, 12]. 


Tunnel amplifier block. The TA used in this radiometer has 
the following parameters: gain per cascade about 11 dB; passband 
2 Af ~ 1 GHz; noise temperature T, = 700 - 800° K. In the case of 


cascade connection of the two TA, the resulting Th ~ 800° K, which 


introduces at the radiometer input (with 28 dB gain of the two- 
cascade PA) less than 1° K excess noise. Circulators with waveguide 
cross section 5 * 23 mn were used as the TA block circulators. Elec= 
tromechanical waveguide shorts (WS, and WS,) are used for receiver 


tuning convenience, and when necessary for increasing the radiometer 
dynamic range, which is achieved by activation of the last stage TA. 


Detector. We noted above that when designing highly sensitive 
radiometer receiver circuits, it is necessary to observe definite 
quantitative relationships between the detected microwave noise fluc= 
tuations and the low-frequency video detector noise fluctuations. 

It follows from (3) that for this purpose, a quite high microwave 
circuit gain (50 4B) 18 required in the direct-amplification radio- 
meters. This leads to several undesirable effects. 
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Thus, because of the high gain in the UHF circuit in a single 
frequency band, additional screening and isolation are necessary in 
the microwave circuit, it becomes more difficult to broaden the UHF 
stage passband, the direct amplification receiver dynamic range de-= 
creases, and the stability of the radiometer as a whole decreases, 
which leads to deterioration of the receiver technical sensitivity. 
Increase of the video detector quality factor is a realistic way to 
decrease the gain in high-sensitivity radiometers. Therefore, in 
developing the radiometer primary attent’on was devoted to the problem 
of increasing the video detector sensitivity or, what is the same 
thing, increasing its quality factor. Among the diodes manufactured 
by Soviet industry, the highest quality factor value in the required 
band is provided by the D609 diode, which was used in constructing 
the radiometer detector section. The parameters of the detector 
section are shown in the table. 


TABLE*® 












Prin? W 
for Af = 1 Hz 


VSWR 





cor 
a 
25 


,635 2,85 
8250 0,84 2,25 
8500 0,61 2,04 
8730 0,64 2,69 
9000 0,72 4,8 
0,88 - 
4,45 - 
4,3 - 





# 
Commas represent decimal points. 


We find the quality factor and sensitivity from the formulas: 


an, VOTE? 
a 5 peo a (6) 


where Ry 4s the video resistance at the operating point; 6 is the 


anode current sensitivity; Ry, 4s the video amplifier equivalent 
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noise resistance; P is the detector sensitivity; k iu the Boltz= 


min 
mann constant; To — 293° K; Af is the video amplifier passband 


(equal to 200 Hg, when measuring the detector parameters). 


Special selection of the D609 diodes is required in order to 
realize this high video detector sensitivity. Figure 5 shows typical 
curves of detector section quality 
factor versus frequency in the 8000 - 
9000 MHz band for various D609 diode 
specimens. The high detector sensitivity 
required additional precautionary meas- 
ures to reduce induction effects in the 
detector section and low frequency pre- 
amplifier circuits. Thus, in order to 
break the chassis currents, the detector Figure 5. Detector 
section was galvanically isolated from section characteristics 
the radiometer waveguide circuit and the 
coaxial cable between the detector and the LFPA was constructed with 
a high degree of shielding (double shield). 





Low frequency preamplifier (LFPA). The LFPA was located in the 
immediate vicinity of the detector section, and was constructed using 
two 6S51=N tube stages. The high input resistance and low equivalent 
noise resistance of the video amplifier provided maximal sensitivity 
of the detector=LFPA system. A filter which cuts off the power sup= 
ply induction at 50 and 100 Hz was included in the amplifier circuit. 
The LFPA was placed in a double steel shield, and was damped meche 
anically with paralon (the latter made it possible to reduce by 
approximately an order of magnitude the microphonic effect of the 
mounting). The radiometer low-frequency block was the ShL-2 low 
frequency unit deveioped at the Special Design Bureau of the Insti- 
tute of Radio Engineering and Electronics of the Academy of Sciences 
USSR with output to an EPP=09 recorder. One of the low frequency 
output channels was connected to an analog-digital converter, from 
the output of which information in binary code was fed to the mag- 
netic recording system. 
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Radiometer microwave block. The entire microwave circuit of 
the radiometer was constructed in a single sealed thermostat, which 
was installed at the front foous of the RT=-22 mirror. This same 
thermostat included the pumping generator (PG), power amplification 
block (PAB), power supply ferrite modulator, and the LFPA. Control 
of the microwave block components was accomplished by remote drives 
from the radiometer control console mounted in the RT-22 control 
room. The temperature inside the thermostat was maintained constant 
at +20° + 0.2°, regardless of the ambient temperature, with the aid 
of a semiconductor thermobattery developed at the Semiconductor 
Institute of the Academy of Sciences USSR. A centrifugal fan was 
mounted in the thermostat to reduce the temperature gradient. Meas- 
urements showed that the overall gain drift of the entire microwave 
circuit of the radiometer together with the LFPA is 4G = 0.5 dB ° 


deg, The gain increases with increase of the temperature inside 


the thermostat. 





In view of the fact that the microwave block was mounted on 
the moving portion of the RT=22 mirror, in constructing the radio- 
meter it was necessary to increase the stiffness of the microwave 
circuit structure. Sealing of the microwave thermostat made it pos- 
sible to eliminate the influence of marine climate on the radiometer 
components. The thermostat design made provision for the possibility 
of operation with a slight differential pressure of the nitrogen used 
to f111 the chamber (0.01 = 0.05 at). 


Operational characteristics of the radiometer. The equipment 
Operated continuously without shutdown for several month-long operat= 
ing cycles. During this time, the equipment operated very stably 
while retaining its basic technical parameters. The sensitivity of 
the radiometer on the antenna remained practically unchanged, AT % 
0.02° K with t # i sea (confirmed by repeated recordings of cali- 
brated sources). The instability of the gain of the entire radio- 
meter channel was recorded in the laboratory, and did not exceed 
1.3% during 30 minutes [a ("step") calibrated signal from a gas- 
discharge noise generator (T, = 25.5°) was applied to the radiometer 
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input to check the gain stability]. The stability of the output 
indicator null decreased somewhat during operation of the equipment 


on the antenna because of the influence 
of atmosphere background fluctuations, 
instability of the power supply, and 
periodic radio noise at the RT=22 mirror 
location. On the average, the insta- 
bility of the radiometer null caused by 
atmosphere background variations during 
good weather does not exceed 0.1° K 
(magnitude of the noise trace for the 
sensitivity of the given radiometer). 

In the case of moderate cumulus cloud 
cover, this instability lies in the 
range 0.1 = 0.25° K. During light rain, 
the instability reaches 1° K, which, for 
the RT=22 mirror corresponds to 1° % 18 
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flux units, 1.2., 18 + 107°° w/m* Hz. 


This large influence of the atmosphere Figure 6. Calibra- 


tion signal 
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on padiometer operation (even when using lobe modulation) is ex- 
plained by the small magnitude of the RT=22 antenna aperture. 


Figure 6 shows the laboratory recording of a 3.2° step with hot 
radiometer input 300° K, and for t = 1 sec. The source shown in 
Figure 7c was recorded under unfavorable atmospheric conditions. 


In conclusion, the authors wish to thank Academician A. 5. 
Severnom and I. @. Moiseyev for the opportunity to use the RT-22 
antenna for the observations, and also the personnel of the Crimean 
Astronomical Observatory of the Academy of Sciences USSR, the Insti- 
tute of Space Studies of the Academy of Sciences USSR, the Sate 
Astronomical Institute im. Shternoderg, and the Moscow State Pegago- 
gical Institute im. V. I. Lenin for their active participation in 
the work using the telescope. 
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IF AMPLIFIER LIMITING FREQUENCY SELECTION IN 
SUPERHETERODYNE MM- AND CM=-BAND RADIOMETER 


Yu. A. Nemlikher, I. A. Strukov and 
L. H. Yudina 


ABSTRACT. A technique is proposed for cal- 
culating the limiting frequencies of an inter- 
mediate-frequency amplifier from the experimen- 
tally obtained frequency dependences of the 
klystron generator noise radiation spectral den- 
sity and IF amplifier temperature. An experi- 
mental relation and technique for measuring the 
amplitude-modulated noise spectral density of 
some klystron types are presented. 


The thermal radar method using equipment installed aboard air- 
planes and artificial Earth satellites (1, 2} is being used more 
and more extensively for studying the natural resources of the Earth. 
Such equipment includes the radiometer with high fluctuation sensi- 
tivity AT, whose magnitude can be calculated using the known formula: 


AT = — ree (1) 


where D is a coefficient which depends on the radiometer construc- 
tional characteristics; T is the receiver noise temperature, includ- 
ing the antenna noise; t is the integrator tine constant; Af is the 
receiver passband up to the first detector. We find from (1. that 
to realize high radiometer sensitivity, it is necessary to reduce 
the receiving system input noise temperature T and broaden its pass= 
band up to the first detector. 
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In prineiple, two schemes for the construction of a radiometer 
for the millimeter and centimeter wavelength bands satisfying these 
requirements are possible: a) the direct amplification radiometer 
scheme; b) the superheterodyne radiometer scheme with shift of the 
signal spectrum into the region of lower (intermediate) frequencies. 
For all its apparent simplicity, the realization of the first radio- 
meter scheme in the subject band encounters difficulties with operat- 
ing frequency increase. With the appearance of boardband low-noise 
transistor amplifiers and resistive transformers using (GaAs) 
Schottky varrier diodes, it is becoming possible to obtain simple and 
highly sensitive input reception devices, which makes radiometers 
of the second type preferable [3]. 


The fluctuation sensitivity of the superheterodyne radiometer 
will be maximal only with proper choice of the IF amplifier working 
frequencies. The lower amplifier frequency boundary is determined by 
heterodyne noise and the upper limiting frequency is determined by 
IF amplifier (IFA) noise. In the general case, the noise temperature 
of the superheterodyne receiver, whose block diagram is shown in 
Figure 1, can be written in the form: 


fo boty [(-A)" 4 ++ Fe —], (2) 


where A is the coefficient characterizing the overall noise of the 
heterodyne klystron in the given regime, which is numerically equal 
to the intermediate frequency value at which the increase of the 
receiver optimal noise temperature, owing to heterodyne noise, is 
equal to unity; m is any positive number m > 0; f is the instan- 


taneous frequency, lying in the low-noise IFA passband; te. is the 
antenna relative noise temperature t, * T,/To» where Tp =u 290° K; 


Lo are the mixer power conversion losses; ty is the relative noise 


temperature cf the mixer itself; Fy is the IFA noise coefficient. 
In order to find the average value of the receiver noise tem- 


perature in the IFA band Af (Figure 2), it is necessary to know the 
form of the approximating curve for F, - 1, under the condition that 
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Figure 1. Block diagram of super= 
heterodyne receiver input section: 





1 — antenna; 2 — resistive fre- 
quency converter; 3 — wideband IF 
amplifier; rae a hetero-~ Figure 2, Illustration for 
y finding average value of equi- 
valent relative noise tem- 
the mixer conversion losses Lo perature 


are constant in the subject fre-= 
quency band Af. Experiment shows that the approximating function 
can be a power-law function: 


rt = (f) +4 3) 


where the coefficient B and exponent k characterize the noise proper- 
ties of the transistor IFA, beginning at some frequency. Provided 
that f << B, the amplifier noise coefficient is constant and equal 
to F. el+t ty. This in correct, since the IFA lower limiting fre- 


quency is bounded by the heterodyne noise, and is selected to be no 
less than 30 - 50 MHz. The coefficient B and exponent k are deter- 
mined not only by the transistors used in the IFA, but also by the 
number of amplifier stages, the amplifier circuit solutions, and so 
on. For example, for one of the amplifiers using Q@T329 transistors, 
B = 570 MHz; t, = 0.5; ks 6. The difference between the true he-= 
havior of the amplifier noise coefficient and that calculated using 
the empirical formula (3) was no more than 10%. 


For the case when m ® 1, by substituting (3) into (2), and 
integrating in the frequency band Af, we can find the average re= 
ceiver noise temperature values: 


a Bt (fet fiemy 
yy te Fagin + ae) + Tre ” 
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0 
where Te mT, Lo (tet). Similarly, when ae » the following formula 


holds for the average noise temperature: 


a™ qyic™ — im ) ap Fa — firk) 
Tay" Lit. T= mh +The | + Tre (5) 


Let us find the limiting IFA frequencies for which the mean 
square of the temperature fluctuations has the minimal value. To 
this end, we substitute into (1), Toy from (4) or (5), respectively, 


and find the partial derivatives with respect to fy and f,- After 


differentiating and some transformations, we have 


mn wee 
4 fa 
_ tbat (5 ar ad acer re 
(fa = ft) ( + hKeUds h—oh i (6) 
fa \t “fa 
4.5/1 4 jitk 4 4,5 lu 
= i is pe = SEs) a, 2 
mals m>0: 2 oe 


itk 
emt (EP) aget E& 1S — 1) -s|~0 
(t+ k) (is — fi) K Lye C= /\ i F esha 


ith 
ety Me) Wd td PY (7) 
im Tamm ~ 5} =o 


(t+ eb fe — fi) ik 
Here, the notations are: 


Thus, after determining experimentally the heterodyne klystron noise 
parameters and measuring the IFA noise temperature in the passband, 
we find the values of m, k, Fo: Then, solving the system of equa- 


tions (6) or (7) for the given receiver noise temperature Tre» we 


determine the limiting IFA frequencies for which the radiometer will 
have the maximal fluctuation sensitivity. 


If we use in the radiometer a mixer of balanced construction 
and it is known that the heterodyne noise is suppressed by a factor 


of 6, then in (4) and (5) we must substitute A : g~1/m in place of 
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the coefficient A. The calculation of the limiting frequencies was 
made on a computer. It was found that at the limiting frequencies, 
the noise temperature increase amounts to 3 dB in comparison with 
its value at the center of the passband. 


Let us return to (2) and (3). After substituting (3) into 
(2), we have 


T= LePo| +)" +(-4)"] + Tree (8) 


This formula (8) describes the superheterodyne receiver noise temperature 
and includes the noise of the components making up the receiver. 
However, this same noise temperature can be ascribed to only a single 
component — the IFA, and all the other components (klystron, mixer, 
antenna) can be considered noiseless. Then, in order to find the 
limiting IFA frequencies for which the radiometer will have the mini- 
mal value of the fluctuation sensitivity At, it is sufficient to find 
the frequencies corresponding to the points of intersection of the 
curve plotted using (8) and the horizontal straight line correspond- 
ing to increase of the noise temperature by 3 dB relative to its 
minimal value in the amplifier passband. The required values of 


A, B, k, m, Tey are obtained from exverimental study of the klystron, 


mixer, and IFA used in the radiometer. Thus, for one of the K-45 
klystrons and an IFA constructed using @1'Z29 transistors, the fol- 
lowing values were obtained: a = 230 MHz; B = 570 MHz; m= 2; k ® 6. 
It is easy to find that the lower amplifier limiting frequency is 

fy = 250 MHz, and the upper is fs = 525 MHz for The #0. If Te 7 


1450° K, we find the corresponding values f= 94 MHz, ff," 776 MHz. 


We see from this example that the IFA band is determined by the 
relationship between the receiver self-noise The and the noise 85 
introduced by the heterodyne and IFA. 


For minimization of the superheterodyne radiometer fluctuation 
sensitivity, it is very important to carry out studies of the noise 
properties of the microwave power generators, which operate as 
heterodynes, parametric amplifier pumping sources, and so on. The 
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actual noisy klystron can be represented in the form of two genera- 
tors: an ideal non=noisy generator at the basic signal frequency, 

and a noise generator. If, in the superheterodyne radiometer, we 

use @ mixer operating in the wideband regime, i.e., one which along 
with the signal frequency also receives the mirror frequency, and 

as the heterodyne we use the subject klystron at the basic frequency, 
then the noise generator can be broken down into two noise generators 
with center frequencies W, ~ Wasp and w, + Warps TeSPE. ively, where 


Wh is the heterodyne frequency and Wap is the intermediate frequency. 


Such a regime is widely used for shifting the spectrum into the videc 
frequency band, where its further processing takes place. 


In the general case, the noise generated by the klystron 
generator can be broken down into amplitude-modulated noise (AM 
noise) plus frequency-modulated noise (FM noise) plus the background 
noise [4, 5]. It was shown in [6] that the spectral density of FM 
noise and AM noise as a function of the intermediate frequency magni- 
tude (IF) is the same, beginning with frequencies higher than 20 MHz. 
Consequently, for IF above 30 MHz, the determination of the magnitude 
and nature of the behavior of the AM noise plus background noise 
makes it possible to evaluate the overall klystron noise radiation. 


Figure 3 shows a block diagram of a measurement setup which 
mak?s it possible to measure the total noise temperature of the 
system and the AM noise, plus background noise of the heterodyne. 
We can show that this is so if we consider that the signal spectrum 
for AM noise: 


ot) a Cy[atn ant + AF sin (om + O)E-+ St sin(on— )¢1, (9) 


and for FM noise and small modulation indices m << 1, which 1s the 
case in the klystron, the signal spectrum 


¢ (1) — Colin ot + Fain (om + A)t— Fain (a —0)¢]. (10) 


Thus, the mixer operating in the wideband regime, in which the sube- 
ject klystron is at the same time both the heterodyne and the noise 
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Signal source, transforms at the IF 
Only the sideband AM noise and the 
background noise, while the FM noise 
Side frequencies are not transformed 
at the IF (they cancel one another). 
The line shown in Figure 3 is used 
as a matched detector section and 





Figure 3. Block diagram of 
permits monitoring the degree of measurement setup: 
matching of the mixer with the 1 — noise generator; 2 — 
circuit. measuring line; 3 — direc- 

tional coupler; 4 — step~ 

less HF attenuator; 5 — 

The measured mixer relative gate; 6 — test klystron; 

. 7 — frequency converter; 

noise temperature [7] is: Ae = converted cutout ’ 
», 47 %,@ 4 equivalent; 9 — low-noise 
: =p [A + (1-H) +. (11) IF amplifier; 10 — preci- 
sion IF attenuator; 11 — 

where N. (w) 48 the spectral density indicator 


of the AM noise plus background noise 

of the klystron; L is the noise signal attenuation in the circuit of 
the subject klystron. In mixers constructed using a Schottky barrier 
diode with small values of the loss resistance R,» the quantity to 


4s calculated from the formula: 


4 
t=tgt—z), (12) 
where ta is the relative noise temperature of the Schottky barrier 


diode (for an ideal diode, t, % 0.5). We see from (11) that the 


measured quantity consists of three noise components: AM noise plus 
background noise of the frubject klystron, thermal noise in the sub- 
ject klystron circuit, and the crystal mixer self-noise. By meas- 

uring L, Lo» and tos we can determine the generator noise radiation 


spectral density. 


The table shows « ‘ectral density N, of the AM noise plus backe- 


ground noise for several klystron types relative to the carrier fre- 
quency power level N, referred to the band Af «= 1 kHs. Measurements 
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NL/Ng =102|—105, 


dB/kHz 


VA 





made at an IF f > 500 MHz chowed that the mixer relative noise tem- 
perature does not vary with frequency and is equal to 


4 4 
tag[p+(t-t)|+s=p+% (13) 
We can see from the table that the AM noise plus background noise of 
the klystrons is quite high, and even at frequencies f ~ 200 MH2 N, 
To ¥ 3° 107 °K. When measuring N, (w) of klystrons as a function 


of frequency, it is necessary to know the mixer conversion losses. 
This quantity is found as 


Ly = WT! (14) 


where N, is the noise generator spectral density at the mixer input; 
No is the noise radiation spectral density at the mixer output 
terminals. The quantity No was measured by replacing the HF noise 


generator plis mixer system by an equivalent IF noise generator with 
known and controllable noise temperature. In this case, the IF 
noise generator had impedance equal to the mixer output impedance. 


We investigated the dependence of klystron noise radiation 
spectral density on the magnitude of the voltage Up on the focusing 


electrode. Figure 4 shows the results of measurement of the nojse 
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Figure 4. Experimental depend- Figure 5. Nature of exp- 
ence of type K-49 klystron erimental dependence of 
noise radiation spectral den- redundant relative noise 
sity on focusing electrode temperature of frequency 
voltage: transformer: 


U. = 1800 V; Unet 


1 — IF 30 MHz; 2 — 60 MHz; 
3 — 100 MHz 


= -170 V; 1 — klystron K-45; 2 — 


spectral density of a type K-49 klystron, referred to the mixer oute 
put, for different IF and constant Lo: We see that the noise gen- 


erated by the klystron depends strongly on U, (the voltage Up en the 
resonator and Unef on the reflector were held constant, and during 


the measurement time the mixer crystal current was held constant) at 
frequencies below f < 30 = 50 MHz, and depends weakly on Up at free 


quencies f > 50 MHz. The klystron noise is minimal for Up correspond= 


ing to the maximal generated signal power. Moreover, during the 
measurement, we observed the hysteresis phenomenon, i.e., the minimum 
of the relative noise temperature shifts in one direction or another, 
depending on how Up changes (from -60 to -100 V or from =100 to -60 


Vv). Therefore, an intermediate frequency fon > 50 MHz should be 
selected in mm- and cm-band superheterodyne receivers. 


Analysis of the experimental results obtained (Figure 5) for t' 
made it possible to obtain the empirical formula (2), which was used 
in calculating the limiting IFA passband frequencies. After the 
limiting frequencies are determined, it 1s necessary to construct 
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the amplifier frequency characteristic in some fashion or other 
(for example, with the aid of filters). 
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STUDY OF SCHOTTKY BARRIER DIODE FREQUENCY CONVERTER 
IN THE SHORT MILLIMETER WAVELENGTH BAND 


V. F. Kolomeytsev, Yu. Yu. Kulikov, A. M. Kupriyanov, 
I. A. Strukov, L. I. Fedoseyev, Yu. B. Khapin and 
Vv. S. Yetkin 


ABSTRACT. We present a technique for and 
results of measurement of Schottky barrier diode 
basic parameter’ .t 1.8 - 2 mm wavelengths, and 
also comparison of these data with the results 
obtained on point-contact silicon diodes. We 
discuss the possibilities for the use of Schottky 
barrier diodes in frequency converters for the 
short millimeter wavelength band. According to 
the measurement results, the relative noise tem- 
perature of the Schottky barrier diode is ap- 
proximately 1/2, and the diode ideality parameter 
is 1.2. The frequency converter relative noise 
temperature is close to 1. The measured conver- 
8ion losses in the radiometric regime were less 
than 13 4B at A = 1.88 mn. 


The absence of short millimeter band UHF leads to the necessity  /88 
for constructing radiometers for this band using the scheme: modu- 
lator-frequency converter=IFAedetector. In this case, the radio- 
meter reliability and stability are determined by the reliability 
and stability of the frequency converter, in which point-contact 
diodes are used [1, 2]. The latter have poor stability, reproduci- 
bility, and reliability, which limits application of short millimeter 


142 


band radiometers for the solution of a whole series of problems of 
radiophysics, radioastronomy, and thermal radar. 


Recently, considerable attention has been devoted, both in the 
USSR and abroad, to study ot majority-carrier diodes (Schottky bar- 
rier diodes fabricated using planar-epitaxial technology). Just as 
the point-contact diodes, the Schottky barrier diodes(1) utilize the 
rectifying properties of a metal semiconductor contact. The use of 
planar-epitaxial technology has made it possible to create a "honey- 
comb" structure with reproducible parameters which do not change 
with time. In the Schottky barrier diodes, there is no minority 
carrier accumulation effect and, consequently, there is no diffu- 
sional capacitance. Thanks to this property, these diodes are fast- 
acting varistors, and may find application in microwave switching 
devices and microwave frequency down-converters. The use of a 
Schottky barrier diode in a 5emm frequency converter was reported 
in [3]. Conversion losses of 5 = 6 dB were obtained. As far as we 
know, no experiments have been conducted on use of Schottky barrier 
diodes in frequency converters at shorter wavelengths. 


The purpose of the present article is to investigate the possi- 
bility of using Schottky barrier diodes in frequency converters in 
the short millimeter segment of the spectrum, and to conduct an ex= 
periment using a Schottky barrier diode and a point-contact diode. 

The experimental conditions were the same in both cases. Diodes with- 
out holders were used in both cases. The studies were conducted in 
the same chamber. The dimensions of the pins for the p/p structures 
and the needle were also the same. No effort was made to optimize 

the characteristics of the frequency converter. A gold-GaAs contact 
of n-type was used as the rectifying contact of the Schottky barrier 
diode, a tungsten-silicon contact was used as the point-contact diode. 





(1)nereaftes, we shall use for brevity the term "Schottky barrier 
diode" to mean "Schottky barrier diode fabricated using planar- 
epitaxial technology". 
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The frequency converter transfer, noise, and impedance charac-~ 
teristics are defined by the Volt-Ampere characteristic. Figure 1 
shows the V-A characteristics of two diodes: Si point=contact and 
GaAs Schottky barrier and their equivalent circuit. The V-A char- 
acteristics of both diod- types are described analytically by the 
expression: 


I= 1[exp (fr) —4], (1) 


where n is a dimensionless coefficient; I, is the saturation current; 


q is the electron charge; k is the boltzmann constant; T is the ab- 
solute temperature. We rewrite (1) in the form: 


i= J, [exp (aU) — 4). (2) £89 


The quantity a is determined from the slope of the semilog V-A 
characteristic: 


eel 
a= 2.308257. (3) 
In the ideal diode case, n # 1 and 


a a SE o/h? = 38,00" for 7-2 300° K. 


Figure 2 shows the diode V-A characteristics on a semilog scale, 
For real diodes, the V-A characteristic slope is a = 30 V-1 for the 





%» 
‘& 
Figure 1. V=A characteristics Figure 2. Semilog V-A char- 
of Si pointecontact and GaAs acteristics of Si point-cone 
Schottky barrier diodes and tact and GaAs Schottky bare 
their equivalent circuit rier diodes 
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Schottky barrier dirxzde, and a = 12 yt for the point-contact diode, 
Thus, the parameter n for the corresponding diodes is: n= 1.28 
for the Schottky barrier diode, n= 3.2 for the point-contact diode. 
The saturation currents for the corresponding diodes are determined 
by the points of intersection of the V-A characteristics with the 


current axis. For the GaAs Schottky barrier diode, I, = 4.5 x 1074 


-4 


A; for the point-contact diode, I, 2 10 A. The diode series ree 


sistance R, can also he determined from the V-A characteristics: 


for the Schottky barrier diode, R, = 15 Ohm, for the point-contact 


8 
diode, R, « 20 Ohm. 


An important characteristic of frequency down-converter diode 
quality is the noise ratio Nps i.e., the ratio of the diode nomir.al 


noise power to the noise power of a resistor whose magnitude is 
equal to Rb + R,- It is not difficult to show that for the Schottky 


barrier diode, the resulting relation for Np has the form: 


al2 & 
Nn = (4) 


For small currents through the diode, when R, >> R,» we can obtain 
Np #®0.5ifn=*i1. With increase of the current through the diode, 


Np Slowly increases because of R, noise. 


In order to exclude the influence of heterodyne noise on the 
frequency converter noise characteristics, the intermediate frequency 
is selected in the 3=cm band. The block diagram of the setup for 
measuring the noise ratio and noise temperature of diodes in this /91 
band is shown in Figure 3. The measurement error amounted to 1° K. 
The noise temperature measurement results of the corresponding 
diodes are shown in Figure 4. The GaAs Schottky barrier diode has 
significantly lower noise temperature in comparison with the 81 
point-contact diode. Another disadvantage of point-contact dioder 
is the considerable noise temperature change with “overheating" of 
the points (shaded region). Curves of the noise ratio of the 


145 





Figure 3. Experimental setup for 
determining diode noise ratio and 
noise temperature 


1 — waveguide chamber; 2 — ex- 
ternal constant bias; 3 — meas- 
uring line; 4 — modulator; 5 — 
precision attenuator; 6 — noise 
generator; 7 — circulator; 8 — 
random signal generator; 9 — 
reference voltage generator; 
10 — tunnel amplifier; 11 — 
detector; 12 — LFA; 13 — super- 
heterodyne demodulator; 14 —— 
recorder 


corresponding diodes as functions of 
the current through the diodes are 
shown in Figure 5. The theoretical 
and experimental noise ratios of the 
Schottky barrier diode agree. 


Figure 4. Noise temperature 
of Si and GaAs diodes versus 
current through the diodes 





Figure 5. Noise ratio 

of Si and GaAs diodes 

versus current through 
the diodes: 


1— theory; 2 — ex=- 
periment 


The basic parameters characterizing the frequency converter 
are the conversion loss L and relative noise temperature t. The 
technique for measuring the relative noise temperature does not 
differ from the teohnique for measuring diode noise ratio. The 
relative noise temperature was found to be equal to t = 1.2 and 
t = 2 = 2.5 for the Schottky barrier and point-contact diodes. 


The conversion losses 
were measured using the 
scheme shown in Figure 6. 
The input signal was the 
noise signal of a GSh=-8 4mm 
tube. Since the spectral 
density of its radiation 
decreases sharply in the 
2=-mm part of the spectrum, 
the tube was calibrated at 
the wavelength A = 1.88 mm 
with the aid of a detector 
receiver. A blackbody in 
boiling nitrogen was used 
as the primary reference. 
The frequency converter 
conversion losses measured 





Figure 6. Experimental setup for 
determining conversion losses: 


1 — frequency converter; 2 — hete 
erodyne; 3 — GSh-8; 4 — modulator; 
5 — precision attenuator; 6 — 
noise generator; 7 — tunnel ampli- 
fier; — reference voltage gener- 
ator; 9 — detector; 10 — LFA; 
11 — superheterodyne detector; 
12 — recorder 


in this fashion in the two-channel regime at A = 1.88 mm were 13 and 
16 dB for the Schottky barrier and point-contact diodes, respectively. 
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INFLUENCE OF PHASE SHIFTER ON FREQUENCY 
DIVIDER CHARACTERISTICS 


Ya. E. Veyber 


ABSTRACT. We study the influence of a 
phase shifter in the feedback loop on circuit 
characteristics for the example of a frequency 
Givider with converter and amplifier. It is 
shown that use of the phase shifter is not ad- 
visable in most cases. We find the external 
signal amplitude at which the synchronism band 
ota a@ maximum for even division ratios 
n2 4. 


INTRODUCTION 


In the last 10 - 15 years, frequency dividers (FD) of the re- 
generative type with multi-order frequency conversion have been 
developed and studied (1, 2, 4, 7 - 9], which have a simple circuit 
(absence of frequency multiplier), wide synchronism band, and ease 
of regulation. This variety of regenerative FD includes the fre- 
quency dividers with converter and amplifier (DCA) [1 - 5], the 
push-pull FD (7, 8], and the push-pull FD with reactive feedback 
(9, 10]. In these FD there are no asynchronous oscillations outside 
the synchronism band [11] because of the smallness of the so-called 
asynchronous current component (Ip4)- Some researchers (3, 6, 12] 


have suggested the introduction of a phase shifter into the feedback 
loop in order to improve the FD parameters significantly. According 
to these investigators, the introduction of "optimal" phase shift 
maximizes the mixer transmission coefficient, which inoreases the 
synchronism band and stability of the division process. Accordi:g 
to Andreyev and Tseitiin (4, 13), FD circuits without phase shift 
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in the feedback loop (i.e., without a phase shifter) are usually 
used, as @ result of which maximal synchronism band and maximal 
output voltage phase stability are obtained. However, there is a 
phase shifter in the circuit ywoposed in [4]. Reference [14] also 
indicates the absence of any marked improvement in the synchronism 
band in the FD with small value of Toi with the introduction of 
complex feedback. 


In the present article, we investigate the influence of a phase 
shifter on FD characteristics and present recommendations relative 
to its use in FD circuits. We study the influence sf the phase 
shifter for the example of a DCA cirecvit whose transistor versicn 
4s shown in Figure 1. The frequency converter (FC) is a bridge con- 
sisting of the diodes Dd, - 


Dy» to which there is ap- 


plied the input signal e 
and the feedback voltage 
u from the cascade ampli- 
fier consisting of the 
transistors T, and To. 


We make-the study using 
the method of slowing 
varying amplitudes, we use 
the symbolic method of 
{15] for writing the equa- 
tions, and we use the apparatus of double Yourter series in the form 
of the modulation characteristics method [11, 16] for calculating 
the current spectrum at the FC output. 





Figure 1. Frequency divider with cone 
verter and amplifier 


FREQUENCY DIVIDER EQUATIONS 


We shall make several assumptions which simplify analysis of 
the circuit (Figure 1). 
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l. We consider that the amplifier operates in the linear re- 
gime and there is no collector current reaction (as, for example, 
in the cascade amplifier). 


2. The inductance of the feedback loop winding has no influ- 
ence on FD bridge operation. The admissibility of this assumption 
is determined by the small magnitude of the modulus of the feedback 
coefficient k for transistor circuits (k + 0, 1). 


3. The stray inductance of the circuit coil is equal to zero. 
For the frequency range 0.1 - 30 MHz in FD circuits, coils are used 
with magnetoelectrically armored cores, having small stray inductance 
even in comparison with the feedback winding inductance. 


4. We assume that in the amplifier, transistors are used with 
limiting current amplification frequency fy» exceeding by 5 = 10 


times the subharmonic frequency fis so that we can ignore the iner- 


tial properties of the current carriers in the transistor. 


Let us formulate the equations of the DCA. Assume that the ex- 
ternal synchornizing signal e # E cos wt, and the feedback voltage: 


um U (1) cov lat + 9 (4), ; (1) 
act on the FC bridge. The complex feedback voltage amplitude: 
U = Ue (2) 


is connected with the complex amplitude of the current through the 
loop I, and the symbolic control impedance Zy (p) by Ohm's law: 


U =1,2,(p). (3) 


The magnitude of the symbolic impedance for the circuit with a 
Single parallel loop is equal to: 


R R 
210) = Te pER® = TEAS cH) 
where Roe is the equivalent active resistance of the loop with time 
constant T at he natural frequency Wos P ® d/dt ia the differential 
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operator; €& = (w/n - Wo) T 1s the generalized loop detuning; k = 
ok « U/U, is the feedback coefficient with modulus k and $5 
UL is the complex amplitude of the loop voltage. Substituting the 


expression presented below for the complex current amplitude (12) 

and also (2) and (4) into (3), separating the real and imaginary 
parts, we obtain a system of two abbreviated equations describing the 
processes in the FD in the synchronism regime: 


TpU =U (Ny, cos np-cosq, + Ny,:lungp-siu p, — 4), 
Tpp = — (Ny, sin np-cosg, — Ny, cosa np-sin g, + 8), (5) 


where N = SRoek is the regeneration coefficient; s, is the equiva- 


lent slope of the characteristics of the FC-amplifier circuit, a 
t 


are the in-phase and quadrature with respect to feedback voltage 
subharmonic current expansion coefficients (EC). 


RESULTS OF FC OPERATION ANALYSIS 


Using the double Fourier series apparatus [11, 16], we deter- 
Mine the instantaneous value of qT, and the complex amplitude of 


a of the current components (with frequency w,* w/n) falling in 


the LC loop passband for the synchronism regime: 


fam = RySq[ fa, n-1 008 ($< +9 —n9) +h, nes 008 (Sa +O+ ng) |, (6) 
I, = (f,cosnp — j/,slnn)e!, (7) 

where 
lea = RSet + hh, wa) (8) 


Here, 5, is the slope of the transistor (amplifier) characteristic; 
n is the division ratio; R, 4s the equivalent resistance of the FC 


load, determined by the amplifier input resistance Ran and the con- 
verter load resistance r, (see Figure 1). In components of the form 


Thy the first subscript denotes the number of the current harmonic 
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through the FC load with ex- 
pansion into a Fourier series 
in the external signal fre- 
quency w, and the second sub=- 
seript is the harmonic number 
with expansion of I, into a 


Fourier series in the cyn- 


chronous frequency w, 3 w/n. 


We shall approximate the Volt- 
Ampere characteristic of the 
FC bridge diodes by a polygon 





Figure 2. Approximation of 
with cutoff (E, > 0 is the modulation characteristic by 
cutoff voltage, ry is the re- polygon: 

The position of the operating 
sistance of the open diode) point on the characteristic is 

determined by the signal ampli- 
in order to take into account tude E; for large feedback 
the influence of the cutoff ae amplitudes a two cute 

Z off points appear; — theo- 
voltage E, on the FD para retical data; 2 — experimental 
meters, since for the small data; 3 — approximation 


FC loads characteristic for 
transistor circuits, the voltages applied to the FC are commensurate 
with the cutoff voltage E,° In this case, the amplitude of the 


first current harmonic I, (u) with frequency w (modulation char- 
acteristic) is equal to 


hud { (9), 1u— 28, <0, 
aya 0,5, u—2E.>0, 
where y, (0) 4a the EC for the first current harmonic [16], and the 
diode bridge dynamic slope Sayn and current cutoff angle 9 are de- 


fined by the equalities 


4 — 28, 
Sayn ° WER PAF C080 = = a. (9) 


In order to facilitate the calculation of the components 143, 
we approximate the modulation characteristic I, (u) by a polygon, 
considering, in so doing, the experimental data (Figure 2). As a 
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result of the calculations, we obtain: 
14,” Sayn Esa [y, (8,) - 7 (6, = 0.417 Sayn UU: 


0.417 Sav, U Cy, (82) - vz (6,0), (10) 


dyn 


where yz (8) 18 the EC for the 2* current harmonic [16]; 5 = 
0.5/1.2 = 0.417 is the slope of the approximating straight line; 
a = U/E and b @ 2 E,/E are the relative values of the feedback and 


bias voltage amplitude; 8, and 8. are the lower and upper current 


cutoff angles, defined by the equalities: 





ens, ===, cus, == oe = cos8, +42. (11) 


With acccunt for (10), the expressions (7) and (8) take the form 





i oe sa re a np — /7, slung) e*, (12) 
where 
Tea = [Fn (4) — Teil yl te (Tne (95) — "fner (Q)l, (13) 
Sq= O11 Saynt 
.The graphs of the EC 2 7 
y (a) for various values 
of the bias b and divi- 44 
sion ration # 2, 4, 5, | il 
constructed using (13) 4 ; Bull 


and (11), are shown in 
Figure 3. We see that 
the bias b (or ampli- 





PAN 

PTAA il 

tt Ni 
PT ON 






tude of the external = S 

signal E) has marked ’ Sra wea 

influence on the nature a ae RE 8 aR A eee 
of the EC curves, which 

confirms the necessity Figure 3. Expansion coefficients y 
for considering the biases b for aivioion ratios, "7" 
cutoff E,. a—ne2b—nehks om—ne 5 


(Figure continued on following page) 
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Let us examine the influence of phase shifts in the feedback 
loop on FC operation. According to (6), the amplitude of the sub- 
harmonic current in the FC load depends on the feedback voltage $, 
and the subharmonic phase w *® does not coincide in the general case 
with the phase ¢. It may be shown that, in the region of stable values 
of the amplitudes a, the signs of the components Ty n-2 and Th ntl 


are opposite. Therefore, according to (6), the maximal and minimal 
subharmonic current amplitudes will be obtained for phases equal to 


rmx Ze Pr mex =F 5 (8 — 1), - 92 aio = Yr min = O. (14) 
The phase angle a of the first loop current harmonic (relative to u), /97 
corresponding to the current extrema, can be found on the basis of 


the equation of FD phase balance, formulated under the assumption of 
absence of phase shifts in the amplifier (Figure 4): 


where % 2 %% — is the phase shift introduced by the phase shifter. 
From (14) and (15), we obtain the following values of the angles «a: 


ay max = 0, w = + 7/2; 
Gy a oe 0, tH = 0, 





Figure 4. Block diagram of DCA for phase 
deviations: 


FC — frequency converter; FS — phase 
shifter 


Consequently, the largest amplitude of the current components with 
frequency w/n at the FC output is obtained with feedback voltage 
phase » = + t/2n, or with load (loop) phase angle a = + 7/2. An 


angle close to this value can be achieved only at the edge of the 


SR NI 
a The phase is indicated relative to the phase of the signal e. 
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synchronism band. When the loop is tuned to resonance with the 
subharmonle frequency (with the feedback loop open), the current 
amplitude will be minimal. To obtain a maximum subharmonic current 
with resonant tuning of the oscillatory system, we can use a coupled 
system of loops, or a phase shifter which introduces the additional 
phase shift ¢ «= + 1/2. We emphasize that the indicated phase shift 
magnitude is optimal only in the sense of obtaining maximum sub- 
harmonic current at the FC output, but the characteristics of the FD 
as a whole may not be improved in this case. 


STATIONARY REGIME CHARACTERISTICS 


The FD behavior in the stationary regime is described by the 
equations: 


Te" Wawa" § = —N7, sinnp, @ = 0, (16) 
tewaere Bett cone, P= +a/2, (17) 


obtained from (5) by the substitutions pU = pd = 0 and o, ©®$s0, 


0/20), We obtain all the basic stationary regime characteristics 
from (16) and (17). 


It 1s well known that for n > 2 in the FD without asynchronous 
component Toi? subharmonics are not excited from the equilibrium 


state [11]; therefore, in place of the threshold with respect to the 
input signal E, 4t is more convenient to use the concept of thresh- 
old regeneration coefficient N,, (suitable for any n), to which there 





corresponds the minimal controlling resistance magnitude, when the 
conditions of exigetence of the stationary regime in the absence of 
detuning, € = 0, are still satisfied: 


4/feme, D = 0, 


Neh * | Arenas, Data, 


(18) 





(ane resonance and phase characteristics are symmetric only for 
the phase shift @ 9, n/a. 
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Here Ye,s max denotes the maximum absolute values of the functions 
Yo and Yg° It is obvious that synchronism is possibile only for 
N > Nen' Figure 5 shows the dependence of Nen on external signal 


amplitude E, plotted from (18) and the EC curves; the dash-dot line 
corresponds to infinite threshold growth for b 2 1, when the FC does 


not operate. The curves show rapid increase of Nen with increase of n, 


the presence of a threshold with respect to external input, and the 
absence of a ceiling for the DCA cir- 

cuit. For n = 2, the presence or N 
absence of the phase shifter in the 

circuit does not alter the FD 

threshold properties, for n > 2 
introduction of the "optimal" phase 

shifter into the circuit reduces the 
threshold magnitude in the best case 

by a factor of 1.5 = 2. 


th 





We note the presence for n > 2 


of the function New (V) for even Figure 5. Threshold re-= 
generation coefficient Nin 
division ratios n and sharp increase versus external signal am- 


of the threshold for odd n in the plitude E 
vicinity of the point b = 0.4 (or 
E# 5 Ey)» corresponding to the midpoint of the modulation character= 


istic. This means that the synchronizing signal amplitude E = SE, 


is optimal for even division ratios n 2 4, since in this case the 
synchronism band reaches a minimum; for odd ratios n 2 3 in the 
region of E values close to 5 E,? subharmonic generation will 
terminate. 


We obtain the resonance a (€) and phase n¢ (§) frequency char- 
acteristica(1) from solution of the stationary regdiine equations (16) 








(1)we shall consider that the amplitudes and phases of the voltages 
Usup and u (see Figure 1) coincide, although in reality, for @» 1/2, 


the phases of these voltages differ by 1/2. 
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and (17). In order to calculate the characteristics for ® = 0 using 
the selected regeneration coefficient N > Nene specifying values of 


cos nd and determining sin ng, we find from the first equation (16) 
the quantity y,. From the Yo. (a) graphs, we find the self-main- 
4 


taine . oscillation emplitude = and the EC Ya? and from the second 


equation (16) we find the value of tne detuning &. We solve (17) 
Similarly. The constructed frequency characteristics for the bias 


b = 0.2 and n = 2, 4, 5 are shown in Figure 662), For the division 
ratio n # 2, the frequency characteristics are reminiscent of the 
characteristics of an oscillatory loop. With increase of n and N 

for ¢ # 0, the DCA frequency characteristics improve (1.e., the 
resonance curves become flatter and flatter, and the phase curves 
become prattically linear), while for @ = 1/2 the characteristics 
deteriorate; the resonance characteristics become sharper and sharper, 
the length of the phase curve linear segment decreases, and the slope 
of the characteristic curve and the phase rate of advance increase. 
With n > 4 for the case @ = 0, and with n 2 5 for the case @ = 1/2, 
the frequency characteristics become multivalued because of the 
oscillatory nature of the EC; excitation of large-amplitude oscllla- 
tions becomes most likely, although in practice we can observe ex- 
citation of small-amplitude oscillations. Transition from one branch 
of the characteristic to the other may be accomplished by jump of 

the phase nd by +7 (Figure 6c), which is undesirable. 


Analysis of the frequency characteristics showed significant 
deterioration with presence of the "optimal" phase shifter in the 
circuit. Aotually, in view of the peaked nature of the resonance 
curves (for ¢ = 7/2), with even small detuning of the loop relative 
to the input signal frequency (for example, aging of the oscillatory 
loop components, operation over a range of temperatures, signal fre- 
quency drift, and so on), marked reduction of the output voltage 
amplitude 1s possible, which may lead to disruption of the operation 





ane characteristics of ou QCA with phase shifter are denoted 
by the subscript "@ = = 90°! 
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r i) =§o* 





Figure 6. Frequency cha.acteristicsa of DCA for bias b @ 0.21 
amns#e2 
(Figure continued on following page) 
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Figure 6. (continued): 
benef 
(Figure concluded on following page) 





Figure 6. (concluded): 
cmne § 
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of the subsequent devices and practical reduction of the operating 
frequency band, although the synchronism band does not change. On 
the other hand, for ¢ =» n/2 for any n, the phase characteristics 
remain nonlinear and the phase advance rate increases significantly 
(for example, for n= 5, N # 40, b = 0.2, the phase advance rate 
increases by 18 times), which leads to a corresponding increase of 
the phase instability. In addition, it is not possible to regulate 
the degree of resonance curve narrowing. 


The synchronism band (4.e., the region of maximal possible de- 
tunings within the limits of which the synchronism regime is main- 
tained, if it was previously established) can be determined on the 
basis of the stationary regime stability conditions, or directly from 
the resonance characteristic, since stability transition takes place 
at the boundaries of the synchronism band. It is convenient to ex- 
press the synchronism band I, in units of generalized detuning g. 


by the known equality: 


a 
& = 0, 
where fy and Q are the natural frequency and equivalent quality 


factor of the FD loop. Figure 7 3 shows curves of DCA synchronism 
band, expressed in units of detuning bo versus regeneration coeffi- 


cient N for division ratios n # 2, 4, 5. The relation bo (N) for 
the case n * 2 18 given by the formula: 


RoaYV—t; Omb<t; O=d,n/2, (19) 


obtained upon substituting into (1€) or (17) the limiting EC values 
Y. "Ys "1+ The graphs show clearly that upon connection of the 


“optimal” phaee shifter into the FD circuit for the range of regene- 
eration values N > (1.5 = 2) Nen used in practice, the synchronism 


band usually decreases, since, in the circuit with the phase shifter, 





. Figures 7a, b are presented for the primary branch of the reson- 
ance characteristic. 
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Figure 7. Dependence of DCA synchronism band, 
expressed in units of detuning bo? on regen= 


eration coefficient N: 


a-—n 2.4; B—ne# 5;1—%20; 2— 
@e« t/2 


the synchronism band is determined not by the quadrature component 
Ys? but rather by the in-phase component Yo which, for the indicated 


values of N, is usually smaller than Ys° 


Thus, in contrast with the assumptions of certain authors [3, 
12], we can consider it to be established that the so-called "opti- 
mal" phase shifter does not increase the synchronism band. The 
nature of synchronism band dependence on the external signal ampli- 
tude (Figure 8) confirms this cor- 
clusion. We see from the curves 
for n 2 4 the synchronism band a 
is maximal for the external sig- : 
nal amplitude E = 5 Eo: 





The amplitude characteris- 
tic, i.e., the dependence of the 
feedback voltage amplitude on Figure 8. Dependence of syn- 
the external signal amplitude Siero eat aapiteddee sig- 
for the resonant loop tuning 1—¢@20; 2—¢6=8 2/2 
case &€ = 0 andn ® 2, 4, 5; 

@ = 0 is shown in Figure 9. 
It is convenient to plot the oharacteristic in the coordinates (a/b, 
1/b), since a/b = (U/E) + (E/2 E,) =» U/2 E, and 1/o = E/2 E,: 
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The values of a and b are determined 
directly from the EC graphs, drawing 
@ horizontal line at the level 1/N 
until intersecting the curve y, (a) 


for ¢ = 0 (or the curve Y, for $ = 


7/2). The amplitude characteristics 
are nearly linear. In the input 
amplitude region 1/b = 2.5, we ob= 
serve disruption of the synchronism 
regime for odd n. Depending on the 
regeneration magnitude, hysteresis 
in the oscillation onset and termina- Figure 9. Amplitude 
tion is possible. The hysteretic part characteristic of DCA 
of the characteristics is located to 

the left of the dash-dot line. The phase shifter does not alter the 
nature of the curves, but for ® = 1/2 they lie considerably higher, 
since in the case the self-oscillation amplitude a is larger. Marked 
increase of the oscillation amplitude may cause deterioration of the 
output voltage form (in the transistorized scheme). 





Thus, summarizing, we note the following: 


1. The subharmonic current amplitude at the FC output is mini- 
mal in the absence of phase shifts in the feedback loop, and is maxi- 
mal when connecting the “optimal" phase shifter into the circuit, 
which introduces the phase shift @ = + 1/2, corresponding to tne / 
feedback voltage phase ¢ = + 1/2 n. 


— 
oS 
nD 


2. The introduction of the phase shifter (¢ = 1/2, n > 2) 
may reduce the regeneration threshold magnitude by a factor of 1.5 
to 2. Here, the synchronism band decreases in most cases, and the 
frequency characteristics deteriorate, which leads to inoreare of 
the phase instability. Deterioration of the output signal shape is 
also possible. These negative phenomena are accentuated with in- 7103 
crease of the division ratio n. 
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3. The use of the phase shifter for frequency division by a 
factor of two is not advisable in general beceuse of the absence of 
any advantages. In those rare cases of FD use when phase stability 
does not play a significant role and, figuratively speaking, the 
simple fact of frequency division is sufficient, for n > 2 a phase 
shifter can be introduced into the circuit in order to lower the 
regeneration threshold. In all other cases of FD use (in devices for 
synchronizing and frequency tracking in phase-sensitive systems) 
circults without. a phase shifter can be recommended in order to 
improve FD phase stability. All that we have said here also applies 
in full measure to frequency dividers with reantive feedback (9, 10], 
so that their application is not advisable. 


4, In this study, we have determined the optimal external uig- 
nal amplitude for which the FD has the maximal synchronism band and 
minimal threshold for n 2 4. We have identified the existence of 
a region of "critical" input amplitudes where generation of odd subhar- 
monics of multiplicity n 2 3 1s not possible. The existence of 
hysteresis in the onset and termination of the oscillations has been 


established. 


The author wishes to thank Prof. 3. I. Evtyanov and V. V. 
Andreyanov, Cand. Tech. Sci., for their review of the manuscript 
and valuable comments. 
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COMB-LINE BANDPASS FILTERS 


Ye. A. Vlasov 


ABSTRACT. A calculation is made of comb-line 
microwave filters with respect to specified losses 
with the aid of a prototype low-frequency filter. 
We examine a technique for calculating shielded 
parallel coupled lines with circular cross section 
conductors. Along with the formulas and graphs, 
we give an example of the calculation of a comb- 
line bandpass filter with respect to specified 
parameters. We examine specific designs of such 
filters in application to their use in spacecraft 
ondoard antenna and feeder devices. 


The radiotechnical system of the modern spacecraft includes 
several radio receiving and radio transmitting devices which can 
operate simultaneously in various frequency channels, providing 
two-way radio communication, transmission of telemetry information, 
and the conduct of various physical measurements. Because of the 
comparatively small dimensions, the antennas of these devices are 
inadequately isolated from one another, whic. leads to the appear- 
ance of mutual interferenoe. For interference suppression, the 
channels of the devices include bandpass filters, each of which 
passes the signal band of its device and suppresses the noise sig- 
nals. In addition, a single antenna 14s often used for simultaneous 
operation of several devices on spaced frequencies. Bandpass fil- 


ters are usually used for combining or separating the signals of 
these frequencies [1]. 


167 


/104 


Comb-line bandpass filters Tet of Tf oh Te gt 
(BPF) [2, 3] have been widely ¢ 
used in the microwave (SHF) 
band. Figure 1 shows the elec- ayy a] y 
trical schematic of such a fil- 
ter. it contains n parallel 
coupled metal bars 1, ...,n 


of length 7 = A/8 or shorter. Figure 1. Electrical diagram 
The bars are grounded at one of th a filter 


end, and loaded by the capici- 


tors ck at the other end. The end bars 1 and n are connected in 

parallel with the matching transformers 0 and n+1, which are also 
grounded on one side, while the other side is connected to the cir- /105 
cuit. Such filters are small, since the bars are shorter than A/8, 

and simple to fabricate. 


Tuning of the filter 1s accomplished by the capacitors, which 
makes it possible to compensate for fabrication errors. The filters 
have low losses, since the use of dielectric materials is minimized. 
The first parasitic passband is located no nearer than the fourth 
harmonic. 


Filter frequency characteristic. The frequency characteristic 
of the ideal bandpass filter is shown in Figure 2. The filter 
passes a signal without loss in the frequency band Af = fo - fp? 
and has attenuation equal to infinity outside this band. However, 
such a characteristic 1s physically unrealizable (4). The real 


filter has some attenuation La in the passband; this attenuation 


depends on the frequency. Usually it 1s approximated in the pass- 
band in the form of the maximally-flat (Figure 3a) or Ch3byshev 
(Figure 3b) characteristic. Outside the passband, the filter 
attenuation also depends on che frequency. 


In addition to the basic passband,microwave filters have para- 
sitic bands which must be considered in their design. Usually, the 
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f, # 
Figure 2. Frequency Figure 3. Maximally-flat BPF char-=- 
characteristic of ideal acteristic (a) and Chebyshev 
BPF BPF characteristic (b) 


filter frequency characteristic is specified for its design: pass- 
band center frequency fy) = (f., + f_p)/2s passband relative width 


we Af/fy = (f,, - f_p)/fo> the pulsation amplitude in the passbsnd 
for the filter with Chebyschev characteristic or the attenuation Li 


at edge of the passband for the filter with maximally-flat frequency 
characteristic, the frequency fo. or fs outside the passband at 


which the attenuation La must be provided. 


Reference [3] presents several graphs for determining the number 
of resonators of filters with Chebyshev or maximally-flat character- 
istic. It 1s convenient to determine the number of resonators for 
any filter with Chebyshev characteristic from the nomogram of 
Figure 4 [5]. 


Example. Let us design a combeline BPF with Chebyshev charac- 
teristic having the following parameters: passband center frequency 
fy # 1222 MHz; passband width w # 0.08; pulsation amplitude in the 


passband Ly «= 0.1 dB; the filter must provide attenuation L, = 50 dB 


at the frequency f. = 1616 MHz; the filter 1s connected to a line 


a 
with characteristic impedance zy # 50 ohm. We determine the normed 
frequency Variable: 
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Qs <2. ja f= 5(42-1), 


in our case, Q= 4 os (iaa3 — 1)=8 . 


From the nomogram (Figure 4), we find that the value n = 3 corre- 
sponds to the values Q = 86, L, = 0.1 4B, L, = 50 4B. 
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Figure 4. Nomogram for determining number n of 
resonators for filter with eee eee eharac teristic 
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Resonator parameters. Historically, the development of micro- 
wave filter theory was preceded by the development of low-frequency 
filter theory. A design technique was developed for the latter, 
and extensive tables of values of the individual parameters were 
compiled. The low-frequency filters were taken as prototypes in 


developing the technique for microwave filter design. Specifically, 
the low-frequency filter with lumped elements is usually used as the 


prototype for the design of comb-line filters [2, 3] (Figure 5). 


ri "9s ta'Sa 


oo om om ol 
7 “af be "hy ; bts 1 "Gn faet ™Saet bin 


Figure 5. Schematic of ladder filter used as 
low-frequency prototype 


Tables of prototype filter parameters for Chebyshev and maximally- 
flat characteristics with number of resonators from 2 to 15 are 
presented in (3, 6]. 


The microwave comb-line filter resonators are characterized by 
the characteristic impedance Za and electrical length 9. Figure 6 


shows coaxial line quality factor as a function of characteristic 
impedance [3]. The optimal quality 
1s obtained for Ye,-Z,=740hm. With 


some approximation, this value of 
the characteristic impedance was 
taken for comb filters with bars of 
rectangular section, and yielded 
quite good results [2 = 4]. Good 
results are also obtained for the 
choice 2a = 74 ohm in the case of 





Figure 6. Coaxial line 
quality factor versus 


circular bars. The electrical characteristic impedance 
length of the comb fiiter resonator 
bars is usually 0 =anlacad, 
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Example. For a filter with n = 3 and La = 0.1, we have the 
following values of the elements: @, = 1; f, = 1.0315; g, = 1.1474; 
B3 = 1.0315; Sy = 1. 


From the graph (see Figure 6), we select fay # 74 ohm, which 
corresponds to Yay * 0.0135 mho. It was previously assumed that 


Za = 50 ohm, i.e., Yq = 0.02 mho. 
- Y, 
y, = 0,67 4, vy = 1/8, 


Determination of partial capacitances. Determinetion of the 
partial capacitances is preceded by determination of the parameters 


which are functions of the electrical length 8: 


Go Ce Fijry 20/ . 
Sa Yo "Ya yr, jon +-(1=1) 


It is convenient to determine these 
parameters with the aid of nomograms and 
graphs. Figure 7 shows the relation: 


yeh 1)= otg 0 + Gcosec?O (2) 


Figure 8 shows the nomogram for deter- 
mining 


4y/ Voy 


CG Gh b ) h 
va Py = o(y4)] ” aa v5 (3) f 





and 
Isat b 
“0 = 9-95) = ayer ve 8 cH) 


for the parameters of the prototype filter 
with n ® 3 and L, # 0.1 dB. In the left 
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WV ie Figure 7. °arameter 


b IX, versus resona- 
tor electrical length 
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Side of the nomogram are plotted curves of the relations: 


for w= 0.1. The right side of the nomogram is used to determine 
the values of these parameters for w «0.01 - 0.1. Nomograms for 
any other values of n, La» w can be constructed using this same 
principle. 





a kis it ar io TT a) a ae ayy 3 


Figure 8. Nomogram for determining the parameters 
G./Vo5 (p and (4, 4,,/%o;) g0%) for a filter with Chebyshev 
characteristic n = 3 and Ly = 0.1 dB 


In view of the fact that the values obtained with the aid of 
the nomograms are normed relative to Yay (normed relative to Ya in 


Mattaei [2]), the formulas for determining the capacitances change 
somewhat : 


SoA VR): (5) 


go? vit pt yttes) +e, (6) 
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Le na rare Apt te — “pha, (7) 
Fo etlt- vA . He — “put gt) 4 Se, (8) 
Sp (1 BH) ° 





es (10) 
Ch. fe SN TY) Fy 4, 
sg, = (Sp) an 
Sooner 2 STATY a Suny, 
by Ve 8 . (12) 


With the aid of the nomogram (Figure 9), we can determine the 
load capacitor capacitances: 


“Ix-fo (13) 
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Figure 9. Nomogram for determining capecitances of 
lumped loading capacitors for several resonator 
electrical lengths: 


1 — 6 = 25°; 2 — 35°; 3 — 45°; 4 — 55° 
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The relation of = @ (f) forY,,;=0,01 mho for several discrete values of 


@ is plotted on the left side. The right side of the nomogram serves 


for determining the quantity cy 


0.01 to 0.02 mho. 


for arbitrary Yas in the limits from 


Example. With the aid of the graph (see Figure 7), we determine 
by/Yay? For 0 = a/4 = 45°, b/Y,,= 1.28. For the found value of bsYay> 


using the nomogram of Figure 9, we determine fae, and ; 


J, J 
32 ig0 os 2 tg 0 
a £ v; gY, 


To the quantity b,/Ya3 = 1.28 corresponds Gees = 0.124 for 
w = 0.1. In order to determine Ge/%ay for w = 0.08, we draw from 
the point Gees = 0.124 on the vertical axis, a sloping ray to 


intersect the point 0 of the horizontal w axis. The point of inter- 
section of the ray with the vertical straight line w = 0.08 is then 
projected to the Oe/% a5 axis, and the point of intersection yields 


the desired value Gea/Yoy # 0.099. We can find the quantity Fi ie 0: 
similarly: 


for wa0,t %tg0 = 0,119, 
a? 


for w= 0,08 tg d =: 0,004, 
a 


The bar self=-capacitances per unit length have the values: 


Go og RETRO (1 1 TMVT-OOTT) a 8,58 = St, 
Eo ay SAT POU (1 — 1,48 + 0,089 — 0,004) + 5,58 = 317 a <, 
St, a SEES (1 — 0,004 — 0,004) me 4,18, 
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The mutual capacitances per unit length are: 
Git ag Set 2 7,5 — 5,58 = 1,02, 


Cis Cs, 370,7-0,0195 
Sit Sen, SRT 20,004 = 0,18, 


From the nomogram of Figure 9, we determine the load condenser capa- 
citances. For f) = 1.223 MHz, @ = 45°,with Yay 7 0.02 the capaci- 


tance cK 2.55 pF, while for Yas = 0.0135, the capacitance ck . 
1.75 pF. 


Determination of filter geometric dimensions. The comb-line 
filter is built structurally in the form of a system of parallel 
resonator bars mounted in a single plane, on both sides of which at 
equal distance from the corresponding bars there are located screen- 
ing plates. By filter geometric dimensions, we mean the cross sec- 
tion dimensions of the individual resonator bars, the distance be- 
tween adjacent bars, and the distance between the screening plates. 
Usually, the bars are circular (Figure 10) or vectangular (Figure 11) 
in cross section. The self and mutual capacitances of the resonator 
bars are shown in the figures. 


~ 


Other conditions ChIP LLL LLL ViZiLcaLeaaad 
btn bt det bnet.n 


being the same, the rz ~ A = t 
structures with circu- On, ese rao 
lar bars are somewhat a ee ee | 
simpler to construct. 


We shali examine the Figure 10. Cross section of filter with 
technique for calcu- circular bara 
lating the cruss sece- 

s ED OL Lk de ddide,  abdide, 
tion geometric dimen- fer fg Sond Eager fnehn 


sions of a filter CEH 
with circular bars. a as z = 
The filver structure ath 


4s broken down iliito 


a number of elemen- Figure 11. Cross section cf filter with 
tary cells equal to rectangular bars 


L76 


the number of bars. Each cell is a Teshaped combination of three 
condensers, one of which represents the self-capacitance per unit 
length of the corresponding bar, and the other two represent the 

coupling capacitances of this bar with the neighboring bars. 


o fud “afer fnat 
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Figure 12. Elementary filter 
cells 


Figure 12 shows the filter 
elementary cells. The uni- 
laterally loaded cells 0 and 
n++ correspond to unilateral 
matching transformers at the 
filter input and output. 


The cells 1 - j = n corre- 
spond to the bilaterally loaded 
resonator bars. With the aid 


of Figures 13 - 15, for each 


yeh elementary cell we deter- 


mine the relative bar diameter 
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Figure 13. Relationship between 
2 (Cy 44) and 1/2 (S/b) for 


several fixed values of a/b: 


0.4; 4 — 0.2; 5 —0.1; 6 — 


d,/d, and half the relative distances from the bar to the two adja- 


cent bars: 
4. Abus and tafe 


Figure 13 shows the dependence of the mutual capacitance 
Cc /e between bars on half the distance between bars 1/2 & /» 
JsJt1 J,J+] 


for several values of the relative bar diameter d/b of the jth cell. 
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Figure 14. Relationship between 
1/2 (Cs/e) and 1/2 (S/b) for 


several fixed values of d/b: 


1 — d/b = 0.7; 2 — 0.6; 3 — 
0.5; 4 — 0.4; 5 — 0.3; 6 — 





Figure 15. Relationship be- 
tween 1/2 (C,/e) and d/b for 


several eae ee of 1/2 
s 


Figure 14 shows the de- 
pendence of half the self- 
capacitance 1/2 C;/e of the 


3° bar on half the distance 


between the bars 1/2 8, y+” for several values of d/b. 
s 


In order to determine the dimensions of the bilaterally loaded 


element, we draw on Figure 13, horizontal lines corresponding to the 
element mutual capacitances Cyi3,578 and Cy. j4r/e- The coordinates 


of the points of intersection of these lines with the lines (1) - 
(5) of constant d/b values are transferred from Figure 13 to Figure 
14, and from them we plot the curves for Cyn, y/e and Cy 41/8 


The sought values of 1/2 $421,5/° and 1/2 8, ,g+1/> are the horigontal 
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coordinates of the points of intersection of the constructed curves 
with the curve corresponding to the sought d,/b value. On the other 


hand, the sum of the vertical coordinates of these points is equal 

to the resonator self-capacitance. In order to find the indicated 
points of intersection, we construct an auxiliary curve through the 
midpoint of the chords joining the points of intersection of the 
curves of constant d/b values with the constant mutual capacitance 
curves. On the graph, we draw a horizontal line corresponding to 
half the rod self-capacitance to intersect the auxiliary curve. 
Through the point of intersection we draw a line parallel to the pre- 
viously mentioned chords. The points of intersection of this line 
with the curves Cy.1,578 and Ds, gei7e lie on the sought d4/b curve. 


The projections of these points on the abscissa axis yield the sought 
values of 1/2 S3-3,4/ and 1/2 S,, jai’ 


The sought values of d,/b are found with the aid of Figure 15. 


On the vertical axis, we ley off the partial self-capacitances 1/2 
C,/e of the bars. Along the horizontal axis, we lay off the values 


of d/b. To find the sought d/bd value, we draw in Figure 14, vertical 
lines corresponding to the found values of Cy. jes and Cy, g+17e find their 
> 


points of intersection with the curves of constant d/b values, and 

transfer them to Figure 15. Through these points in Figure 15, we 

draw curves of the found values of 1/2 8,_, ,/b and 1/2 85 443 /> 7113 
3 » 


On the vertical axis of the figure, we lay cff the previously found 
self=-capacitance half-values, which in sum yield the self-capacitance 
Cy /e of the jth bar. Through these points we draw horizontal lines 


to intersect the corresponding 1/2 $5.3,5/° and 1/2 Sy y41/> curves. 


The points of intersection are projected to the horizontal axis of 
the figure to a single point corresponding to the sought dy /b value. 


This technique applies to the bilaterally loaded bars. 
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The unilaterally loaded bars, specifically the matching trans- 
formers at the filter input and output, have some finite mutual 
coupling capacitancs Cy 17€ with the neighboring bar of the filter 


structure on one side, and the coupling capacitance Cy o/£ = 0 on 
3 


the other side. Theoretically, zero coupling capacitance is obtained 
with Sp) o/b = ™ ; in practice, go0d results are obtained with So o/>2 
3 3 


2 1.5.In calculating the cell with unilaterally loaded bar, we take 
the vertical line 1/2 S/b = 0.75, as the line of zero capacitance. 


Example. We determined previously the scolf and mutual capaci- 
tances of the resonator bars. The subject filter can be represented 
in the form of five cells, and, in view of symmetry, cells 1 and 5, 

2 and 4 are pairwise analogous. In Figure 13, we draw two horizontal 
lines corresponding to the previously determined values of Cy 37° oCs iy 


/e #= 1.92 and Cy a/e * Cy 3/e = 0.48, to intersect the lines of constant 
3 


d/b values (curves 1 - 5). We transfer the points of intersection 
to Figure 14, and draw through them curves corresponding to these 
capacitances. The first filter cell is a unilaterally loaded bar 

with the partial self-capacitance Cy/é = 5.58, and the partial 


coupling capacitances Cy o/& = 0, Cy 1/£ = 1.92, In Figure 14, we 
s 9 


make an additional construction: the points of intersection of each 
of the constant d/b lines (1) - (8) with the curve a Cy 1/8 = 1.92 


and the line 1/2 S/b = 0.75 are joined by chords. Through the chord 
midpoint, we draw the auxiliary curve b. Through the point on the 
vertical axis of the figure corresponding to the value 1/2 Cy /es 


# 5.56/2 = 2.79, we draw a horizontal line to intersect the constructed 
auxiliary curve b. Through this point of intersection, we draw a line 
parallel to the previously constructed chords to intersect the lines 

a (Cy | /E = 1.92) and 1/2 S/b = 0.75. The pointe of intersection have 
the coordinates (2. 26; 0.1) and (3.32; 0.75). The sum of the vertical 
coordinates 2.26 + 3.32 = 5.58 « Co/e, which was to be proved, i.e., the 
found points lie on the curve of the sought dy /b value. The 
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horizontal coordinates yield the half-distances: 


To find the dimension d,/bs we find in Figure 14 the points of 


intersection of the vertical lines corresponding to the values 
1/2 So 0/> = 0.75 and 1/2 85 ,/? = 0.1 with the curves of constant 
i] , 


d/b values, and the found points are transferred to Figure 15. In 
this figure, half the bar self-capacitances are plotted along the 
vertical axis, and the a/b values are plotted along the horizontal 
axis. Using the transferred points, we draw curves of constant 

1/2 S/b = 0.1 (curve 2) and 1/2 S/b = 0.75 (curve 5) values. Then 

we draw horizontal lines of the previously found valuee of 1/2 Cb/e = 
3.32 to inte~sect the curve 5, and 1/2 Co/e = 2.26 to intersect the 


curve 2. The points of intersection are projected onto the horizontal 
axis to a single point d,/b «© 0.49. In view of symmetry, cell 5 
has the same dimensions. 


Cell 2 has the coupling capacitances C, 1/€ = 1.92; Cy 2/e = 
) 3 


0.48, curve b. The auxiliary constructicn is mede similarly, and 
through the chord midpoints we draw the auxiliary curve da; through 
the point of intersection of the auxiliary curve d with the lines 
1/2 Ci/e # 3.17/2 = 1.585, a straight line is drawn parallel to the 


chords. The points of intersection of the parallel line with the 
curves a — Cy /& = 1.92, and b -— C, o/e = 0.48 have the coor- 
> > 


dinates (1.28; 0.07) and (1.88; 0.26). Hence the half-distances 
to the neighboring bars are: 


1/g Se, ,/o = 0,07, "/y Sis/b = 0,26. 


In order to determine the bar diameter d,/b; the points of 


intersection of the vertical lines 1/2 S/d = 0.07 and 1/2 S/b ® 0.26 
with the lines of constant d/r values (1 - 8) (Figure 14), are trans- 
ferred to Figure 15, and through them we draw the curves 1 and 3 
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corresponding to these values 1/2 S/b = 0.07 and 1/2 S/b = 0.26, and 
find the points of intersection of these curves with the horizontal 
lines 1/2 Ci/e = 1.88 and 1/2 Ci/e = 1.28, respectively. The found 


points, projected onto the horizontal axis, yield the sought value 
d,/b = 0.32. Cell 4 has similar dimensions. 


Cell 3 has the same coupling capacitances C, /é = C, 3/¢ = 0.48. 
3 2 


which simplifies the calculation: in Figure 14 we draw the horizontal 
line 1/2 Co/e = 4,13/2 = 2.065 to intersect the curve b (Cy 2/e = 
’ 


0.48). The horizontal projection of the point of intersection yields 
the sought value 1/2 S) o/b = 1/2, 85 3/> # 0.27. The points of in- 
3 a 


tersection of the vertical 1/2 S/b = 0.27 with the lines (1 - 8) of 
constant d/b values (see Figure 14) are transferred to Figure 15, 
and through them we draw the corresponding curve 4. The horizontal 
coordinate of the point of intersection of this curve with the hori- 
zontal line 1/2 Co/e = 2.065 yields the sought value d,/b = 0.36 


In determining the filter cross section dimensions, the relative dis- 
tance between the neighboring bars is obtained by summing like half- 
distances found for each of the bars. The basic dinerisions of the 
filter are summarized in the following table. 







+ = 0,36 + = 0,32 a = 0,40 


Se 15 SF 0.07 | 0.881 Fata zd =f 21,5 
i 


Note: Commas represent decimal points. 
Figure 16 shows the filter cross section. 


Construction of lumped loading condensers. The filter resonator 
loading condenser capacitance can be determined with the aid of the 


nomogram of Figure 9. In the frequency range in question, the 
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condenser capacitance for 
the various filters does 
not exceed 10 pF. Usually 
such a capacitance is pro- 
vided by a flat condenser 
with two plates, one of 
which is a metal capaci- 
tive plate mounted on 

the end of the resonator 
bar, and the other is a 
corresponding metal block 
which is in electrical Figure 16. Combeline filter (top 
contact with the upper view) 

and lower metal plates 
(Figure 17). Usually, 
these metal blocks are made 
movable. This makes it 
possible to regulate load- 
ing condenser capacitance 
in the tuning process, 
thereby compensating for 
any possible small devia- 








Figure 17. Cross section of comb- 
line filter with circular bars, cal- 
tions in the filter culated in the example 


dimensions. 


Figure 18 shows a comb filter loading condenser unit [2]. In 
this unit, tuning of the loading condenser is accomplished in two 
steps: rough tuning is first accomplished by displacing the metal 
block, after which the position of the block is fixed, and fine 
filter tuning 1s accomplished with the aid of trimming screws. In 
other words, this process can be termed two-step tuning. 


Figure 19 shows a comb filter loading condenser unit which 
permits accomplishing one-step tuning [9]. Here, a system of 
threaded joints, including a screw mounted on a block and a screw 
mounted on a special bracket (yoke), which in turn is rigidly 


183 


7115 











By J 
aca ora 





CZ LLLLLL LL a (_——$——— 


Figure 18. Resonator bar Figure 19. Resonator bar load- 
loading condenser unit with ing condenser unit with "single- 
"two=step" tuning: step" tuning: 
1 — trimming screws; 2 — 1 — trimmer block; 2 — sleeve; 
trimming block; 3 — capaci- 3 — screws; 4 — capacitive 
tive plate plate 


attached to the filter case, is 


a4 


sf ft 
used to displace the block with a ce an ae ae i 
small step and fix its position. The m 
screws are made movable relative to Z i 
one another by means of a4 metal UW 
sleeve with internal thread. Mobil- a 


ity of the entire unit is achieved by In/ Out In¢ 
making the screws with different 

pitch or different direction of the Figure 20. Electrical 
threading. As the sleeve is rotated, cirouit of filter for fre- 


quency summation with 
the metal block displaces, and the bandpass elements 


displacement step is equal to the sum 

of the thread pitches when the threads are made in opposite direc- 
tions, and equal to the difference of the pitches when the screws 
have threads in the same direction. 


Frequency separation using comb-line bandpass filters. Band- 
pass filters have found wide application as components of circuits 
for frequency division and summation. Figure 20 shows the electrical 
circuit of such a fiiter [10]. The filter has nine bars. Bars 2, 

3, 4 and 6, 7, 8 are bandpass element resonators, and bars 1, 5, 9 
are matching transformers. The element including resonators 2, 3, 4 
is desigied to pass a signal with frequency fii the element including 
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resonators 6, 7, 8 is deaigned to pass a signal with frequency fo: 


In this filter, the transformer 5 1s loaded on both sides. Here, 
it is important that the mutual coupling capacitances of transformer 


5 with bars 4 and 6 not very along its length, 1.e., that the parallel 


coupling segment lengths be the same on both sides. This implies 
that the transformer 5 and resonators 4 and 6 associated with it must 
have the same geometric length. 


As we have noted previously, the resonator bars in the comb-line 
filter may have electrical length 6 <¢7/4. This possibility permits 
making the bars of the subject separation filter of identical geo- 
metric length 2. The bar electrical length is determined by the 
geometric length and operating frequency: 


==, (14) 
= ly (15) 


where A = C/2nf is the wavelength. Since ty = Los the bar electrical 


lengths are connected with one another by the relation: 


aah, (16) 
In designing the filter, we first select the resonator bar 
electrical length of one of the bandpass elements, determine the bar 
geometric length, and use the above relation to determine the elec- 
trical length of the bars of the second element. The further calcue 
lation is analogous to that of conventional comb-line bandpass ele- 


ments. The results of an experiment with a separation filter designed 


using this technique for 8, ./0, = 0.8 showed good agreement between 


the characteristics obtained and those assumed in calculating the 
bandpass elements. 
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